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~<XtS*  ABSTRACT 


The  primary  purpose  of  the  Transonic  Turbine  Test  Rig  of  the  United 
States  Naval  Postgraduate  School,  Monterey,  California,  is  to  permit 
investigation  of  turbine  performance  and  flow  conditions  for  transonic 
rotor  inlet  velocities.  This  turbine  also  has  provisions  for  investi- 
gation of  the  effects  on  performance  of  widely  varied  axial  and  tip 
clearances.  High  turbine  total-to-total  pressure  ratios  are  achievable, 
since  an  ejector-type  exhauster  assembly  can  maintain  a  partial  vacuum 
in  an  airtight  turbine  discharge  hood.  This  study  is  concerned  primarily 
with  initial  turbine  performance  tests— of  the  nature  of  suitability 
tests— both  with  and  without  the  rotor  installed.  The  installation  is 
described,  and  results  of  flow  rate  calibration  and  labyrinth  seal 
leakage  tests  are  included,  along  with  turbine  performance  test  results. 
Al**j presented  are  a  simplified  three-dimensional  turbine  performance 
prediction  and  a  one-dimensional  analysis  of  exhauster  performance.  The 
turbine  performance  prediction  method  uses  the  so-called  equivalent 
enthalpy  (and  corresponding  equivalent  total  thermodynamic  properties) 
to  permit  analysis  of  the  relative  flow  through  the  rotor  in  a  manner 
analogous  to  that  used  for  absolute  flows. 


TABLE  OF  CONTENTS 


Section 

Page 

1. 

Introduction 

19 

2. 

Installation 

20 

3. 

Turbine  Performance  Prediction 

38 

4. 

Exhauster  Performance  Analysis 

79 

5. 

Determination  of  Flow  Rates 

96 

6. 

Force 

and  Moment  Measurement  Systems  Calibration 

114 

7. 

Turbine  Tests 

118 

8. 

Discus 

ssion  and  Recommendations 

145 

Bibliography- 

148 

Appendix 

I 

Formula  Development  for  Turbine  Performance 
Prediction 

149 

II 

Computation  of  Flow  Angles  and  Loss  Coefficients 

162 

III 

Computer  Program  for  Turbine  Performance 
Prediction 

176 

IV 

Computer  Program  for  Exhauster  Analysis 

202 

V 

Computer  Programs  Pertinent  to  Flow  Rate 

216 

Measurements 
VI   Computer  Program  for  Turbine  Test  Data  Reduction    236 


LIST  OF  TABLES 

TABLE  Page 

I.   Performance  Analysis  Input  Values  for  which  Results         64 
are  not  Available 

II.   Summary  of  Performance  Analysis  Contents  of  Figs.  19-30     67 

III.   Flow  Equation  Constants — Square-Edged  Orifices,  102 

Transonic  Turbine  Test  Rig 

IV.   Flow  Nozzle  Discharge  Coefficients,  Transonic  Turbine      107 
Test  Rig 

V,   Summary  of  Formulas  for  Determining  Flow  Rates,  113 

Transonic  Turbine  Test  Rig 

VI.   Percent  Tolerances  of  Flow  Equation  Variables  and         115 
Parameters — Square-Edged  Orifices,  Transonic  Turbine 
Test  Rig 

VII,   Summary  of  the  Turbine  Test  Results  Presented  in  135 

Figures  41  through  49 

VIII.   Calculation  of  Predicted  Stator  Gas  Outlet  Angles         167 

IX.   Calculation  of  Predicted  Rotor  Gas  Outlet  Angles  169 

X.   Calculation  of  Predicted  Rotor  Loss  Coefficients  172 


LIST  OF  ILLUSTRATIONS 

Figure  Page 

1.  Piping  Installation,  Transonic  Turbine  Test  Rig  21 

2.  Turbine  Hood  and  Exhauster  Assembly,  Transonic  Turbine  22 
Test  Rig 

3.  Flow  Measurement  Nozzle,  Transonic  Turbine  Test  Rig  24 

4.  Labyrinth  Seal  Configurations,  Transonic  Turbine  Test  Rig    25 

5.  Exhauster  Configuration,  Transonic  Turbine  Test  Rig  26 

6.  Turbine  and  Plenum  Assembly,  Transonic  Turbine  Test  Rig  28 

7.  Turbine,  Transonic  Turbine  Test  Rig  29 

8.  Sta tor-Plenum  Assembly,  Transonic  Turbine  Test  Rig  30 

9.  Assembly  Details,  Stator-Plenum  Assembly,  Transonic  31 
Turbine  Test  Rig 

10.  Gas  Path,  Transonic  Turbine  34 

11.  Stator  Blade  Profiles  35 

12.  Rotor  Blade  Profiles  36 

13.  Exhauster  Ejector  Head  37 

14.  Turbine  Expansion  Process  41 

15.  Coordinates,  Streamlines,  and  Angles  used  in  the  42 
Performance  Analysis 

16a.   Method  Used  for  Correcting  Streamline  Locations,  58 
Performance  Analysis 

16b.   Method  for  Approximating  Slopes— Performance  Analysis  60 

17.  Method  for  Determining  Rate  of  Change  of  Gradients  61 

18.  Schedule  of  Values  of  Input  Variables  used  in  the  66 
Performance  Analysis 

19.  Variation  of  Referred  Flow  Rate  with  Stator  Pressure  68 
Ratio—Performance  Analysis 

20.  Variation  of  Area  Restriction  Factor  with  Non-Dimensional    68 
Radius — Performance  Analysis 

21.  Variation  of  Stator  Exit  Mean  Streamline  Axial  Velocity  69 
with  Pressure  Ratio— Performance  Analysis 


Figure  Page 

22.  Radial  Distribution  of  the  Ratios  of  Stator  Exit  Axial      69 
Velocity  and  Tangential  Velocity  to  Mean  Streamline 

Axial  Velocity 

23.  Variation  of  Stator  Exit  Pressure  Ratio  with  Non-Di-        70 
mensional  Radius — Performance  Analysis 

24.  Variation  of  Increases  in  Non-Dimensional  Entropy  70 
with  Rotor  Exit  Non-Dimensional  Radius — Performance 

Analysis 

25.  Variation  of  Equivalent  Enthalpy  with  Rotor  Exit  Non-       71 
Dimensional  Radius — Performance  Analysis 

26.  Variation  of  Rotor  Exit  Axial  Velocity  Ratio  with  71 
Non-Dimensional  Radius — Performance  Analysis 

27.  Variation  of  Rotor  Exit  Referred  Axial  Velocity,  72 
Stator  Exit  Mean  Streamline  Mach  Number,  and  Head 
Coefficient  with  Referred  RPM — Performance  Analysis 

28.  Variation  of  Turbine  Efficiency  with  Head  Coefficient —     73 
Performance  Analysis 

29.  Variation  of  Theoretical  Degree  of  Reaction  with  74 
Head  Coefficient — Performance  Analysis 

30.  Variation  of  Referred  Rotor  Moment  with  Referred  RPM —      74 
Performance  Analysis 

31.  Possible  Exhauster  Flow  Processes  81 

32.  Variation  of  Turbine  Total  to  Total  Pressure  Ratio         94 
with  Inlet  Total  Pressure — Exhauster  Analysis 

3$.   Variation  of  Turbine  Inlet  Total  Temperature  Required       94 
to  Maintain  Specified  Outlet  Temperatures  with  Inlet 
Total  Pressure — -Exhauster  Analysis 

34.  Shaft  Seal  Flow  Function  106 

35.  Flow  Nozzle  Discharge  Coefficient  107 

36.  Plenum  Labyrinth  Seal  Discharge  Coefficient  109 

37.  Least  Squares  Parabola  for  the  Plenum  Labyrinth  Seal       110 
Flow  Function 

38.  Configuration  for  Calibration  of  Force  and  Moment         117 
Measurement  Systems 

39.  Illustration  of  the  Sta tor-Plenum  Assembly  for  Momentum    122 
Equation  Considerations 


Figure  Page 

40.  Radial  Equilibrium  of  a  Fluid  Particle  126 

41.  Experimentally  Determined  Values  of  Stator  Efficiency —    136 
TTTR 

42.  Experimentally  Determined  Values  of  Stator  Gas  Efflux      136 
Angle-- -TTTR 

43.  Experimentally  Determined  Values  of  Turbine  Efficiency —    137 
TTTR 

44.  Experimentally  Determined  Values  of  Referred  Flow  Rate —    138 
TTTR 

45.  Experimentally  Determined  Values  of  Theoretical  Degree     138 
of  Reaction— TTTR 

46.  Experimentally  Determined  Values  of  Referred  Rotor        139 
Moment — TTTR 

47.  Experimentally  Determined  Values  of  Rotor  Efficiency —     140 
TTTR 

48.  Experimentally  Determined  Values  of  Referred  Stator        141 
Axial  Velocity— TTTR 

49.  Experimentally  Determined  Values  of  Referred  Rotor        142 
Axial  Velocity— TTTR 

50.  Illustration  of  Boundary  Layer  Effects  at  the  Exit        156 
Plane  of  a  Blade  Row 

51.  Thermodynamic  Process  through  a  Blade  Row  156 

52.  Energy  Parameter,  H***,  as  a  Function  of  Pressure  Ratio     161 

53.  Blade  Geometry  used  in  Gas  Angle  Predictions  162 

54.  Variation  of  Predicted  Stator  Gas  Outlet  Angle  with        168 
Mach  Number 

55.  Variation  of  Predicted  Stator  Gas  Outlet  Angle  with        168 
Radius 

56.  Variation  of  Predicted  Rotor  Gas  Outlet  Angle  with        171 
Mach  Number 

57.  Variation  of  Predicted  Rotor  Gas  Outlet  Angle  with        171 
Radius 

58.  Variation  of  Rotor  Loss  and  Velocity  Coefficients  with     175 
the  Ratio  of  Incidence  Angle  to  Positive  Stalling 

Incidence  Angle 


TABLE  OF  SYMBOLS 

Symbols 

2 

A  Cross-sectional  area  (in.  ) 

A  Constant  dependent  on  square-edged  orifice  geometry  (dimension- 
less). 

a  Exit  opening  of  a  cascade  (in.). 

B  Constant  dependent  on  square-edged  orifice  geometry  (in.  ). 

B  Subscripted,  one  through  three;  see  Eqs.  (122)-(124). 

b,  Blade»s  departure  from  being  straight-backed  (in.).;  see  Fig.  53. 

b  Between-teeth  chamber  heighth,  labyrinth  seals  (in.). 

C  Constant  in  square-edged  orifice  equation  (("R/in.Hg-in.H^O)^  lbm/sec.) 

C-j^  Conversion  factor,  2gJ  (ft2-lbm/sec2-BTU). 

C2  Conversion  factor,  2gjc  (ft2/sec2-°R). 

c  Chord  of  blade  (in.). 

Cp  Specific  heat,  constant  pressure  (BTU/lbm-°R). 

D  Diameter  (in.). 

E  Kinetic  energy  (ft-lb). 

e  Curvature  of  the  back  of  a  blade;  see  Eq,  (288)  and  Fig.  53. 

F  Force  (lb). 

f  Correction  factor,  see  Eq.  (200). 

g  Universal  gravitational  constant  (32.174  lbm-f t/lb-sec  ) . 

H  Total  enthalpy  (BTU/lbm). 

H  Energy  parameter,  boundary  layer  (dimensionless);  see  App.  I. 

h  Static  enthalpy  (BTU/lbm). 

h  Blade  heighth  (in.). 

b^  Differential  pressure  (in.  fi^O). 

HP  Horsepower. 

I  Integrand. 
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Symbols 

i  Incidence  angle  (deg.  or  radians). 

is  Positive  stalling  incidence  angle  (deg.  or  radians). 

i  Unit  vector. 

J  Conversion  factor  (778.16  ft-lb/BTU). 

j  Dimension  for  determining  blade-back  curvature  (in.);  see  Fig.  53. 

K  Discharge  coefficient,  flow  metering  or  labyrinth  seal  (dimension- 
less). 

k  Tip  clearance  (in.). 

k-  Head  coefficient  (dimensionless);  see  Eq.  (  81). 

L  Distance  between  stations  ahead  of  and  after  a  rotor  (in.). 

M  Mach  Number  (dimensionless). 

M  Moment  (ft-lb  or  in-lb). 

m  Mass  flow  rate  (slugs/sec). 

m  Exponent  used  in  boundary  layer  calculations  (dimensionless); 
see  p.  45  . 

N  Rotational  speed  (RPM). 

n  Number  of  throttlings,  labyrinth  seals  (dimensionless). 

n  Polytropic  exponent  (dimensionless);  see  pp.  79-80, 

n  Outward  normal  (unit) vector. 

P  Pressure  (in.  Hg.  abs). 

p  Pressure  (psia). 

R  Gas  constant  for  air  (53.345  ft-lb/lbm-  R). 

R  Reaction  force  (lb). 

Re  Reynolds  Number  (dimensionless). 

r  Labyrinth  seal  overall  (inlet  total  to  exit  static)  pressure 
ratio  (dimensionless). 

r  Radius  (in.). 

r*  Theoretical  degree  of  reaction  (dimensionless);  see  Eq.  (80). 
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Symbols 

s  Entropy  (BTU/lbm-°R) . 

s  Non-dimensional  entropy  (s/c  )  or  entropy  increase  (  As/c  ). 

s  Cascade  blade  pitch  (in.). 

s  Surface  area  (in.  ). 

s  Labyrinth  seal  between-teeth  spacing  (in.). 

T  Temperature  (PR). 

t  Temperature  (°F). 

t  Blade  thickness  (in.). 

t  Labyrinth  seal  tooth  thickness  (in.). 

te  Blade  trailing  edge  thickness  (in.). 

t*    Projection  of  Blade  trailing  edge  thickness  on  the  exit  plane  of 
the  blade  row  (in.). 

U  Peripheral  velocity  (ft/sec). 

u  Velocity  within  a  boundary  layer  (ft/sec). 

V  Absolute  velocity  (ft/sec). 
W  Relative  velocity  (ft/sec). 
W  Flow  rate  (lbm/hr). 

w     Flow  rate  (lbm/sec). 

Wf    Fraction  of  the  total  flow  rate  which  passes  between  the  hub 
streamline  and  any  other  streamline  (dimensionless). 

\frf       Reference  flow  rate  parameter  (in.2)j  see  Eq.  (16). 

X     Non-dimensional  radius,  r/rm,  where  rm  is  the  mean  streamline 
radius  (dimensionless). 

X     In  square-edged  orifice  equations,  Re(10~°),  (dimensionless). 

X     In  exhauster  performance  analysis,  the  axial  distance  down- 
stream of  the  inlet  to  the  diverging  portion  of  the  exhauster 
piping  (in.). 

Xe    See  Eq.  (268),  (dimensionless). 

Y  Non-dimensional  axial  velocity  ratio,  vaAam>  where  V-  is  the 
mean  streamline  axial  velocity.  ■ 
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Symbols 

Y^    In  flow  measurement  equations,  expansion  factor  (dimensionless); 
see  Eq.  (159). 

Y  Pressure  loss  coefficient  used  in  determining  the  loss  coefficient, 

Z  Number  of  blades  (dimensionless). 

zf  Absolute  viscosity  in  centipoises. 

z  =  lOOz  (centipoises). 

z  Axial  distance  in  cylindrical  coordinates  (in.). 

Greek  Letters 

0<    Discbarge  coefficient,  labyrinth  seals  (dimensionless). 
(X         Stator  gas  flow  angle  (deg,  or  radians) 

ex"    Coefficient  of  thermal  expansion,  flow  metering  (dimensionless). 
fl         Rotor  gas  flow  angle  (deg.  or  radians). 
/}         Rotor  geometric  blfcde  angle. 


A 


Ratio  of  orifice  or  nozzle  minimum  diameter  to  upstream  pipe 
diameter,  in  flow  measurement  (dimensionless). 


T    Specific  heat  ratio,  c  /cv  (dimensionless). 

ST 

T*  Velocity  carryover  factor,  labyrinth  seals  (dimensionless). 

%  Labyrinth  seal  tooth  radial  clearance  (in.). 

8  Boundary  layer  thickness  (in.). 

6  Referred  pressure  pt  /p   ,,  Pstd  =  14.7psia  (dimensionless). 

S        Boundary  layer  displacement  thickness  (dimensionless). 

Boundary  layer  energy  thickness 
£    Less  coefficient  (dimensionless). 

f    Correction  factor  for  Reynolds  Number,  square-edged  orifices 
(dimensionless) • 

i{  Efficiency  (dimensionless). 

71  Non-dimensional  distance  from  the  wall  in  a  boundary  layer  (y/6  )• 

0  Referred  temperature,  Tt  A  td,  Tgtd=518.4  R  (dimensionless). 

K  Curvature  factor  (dimensionless);  see  Eq.  (258). 
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Greek  Letters 

X    Angle  of  flow  in  a  meridional  plane  (deg.  or  radian). 

X     Factor  used  in  estimating  blading  secondary  loss  coefficients 
(dimensionless);  see  App.  II, 

\«;x  Pressure  loss  coefficient,  blading  leaving  loss  (dimensionless). 

f  Area  restriction  factor  (dimensionless);  see  Eqs.  (21)  and  (187). 

/o  Density  (slugs/ft3). 

$  Non-dimensional  flow  function;  see  Eq.  (10). 

(p    Labyrinth  seal  pressure  ratio  function  (dimensionless);  see 
Eq.  (164). 

(Py    Absolute  flow  velocity  coefficient,  V/Vjjj,  (dimensionless). 

"Ij/    Relative  flow  velocity  coefficient,  W/W<j>H,  (dimensionless). 

or   Angular  velocity  (radians/sec). 

Subscripts 

a     Ambient  conditions. 

a     Axial . 

ACT   Actual  (or  computed)  value. 

b     Refers  to  flow  through  stator  blade  row,  exhauster  performance 
analysis. 

c     Subsonic  compression,  in  exhauster  performance  analysis. 

c     In  turbine  equations,  refers  to  the  flow  director  ring  shrouding 
the  rotor. 

d  In  exhauster  performance  analysis,  turbine  hood  conditions. 

E  Refers  to  so-called  equivalent  thermodynamic  properties. 

e  In  exhauster  performance  analysis,  supersonic  expansion. 

e  In  turbine  equations,  entrance  station  of  plenum. 

ef f  Effective . 

f  Due  to  or  corrected  for  friction. 

H  Hub. 
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Subscripts 

I     In  the  vicinity  of  the  cover  plate  (see  Fig.  6). 

i     Used  with  efficiency,  1\  ,  to  designate  overall  turbine  (internal) 
efficiency. 

i     Used  with  M  to  denote  the  moment  about  the  turbine  centerline 
acting  on  the  stator-plenum  assembly. 

is  Value  for  isentropic  expansion  from  inlet  conditions. 

ft  Refers  to  the  turbine  plenum  labyrinth  seal. 

jl  In  exhauster  analysis,  downstream  of  a  shock. 

m  Mean  streamline. 

N  In  exhauster  analysis,  refers  to  ejector  nozzle  flow. 

n  Refers  to  flow  measurement  nozzle. 

o  Denotes  inlet  conditions. 

o     Used  with  the  coefficient  of  thermal  expansion  to  denote  the 
square-edged  orifice  value. 

p  Used  with  force,  F,  to  denote  pressure  forces. 

p  Plenum  condition. 

p  Used  with  efficiency,  y\   ,  to  denote  polytropic  efficiency. 

p  Used  with!  loss  coefficient,  £  ,  to  denote  profile  loss  coefficient, 

R  Relative  flow  value. 

R  Rotor  value. 

r  Radial . 

ref  Referred  value. 

req  Required. 

s  Stator. 

s  In  flow  rate  equations,  refers  to  shaft  labyrinth  seal. 

s     As  second  subscript  denotes  value  for  isentropic  expansion  from 
equivalent  total  conditions. 

s     Used  with  loss  coefficient  to  denote  the  secondary  loss  coef- 
ficient. 
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Subscripts 

T  Tip. 

T  In  exhauster  performance  analysis,  turbine. 

TH  Theoretical  value,  i.e.,  for  isentropic  expansion, 

u  Tangential . 

u  In  exhauster  performance  analysis,  upstream  of  a  shock. 

z  Refers  to  the  axial  direction,  cylindrical  coordinates* 

8  Refers  to  angle  coordinate  in  cylindrical  coordinates 

1  Station  between  the  stator  and  the  rotor. 

1  Station  at  exhauster  nozzle  exit  plane. 

1  Streamline  number  one  (hub). 

2  Station  at  the  rotor  exit. 

2  Station  at  the  inlet  to  the  diverging  portion  of  the  exhauster. 

2  Streamline  number  two  (between  hub  and  mean  streamlines). 

3  Station  at  the  exit  of  the  exhauster  piping. 

3  Streamline  number  three  (mean). 

4  Streamline  number  four  (between  mean  and  tip  streamlines). 

5  Streamline  number  five  (tip). 

Superscripts 

*  Denotes  choke d-f low. 
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1.  Introduction, 

Little  qualitative  information  is  available  concerning  turbine 
performance  under  conditions  of  transonic  rotor  entrance  velocities  and 
extreme  total  to  static  turbine  pressure  ratios.  Additionally,  axial  and 
radial  clearances  in  turbines  for  use  in  rocket  applications  and  gas 
turbine  plants  are  often  subject  to  wide  variations  because  of  extreme 
temperature  variation. 

Information  regarding  radial  clearance  effects  on  efficiency  avail- 
able to  designers  indicates  that  tip  clearances  should  be  limited  to  about 
one  percent  of  blade  heighth.  However,  the  smaller  a  turbine,  the  more 
expensive  it  becomes  to  achieve  the  close  tolerances  required  to  maintain 
this  percent  tip  clearance.  Also,  temperature  effects  become  more  critical 
for  smaller  turbines,  with  respect  to  maintaining  the  desired  radial 
clearance  over  a  range  of  operating  conditions. 

Much  less  information  is  available  to  designers  on  the  effects  of 
axial  clearances  on  performance.  Methods  of  analysis  of  turbine  perform- 
ance generally  used  in  turbine  design  assume  axisymmetric,  uniform  flow 
conditions  between  blade  rows,  but  do  not  consider  the  axial  distance 
required  to  permit  the  flow  to  achieve  such  stabilized  conditions. 
Consequently,  axial  clearances  between  blade  rows  are  dictated  primarily 
by  space  and  weight  considerations,  rather  than  by  knowledgeable  consider- 
ations of  axial  clearance  effects.  It  is  even  conceivable  that  optimum 
axial  clearances,  from  a  performance  viewpoint,  may  be  less  than  normally 
used  for  some  types  of  blading. 

The  Transonic  Turbine  Test  Rig  installed  at  the  Propulsion  Laboratory 
of  the  Department  of  Aeronautics,  U.S.  Naval  Postgraduate  School,  Monterey, 
California,  was  designed  by  Professor  M.  H.  Vavra  of  that  department  to 
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permit  investigation  of  the  above-stated  problems.  This  thesis  is 
concerned  with  describing  the  installation  and  its  capabilities,  calibration 
tests,  and  initial  performance  tests  of  the  turbine  rig,  with  a  stator 
consisting  of  converging  nozzle  type  blading  and  an  impulse  type  blading 
rotor. 

Funds  for  designing  and  building  the  Transonic  Turbine  Test  Rig,  as 
well  as  operating  funds,  were  furnished  by  the  Bureau  of  Naval  Weapons 
(RAPP-14).  The  generosity  of  that  Bureau  is  greatly  appreciated.  The 
author  is  most  appreciative,  too,  of  the  guidance  and  counsel  provided 
by  Professor  Vavra,  and  of  the  assistance! in  performing  experiments  so 
willingly  given  by  Mr.  L.  T.  Clark  of  the  Department  of  Aeronautics. 
2.  Installation. 

The  Transonic  Turbine  Test  Rig  (TTTR)  installation  is  shown  schemati- 
cally in  Fig.  1,  whereas  Fig.  2  is  a  photograph  of  the  turbine  hood  and 
exhauster  assembly  (with  the  exhauster  rotated  180°  from  its  actual 
configuration  in  the  TTTR  installation) .  The  turbine  proper  is  a  single- 
stage,  impulse-type  design.  The  salient  components  of  the  installation 
shown  in  Fig.  1  are  thes 

1.  | flow  measurement  nozzle; 

2.  removable  standard  orifice  installation; 

3.  plenum  labyrinth  leak  rate  measurement  standard  orifice  (two- 
inch  pipe  orifice,  Fig.  1); 

4.  turbine  plenum  labyrinth  seals; 

5.  shaft  labyrinth  seals; 

6.  turbine  hood;  and  the 

7.  turbine  proper  and  turbine  plenum. 

This  turbine  can  be  operated  exhausting  either  to  atmospheric 
pressure  or,  with  the  hood  and  exhauster  installed,  to  a  partial  vacuum. 
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Notes: 
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tios greater  than  $tl. 
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Supply  air  is  furnished  by  the  Allis  Chalmers,  VA  312  compressor  installed 
at  the  Propulsion  Laboratory.  Supply  air  temperature  is  regulated  by 
means  of  an  aftercooler  installed  downstream  of  the  compressor.  The 
maximum  temperature  of  air  supplied  by  the  compressor  is  300  F.  However, 
since  the  turbine  has  aluminum  blading,  a  limiting  value  of  200° F  has 
been  established  for  turbine  inlet  total  temperature.  The  maximum 
turbine  inlet  total  pressure  attainable  is  about  43  psia.  Exhausting  to 
the  atmosphere,  the  total  to  static  pressure  ratio  possible  is  about  2.8, 
and  a  maximum  total  to  total  pressure  ratio  of  about  5.6  is  anticipated 
when  exhausting  to  a  partial  vacuum  (see  Section  4,  Exhauster  Performance 
Analysis).    Maximum  turbine  RPM  is  20,000  and  maximum  horsepower  is  200. 

The  flow  rate  through  the  turbine  blading  during  operation  is  the 
flow  measured  by  the  flow  nozzle  less  the  plenum  labyrinth  seal  leakage. 
The  flow  nozzle,  with  pressure  and  temperature  tap  locations  indicated, 
is  shown  in  Fig.  3,  and  the  plenum  labyrinth  seal  configuration  is  shown 
in  Fig.  4.  Since  there  is  insufficient  piping  upstream  of  the  flow  noz- 
zle to  ensure  undisturbed  flow  through  the  nozzle,  calibration  of  the  flow 
nozzle  was  required.  The  removable  standard  orifice  assembly  is  provided 
for  that  purpose.  Plenum  labyrinth  seal  leak  rate  can  not  be  measured 
during  turbine  operation j  hence,  the  two-inch  pipe  orifice  is  provided 
for  determining  this  flow  rate  in  separate  tests.  During  such  tests,  all 
flow  to  the  plenum  passes  through  the  two-inch  pipe  and  the  only  exit  for 
the  flow  is  through  the  plenum  seals,  since  the  stator  exit  is  blocked  by 
a  special  cover  plate. 

In  analyzing  exhauster  performance,  the  flow  rate  to  the  exhauster 
ducting  (see  Fig.  5)  consists  of  the  flow  from  the  turbine  hood  and  the 
flow  through  the  exhauster  nozzle.  The  former  consists  of  the  flows 
through  the  turbine  blades,  through  the  plenum  labyrinth  seals,  and  through 
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the  shaft  labyrinth  seals  (from  the  atmosphere  to  the  partial  vacuum  in 
the  turbine  hood).  There  are  no  provisions  for  measuring  the  leakage 
through  the  shaft  seal,  and  this  flow  rate  must  be  estimated.  Also,  flow 
through  the  exhauster  nozzle  can  not  be  measured  directly.  However,  with 
the  entire  flow  of  the  compressor  routed  to  the  turbine  rig,  the  flow 
through  the  exhauster  nozzle  is  the  flow  rate  delivered  by  the  compressor, 
determined  as  specified  by  Vavra,  [9]     minus  the  flow  rate  measured  by  the 
flow  measurement  nozzle.  Additionally,  since  the  flow  through  the  exhauster 
nozzle  will  normally  be  choked,  the  flow  rate  can  be  determined  quite 
accurately  by  measuring  the  total  conditions  ahead  of  the  nozzle. 

Both  the  removable  and  the  two-inch  pipe  orifice  assemblies  conform 
to  specifications  for  standard  square-edged  orifice  installations  as  set 
forth  by  Stearnes,  et.  al,  [2]   Detailed  information  regarding  specifi- 
cations applicable  to  these  assemblies,  including  pertinent  dimensional 
information,  is  given  by  Eckert.  [6] 

The  turbine  proper  and  the  turbine  plenum  are  shown  schematically  in 
Fig.  6.  The  rear  support  and  front  hanger  shown  are  attached  in  turn  to 
a  cradle  flexure  assembly  affixed  to  the  rig  table.  The  set  screw  at  the 
rear  support  is  provided  for  setting  axial  clearances.  Figs.  7  and  8  are 
photographs  of  the  turbine  with  and  without  the  rotor,  respectively.  Fig. 
9  shows  the  manner  in  which  the  floating  armature  assembly  fits  into  the 
plenum.  As  shown  in  Fig.  6,  the  following  instrumentation  is  installed 
in  the  turbine  and  plenum  assemblys 

1.  Static  pressure  (p  )  and  static  temperature  (T  )  in  the  plenum; 

2.  Total  pressure  (pt  ),  static  pressure  (p0)  and  total  temperature 
(Tt  )  ahead  of  the  statorj 

^Figures  in  brackets  refer  to  bibliography  entries  of  Section  9. 
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3.  Static  pressures  at  the  hub  (piH)  and  at  the  tip  (piT)  of  the 
stator  in  the  exit  plane;  and 

4.  Reluctance  gages  for  measuring  the  axial  force  and  the  moment 
about  the  axis,  acting  on  the  stator-plenum  assembly. 
The  air  dynamometer  used  for  measuring  rotor  torque  is  shown  in  Fig.  7. 

Spoked-wheel  type  flexures  are  installed  at  both  the  front  and  rear 
bearings  of  the  rotor  bearing  housing  for  measuring  the  axial  force  on 
the  rotor.  Strain  gages  are  mounted  on  both  sides  of  each  of  the  three 
webs  of  these  flexures.  Pressure  taps  are  installed  around  the  sides  and 
in  the  top  of  the  turbine  hood  for  measuring  turbine  total  discharge  pres- 
sure (pt(j)  during  operation  with  the  hood  and  exhauster  installed. 

The  thrust  and  torque  flexure  shown  in  Fig.  6  was  included  to  permit 
determination  of  the  pressure  (p^T)  acting  on  the  cover  plate  concentric 

with  and  inside  the  stator  ring.  The  axial  force  acting  on  the  cover 

U  1  0      & 

vt    plate  will  be  determined  by  a  flexure  in  the  position  shown.  The  pressure 

p  (PlT)  can  tnen  be  determined  from  this  force  and  the  pressure  in  the 

cavity  behind  the  cover  plate  (pih).  However,  the  flexure  shown  has  not 
been  instrumented,  since  stator-only  tests  showed  it  to  be  too  critical 
in  wall  thickness  and  axial  groove  size,  with  respect  to  desired  sensiti- 
vity and  adequate  strength.  In  the  future,  this  set-up  will  be  replaced 
by  a  spoked-wheel  type  flexure. 

The  instrumentation  of  the  stator  and  plenum  assembly  was  designed 
to  permit  using  the  momentum  and  moment  of  momentum  equations  to  deter- 
mine stator  performance.  For  accuracy  in  such  performance  measurement, 
the  pressures  at  the  hub  and  at  the  tip  of  the  stator  exit,  and  the  pres- 
sure acting  on  the  cover  plate  must  be  known,  along  with  the  reaction 
force  and  moment  acting  on  the  assembly.  Additionally,  the  average  of 
the  hub  and  tip  pressures  is  a  better  value  for  stator  exit  pressure  (p-^) 
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than  either  value  by  itself.  During  stator  tests  without  the  rotor,  the 
pressure  acting  on  the  cover  plate  was  measured  directly  by  using  taps 
installed  inside  a  specially  fabricated  cylinder  placed  over  the  rotor 
bearing  housing.  However,  such  pressure  could  not  be  determined  for  the 
turbine  tests  conducted  to  date. 

The  stator-plenum  assembly  can  be  moved  axially  a  distance  of  1.2  in. 
to  achieve  various  stator  to  rotor  axial  clearances.  Moreover,  the  rotor 
can  be  moved  axially  approximately  0.25  in.  for  increased  range  and  flexi- 
bility in  setting  axial  clearances. 

Fig.  10  shows  the  turbine  gas  path  and  pertinent  dimensions  for  maxi- 
mum rotor  blade  heighth.  The  basic  configuration  shown  is  that  for  use 
with  the  converging  nozzle  stator  guide  vanes.  The  dashed  lines  indicate 
the  gas  path  when  the  second  set  of  guide  vanes,  of  converging-diverging 
nozzle  type,  are  installed.  It  is  evident  from  Fig.  10  that  a  range  of 
tip  clearances  from  the  minimum  practicable  to  about  0.6  in.  are  attain- 
able merely  by  reducing  the  rotor  outer  diameter.  Additionally,  through 
the  use  of  inserts  between  the  wheel  outer  diameter  and  the  case,  various 
gas  path  flare  angles  can  be  achieved,  along  with  varied  radial  clearances 
and  radial  overlaps  at  the  tip,  obtainable  by  varying  axial  clearance. 

The  blading  profiles  of  the  two  available  stators  are  shown  in  Fig. 
11,  and  the  profiles  of  the  two  available  rotors  are  shown  in  Fig.  12. 
The  area  ratio  —  exit  opening  to  throat  —  of  the  converging-diverging 
blading  is  1.25,  corresponding  to  an  isentropic  Mach  Number  of  1.6,  One 
rotor  blading  is  double-circular-arc  profile,  while  the  other  is  a  profile 
with  gradually  changing  curvature. 

The  exhauster  head  (nozzle  plus  turbine  hood  flow  inlet  section)  of 
the  exhauster  assembly  shown  in  Fig.  5  is  shown  in  Fig.  13.  Axial 
locations  of  pressure  taps  in  the  ducting  are  indicated  in  Fig.  5.  At 
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each  location,  four  taps  are  equally  spaced  circumferentially  and  commoned 
for  obtaining  an  average  reading.  Fig.  13  shows  also  the  locations  of 
pressure  taps  along  the  outer  contour  of  the  exhauster  head;  at  each 
location,  three  taps  are  equally  spaced  circumferentially  and  commoned. 
The  location  of  the  total  pressure  and  total  temperature  tap  upstream  of 
the  nozzle  has  not  yet  been  determined.  Consequently,  the  location  shown 
is  approximate  only. 

All  pressures  are  measured  in  inches  of  mercury  against  atmospheric 
pressure  as  reference,  using  96  in,  manometers  with  0.1  in.  graduations. 
Differential  pressures  across  the  flow  nozzle  and  orifice  taps  are 
measured  in  inches  of  water.  Temperatures  are  measured  using  an  ice  bath 
as  reference. 
3.  Turbine  Performance  Prediction. 

General 

The  basic  method  used  for  turbine  performance  prediction  is  that  given 
by  Vavra.   This  method  is  a  simplified  three-dimensional  analysis  based 
on  the  equation  of  motion  and  the  equation  of  continuity,  with  solutions 
obtained  at  between-blade-row  locations.  In  contrast  to  the  approach  of 
Vavra,  the  equations  of  motion  and  continuity  for  relative  flows,  rather 
than  for  absolute  flows,  were  used  in  rotor  calculations.  This  procedure 
simplifies  computations,  since  the  predicted  relative  rotor  gas  exit 
angles  can  be  used  directly  in  the  equation  of  motion.  Therefore,  itera-. 
tion  based  on  change  in  total  enthalpy  across  the  rotor  is  not  required. 
Also,  loss  coefficients,  estimated  using  the  method  of  Ainley,  [4,5]  are 
obtained  in  the  same  manner  for  the  rotor  as  for  the  stator.  Consequently, 
their  application  must  be  in  a  manner  analogous  to  that  used  for  the  stator. 

^Vavra,  M.  H.,  Aero-Thermodynamics  and  Flow  in  Turbomachines  (New 
York,  Londons  John  Wiley  and  Sons,  Inc.,  1960)  Chap.  16. 
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Proper  application  is  achieved  by  introduction  of  the  so-called  equivalent 
enthalpy  (defined  later).  The  assumptions  made  in  the  performance  analysis 
are  set  forth  below. 

1.  An  infinite  number  of  blades  is  assumed  so  that  downstream 
conditions  do  not  influence  conditions  upstream, 

2.  At  the  stations  of  interest,  axisymmetric  stream  surfaces  exist, 
so  that  all  derivatives  with  respect  to  0,  i,  e.,  in  the  peripheral 
direction,  are  zero. 

3.  The  flow  is  adiabatic  and  steady  so  that  along  any  given  stream- 
line there  is  no  change  of  total  enthalpy  in  the  stator  and  no  change  in 
relative  total  enthalpy  in  the  rotor.  (Note,  however,  that  since  cylin- 
drical coordinates  are  to  be  used,  gradients  of  total  and  relative  total 
enthalpy  can  exist  in  both  the  radial  direction  and  the  axial  direction, 
due  to  streamline  curvature,  depending  on  inlet  conditions.) 

4.  Since  the  equations  are  to  be  solved  for  between-blade-row 
locations,  the  increase  in  entropy  due  to  friction  within  the  blade  rows 
can  be  determined  using  the  loss  coefficient  and  the  inlet  and  outlet 
conditions.  However,  it  is  necessary  to  assume  that  entropy  changes 
along  the  streamlines  between  blade  rows  are  negligible. 

Development  of  Formulas 
The  basic  equations  of  motion  for  absolute  and  relative  flows  are: 
7H- 7x  (7x7)  +  TVs  (1) 

VHR  =  Wx(Vx#+Z£r)+T\7s  (2) 

Total  relative  enthalpy,  HR,  can  be  written: 

HR-K+f-^=K+Jf  -£         U-wru        (34) 
H,  =  tv+  4^  -  -£  Ob) 
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The  subscripts  indicates  1  —  ahead  of  the  rotor;  2  —  after  the  rotor; 
TH  —  the  theoretical  value;  and  2*  —  the  conditions  after  the  rotor  for 
isentropic  expansion  from  TtE»He/cp  t0  P2*  For  ac*iabatic  flow,  the 
total  relative  enthalpy  remains  constant  along  any  given  streamline.  Hence, 
equivalent  enthalpy,  defined  as: 

HE=^+  ^  =HR4  -£  (4) 

must  also  remain  constant  along  each  streamline.  It  follows  from  (3)  and 
(4)  that?  + 

bi  -  n2+-  -f-   -  h,  ■*-  -f g (5) 

It  is  evident  that  the  use  of  equivalent  enthalpy  permits  analysis  of  the 
relative  flow  through  the  rotor  in  a  manner  analogous  to  that  used  for 
absolute  flows,  including  the  use  of  normal  thermodynamic  relations. 
Velocity  coefficients  are  defined  as: 

<fc=Hb 

Then,  in  terms  of  the  loss  coefficients: 


for  absolute 

flows 

(6) 

for  relative 

flows 

(7) 

for  absolute 

flows 

(8) 

for  relative 

flows 

(9) 

Fig.  14  shows  the  relations  between  thermodynamic  properties,  velocities 
and  loss  coefficients  in  the  form  of  a  typical  expansion  through  a  turbine, 
The  expansion  process  is  shown  for  three  different  streamlines  in  order 
to  emphasize  the  variation  of  properties  which  have  to  be  considered  in 
the  analysis. 

Figs.  15a  -  15  c  show  the  coordinate  system  and  sign  conventions 
used.  Station  subscripts,  zero,  one  and  two  refer  to  ahead  of  the  stator, 
between  the  stator  and  rotor,  and  after  the  rotor,  respectively. 

The  equation  of  continuity  is  used  in  its  non-dimensional  flow  func- 
tion —  (J)  —  form.  Written  in  differential  form: 
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&fn = ^ = dA/  ^  (tf  - 


(10) 


(11) 

for  relative  flows.  Thens 

dA-fZadr  (12) 

Here,  fa»  is  the  blade  exit  opening  (see  Appendix  II)  and  the  factor,  f  , 
is  an  area  restriction  coefficient  which  accounts  for  loss  of  effective 
area  due  to  the  boundary  layer  on  both  sides  of  the  flow  channel.  This 
factor  is  further  defined  below.  Eq.  (12)  can  be  made  non-dimensional  by 
multiplying  and  dividing  by  amr ,  the  product  of  mean  streamline  exit 
opening  and  radius.  Hence  % 

d/\=  I  2  -§-dXamrm  (13) 

Eqa,(10)  and  (11)  in  integral  form  become,  respectively? 
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Stator 


For  wr£,  the  reference  flow  rate,  in  square  inches,  'a1  and  rm  must  be  in 
inches  and  pt  in  psia..  It  follows  from  (14)  and  (15)  that  continuity 


is  satisfied,  provided  thats 
,*r 


a-f$dX 
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J Rotor 


(17) 


Rotor 


Eq.  (17)  neglects  the  flow  area  available  because  of  tip  clearance.  The 
differential  area  due  to  tip  clearance  is  J 
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dA  =  Zi\r  dr 


(18) 


Letting  rc  denote  the  inner  radius  of  the  flow  direction  ring  around  the 


rotor,  the  integral  form  of  (15)  becomes: 
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(19) 


The  2  in  the  denominator  of  the  integrand  is  required  because  the  entire 
annular  area  rather  than  a  between-the-blade  area  results  from  integration 
of  (19).  For  small  tip  clearances,  sufficient  accuracy  should  result 
from  taking  X  of  the  integrand  as  a  constant  equal  to  Xj.,  Xc,  or  the 
average  of  the  two.  Considering  tip  clearance  area,  and  assuming  that 
the  values  of  5),  ^  ,  Pte,  and  Tt„  in  the  clearance  area  are  approximately 
the  same  as  those  for  the  rotor  tip,  (17)  becomes: 


^  J  Rotor 


(20) 


*r£ 


In  addition  to  the  requirement  that  overall  continuity  be  satisfied  by 

the  stator  and  rotor  flow,  it  is  necessary  that  the  flow  rate  between  each 

two  adjacent  streamlines  remains  constant. 
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The  term  f  in  (12)  is 


£   a 


H*"- 


H' 


(21) 


-1  +Sp 

Here,  gp  is  the  profile  loss  coefficient.  For  the  stator,  using  £p  in 
(21)  corresponds  quite  closely  to  the  often-used  practice  of  assuming  that 
one-half  the  total  loss  occurs  to  the  throat.  However,  for  the  rotor, 


Vavra,  M.  H.,  Problems  of  Fluid  Mechanics  in  Radial  Turbomachines 
( Rhode-Saint-Gene se,  Belgium?  Von  Karman  Institute  for  Fluid  Dynamics,  1965) 
VKI  Course  Note  55b,  pp.  G46-50. 
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there  is  no  such  direct  correlation.  In  (21),  jr*"*"  is  the  energy  para- 
meter defined  as  the  energy  thickness,  £***,  divided  by  the  displacement 
thickness, £  .  These  thicknesses  are  given  by: 


^-'H'-WfTjJk,*,] 


S*=  SJi-0- 


n 


U-Xtn  ) 


dn] 


where : 


*«='-(i 


r-i 
r 


(22) 
(23) 

(24) 


In  (22)  and  (23),  rj  ,  is  the  normally-used  non-dimensional  distance  from 
the  wall,  y/s    ;  also,  m  ■  l/n,  where  n  is  the  velocity  profile  exponent 
in:  (u/U)  =  (y/g  )Vn. 

The  development  of  (1)  and  (2)  into  usable  form  is  quite  lengthy. 
Consequently,  this  development  is  included  in  Appendix  I.  It  is  evident 
that  these  two  equations  differ  only  by  the  term  of  (2),  (Wx2tiir),  which 
represents  the  Coriolis  acceleration.  Consequently,  the  derivation  given 
in  App.  I  is  for  (2).  The  form  of  (1)  used  is  obtained  from  the  equation 
developed  from  (2)  by  substituting?  H  for  HE;  OCj  for  /&%   ;  and  U  =  0.  Also 
included  in  App.  I  are  the  derivations  of  (21)  -  (23),  along  with  a  plot 
of  H***  for  m  =  0.15.  The  final  form  of  (1)  is: 


—  ^  <•:  *,  ,  c,  cos1*,    dH 


c,  H  cos  *<*i 


iT/i' 


Y*Vl 


—  Sin  a<[ 


Ids* 
TJX, 


(25) 


and  that  of  (2)  is: 


^«v4f^-ts'^ 


ds? 


dt   (26) 
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Here; 


I  -  Va/Vam  (or  Wa/Wam) 

X  =  r/rm 

k  -   4  to  6,  streamline  curvature  factor 

8r  =  streamline  displacement 

L  ■  length  over  which  streamline  displacement  occurs 

s*  -  s/cp,  non-dimensional  entropy 

U  =  wheel  speed,  ft/sec. 

Ci  =  2gJ  -  (2) (32. 174) (778. 16)  ft2lbm/sec2BTU,  converts  H 
to  ft2/sec2  from  BTU/lbm. 
In  the  curvature  term  (±  f<rm  -^r) ,   the  plus  sign  is  for  curvature  whose 
radius  of  curvature  points  in  the  positive  radial  direction  (see  Fig.  15a, 
p.  42). 

Total  temperature  can  be  used  in  (25)  and  (26)  by  substituting  C2Tt 
for  C^H,  where:  C2  =  2gjcp.  The  entropy  term  in  these  equations  is 
evaluated  along  the  line  of  constant  pressure  for  the  blade  exit.  The 
required  formula  is:4 


S*  = 


s,-sc 


h 


• 

H  - 

vz       1 

ZgJ 

H  - 

(27) 


Differentiating  with  respect  to  the  non-dimensional  radial  coordinate: 
ds*_  d 


dX   dX 


in 


f  C 


C,  M  coslc< 


Y2  VC 


(28) 


C.Ha»fe(l-f) 

Since  Y2  appears  in  terms  on  the  right  side  of  (25)  and  (26),  these 

equations  must  be  solved  by  successive  approximations.  For  a  solution  to 
be  valid,  it  must  satisfy  both  overall  continuity  and  between-streamline 
continuity,  as  previously  mentioned.  Overall  continuity  is  controlled 

Vavra,  M.  H.,  Aero-Thermodynamics  and  Flow  in  Turbomachines  (New 
York,  London:  John  Wiley  and  Sons,  Inc.,  I960)  pp.  445-447. 
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primarily  by  the  axial  velocity,  while  between-streamline  continuity  is 
basically  a  function  of  the  streamline  positions  assumed.  In  integral 
form,  (25)  and  (26)  are  of  the  forms 

inY*=  TldX  +inCz  (29) 

where  XQ  is  arbitrary,  and  ln(C2)  is  the  constant  of  integration.  With 
the  boundary  condition,  Y  ■  1.0  at  X  »  1.0; 

o=l  ux  +AC4 

Then: 

>fnC2  =-('ldX 
and: 

Ma=  fidx-f'rdx 

Eq.  (29)  then  becomes: 

M2=jfldX  <3°) 

In  solving  (25)  and  (26),  the  sign  convention  for  the  flow  angles 
must  be  carefully  observed  (i.e.,  angles  are  positive  in  the  direction 
of  rotation).  In  addition,  when  integrating  graphically  or  numerically 
according  to  (30),  the  direction  of  integration  must  be  kept  in  mind. 
Integrating  toward  X  (1.0,  the  summation  must  be  multiplied  by  -1.0.  To 
facilitate  solution,  (30)  can  be  expanded  in  an  infinite  series: 

Y=e^'ldX  =  |  +  n^+-gi  +  ^  +  iS5r+  -  i  H-J,"IdX  (31) 

Method  of  Solution 
Loss  coefficients  and  gas  flow  angles  were  determined  for  both  the 
stator  and  rotor  using  the  method  of  Ainley.  [4,5]   In  exception  to  that 
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method,  however,  a  blockage  factor  due  to  trailing  edge  thickness  was 
included  in  angle  determinations.  Also,  leaving  losses  were  computed 
using  the  procedure  given  by  Vavra.   Computation  of  flow  angles  and  loss 
coefficients  is  given  in  Appendix  II.  Loss  coefficients  were  assumed  to 
be  independent  of  Mach  Number  and  independent  of  radius,  except  for  the 
influence  of  radial  variation  of  incidence  angle.  Angles  were  determined 
for  Mach  Numbers  of  0.5  and  1.0.  Since  only  the  mean  streamline  angles 
could  be  determined  for  M  =  1,  constant  differences  between  values  of 
flow  angles  for  the  mean  streamline  and  those  for  the  hub  and  tip  stream- 
lines were  assumed  for  all  Mach  Numbers.  Such  differences  were  taken  as 
those  for  M  ■  0.5.  Smooth  curves  with  inflection  points  at  M  =  0.75  were 
drawn  between  the  values  of  angles  at  M  =  0.5  and  M  ■  1,  following  Ainley. 
[4,5]   Only  the  converging  guide  vane,  full  rotor  configuration  has  been 
analyzed. 

Since  the  actual  curvature  of  streamlines  is  quite  uncertain  at 
best,  the  curvature  was  assumed  to  be  zero.  The  first  terms  of  (25)  and 
(26)  then  reduce  tos  cos2(  )  ds*/dX.  When  this  term  is  combined  with 
the  sine  squared  part  of  the  final  terms,  (25)  and  (26)  become: 

+  l'~  tfC  J  dX  (32) 


V? 


LLU*  cosfr    ,    C  cos1/?,,    djhb  j,  (1  -  CHtCo&l  \    ds? 
OC"  t%l  dX,l        £<,      j    dX, 

These  two  equations,  along  with  the  continuity  equation  in  the  forms  of 
5Ibid.,  Eq.  5(63),  p.  89. 
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(33) 


(16)  and  (20),  are  the  major  equations  to  be  solved  in  the  analysis. 

Conditions  ahead  of  the  stator  were  assumed  ass  constant  total 
enthalpy;  constant  entropy 5  and  uniform  velocity  with  neither  a  tangential 
nor  radial  component.  It  follows  from  the  adiabatic  flow  assumption  that 
total  enthalpy  is  constant  at  the  stator  exit  also,  and  dH/dX  at  station 
one  is  then  zero.  Similarly,  due  to  the  assumption  of  constant  entropy 
at  the  stator  inlet,  the  only  gradient  of  entropy  existing  after  the 
stator  is  that  due  to  friction  through  the  blading,  as  given  by  (28).  As 
a  result  of  the  assumed  axial-component-only  velocity  ahead  of  the  stator, 
the  loss  coefficient  for  the  stator  becomes  constant  at  all  radii  and 
equal  to  the  value  for  zero  incidence  angle. 

The  analysis  was  programmed  for  the  CDC  1604  computer,  using  FORTRAN 
60.  The  program  is  included  as  Appendix  III.  The  procedural  steps  accom- 
plished by  the  program,  PERFORM,  are  set  forth  in  succeeding  paragraphs 
of  this  section. 

Five  streamline  radial  positions  were  used,  including  the  hub  and 
tip  positions.  Inputs  for  the  program  are  listed  below.  Except  as  indicated, 
similar  data  is  read  into  the  computer  for  both  the  stator  and  rotor. 

1.  Radii  of  the  five  streamlines,  in  inches. 

2.  Values  of  gas  efflux  angles  in  radians  for  the  mean  streamline, 
and  for  Mach  Numbers  of:  0.5,  0.7,  0.75,  0.8  and  1.0. 

3.  Values  of  the  differences  in  flow  angles  from  the  mean  stream- 
line value  for  the  hub,  mean  (zero  difference)  and  tip  streamlines. 

4.  Rotor  loss  coefficients  and  rotor  profile  loss  coefficients  for 
five  values  of  the  ratio  of  incidence  angle  to  stalling  incidence  angle 
(i/is)s  6,  0.5,  1.0,  -0.8,  and  -1.5. 

5.  Values  of  the  throat  opening  dimension  for  the  hub,  mean  and 
tip  streamline  radii. 
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6.  The  number  of  sets  of  values  given  in  (7.)  below  for  which 
computations  are  to  be  accomplished. 

7.  Conditions  for  which  the  analysis  is  to  be  done:  Mach  Number  at 
the  mean  streamline  at  the  stator  exit  (M^  ) ;  stator  inlet  total  tempera- 
ture  (°R)  and  total  pressure  (psia.);  Mach  Number  assumed  for  the  mean 
streamline  at  the  rotor  exit;  and  turbine  RPM. 

Considerable  use  of  parabolic  equations  is  made  to  facilitate  compu- 
tation of  values  required  in  the  analysis.  Subroutine  PARAB  determines 
the  required  constants,  a,  b,  and  c,  for  the  general  equation: 

y  -  a   +  bx-tvcx2  (34) 

This  subroutine  solves  three  equations  in  three  unknowns  using  appropriate 
determinants.  Gas  angles  for  the  mean  streamline  at  the  stator  and  rotor 
exits  are  represented  by  (34)  for  two  ranges  of  Mach  Number:  0.5  to  0.75; 
and  0.75  to  1.0.  For  a  given  (input)  Mach  Number,  the  mean  streamline 
angle  can  be  determined.  Using  the  angle  differences  between  the  mean 
streamline  and  the  hub  and  tip  values,  the  flow  angles  for  the  hub  and 
the  tip  streamlines  are  computed.  Then,  the  parabola  for  flow  angle  as 
a  function  of  the  non-dimensional  radial  coordinate,  X,  is  found. 
Constants  of  Eq.  (34)  are  also  determined  for:  rotor  total  and  rotor 
profile  loss  coefficients  as  functions  of  i/is,  for  i/is  =  0  to  1  and 
i/is  s.  0  to  -1.5;  and  the  ratio  of  throat  dimension,  a,  to  the  mean 
streamline  value,  am,  as  a  function  of  X  for  both  the  stator  and  the 
rotor. 

Subroutine  STATOR  solves  (32),  without  the  dH/dX  term,  by  successive 
approximations.  Since  part  of  the  input  is  the  Mach  Number  of  the  flow 
for  the  mean  streamline  at  the  stator  exit,  and  since  the  flow  angle  for 
this  streamline  is  known,  the  axial  velocity  can  be  determined  from: 
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Tim"  \c 


^-1*1 


(35) 


and: 


L=WUi  (36) 


^i  =  V.„  "^  <*,.  (37) 


(It  should  be  noted  that  mass  flow  rate  is  yet  to  be  determined.)  In  the 
iterative  procedure,  the  first  approximation  is  obtained  with  ds*/dX  -  0. 
In  the  second  and  succeeding  approximations,  the  values  of  Y  determined 
in  the  preceding  iteration  are  used  to  find  ds*/dX.  In  determining  values 
of  Y,  using  (31),  the  series  through  the  third  power  is  used.  When  final 
values  of  Y  are  found,  conditions  at  all  streamlines  are  given  by: 

Va,  =  VaimY,  '  (38) 

v,  =  ~^r  <39> 

*  = T*  -  iAt  (40) 

T.u^-'fcf-  (41) 

T 

P.=  S.Hj?rr  (42) 

Vtt,=  l/.tanrf,  (43) 

Subroutine  FLOWREF  is  called  to  solve  (16)   for  the  stator.  This 
subroutine  checks  the  pressure  ratio,  p/pt  ,  to  determine  whether  or  not 
the  flow  is  choked  for  each  streamline.  For  choked  flow,  the  flow  func- 


tion, $  ,  is  set  equal  to  $*,  where: 
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rf=i 


$  =\rTTl       J  Tfl  (44) 

This  equation  is  obtained  from  the  expression  for   Q>   in  (10),  with: 

i  =  ten)*"  <«> 

Values  of  H    are  computed  using  the  expression  derived  in  App.  I: 


H 


Xe  ■  Xe  Xe       A         Xe 


+    .         .     4-     .        .      +■  ^^^r-r-  + 


f 


Xc-I  Jm-H  5m  fl  7m -H  9m  +  l    f    ||m-*-. 

~T  ~      i      T"xT  ;      xT~     xT~     xT 


+  -5-^ »-     -         ,      4-   .,    '    ,      + 


Xc-I     m  +-  I     3m  +  l     5m +1    7mfl    9rrv 
The  area  reduction  coefficient,  |  ,  is  then  found  from  (21).  Integration 
is  accomplished,  and  between-streamline  flow  fractions  are  determined. 
With  Wf  defined  as  the  fraction  of  the  total  flow  which  passes  between  the 
hub  and  any  other  given  streamline,  there  is: 

»  =  V*dX (47) 

*     fr*dx 

Here,  1^  is  the  value  of  the  integrand  of  (16).  The  actual  flow  rate  in 
lbm/sec  is  given  by: 

*=**^/F£t«*jx  (48) 

It  would  be  desirable  to  have  one-quarter  of  the  total  flow  pass 
between  each  pair  of  adjacent  streamlines.  However,  this  condition  is 
not  absolutely  necessary,  since  for  this  analysis  zero  curvature  was 
assumed  and  conditions  ahead  of  the  stator  are  completely  uniform.  As  a 
result,  the  conditions  at  the  stator  exit  can  be  used  to  establish  both 
overall  and  between-streamline  continuity  requirements.  Slight  differences 
from  multiples  of  one-quarter  for  the  flow  fractions  given  by  (47)  are 
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then  tolerable.  It  remains  mandatory,  though,  that  the  between-streamline 
flow  be  the  same  for  the  assumed  streamlines  at  the  rotor  exit  as  for  the 
corresponding  streamlines  at  the  stator  exit. 

Subroutine  R0T0R1  is  used  to  convert  the  absolute  flow  conditions 
ahead  of  the  rotor  to  relative  conditions.  For  angles  and  velocities: 

U|  =  (60)  (12) 5   f'  in  inches  (49) 

Wu,  =  \/«|-U,  (50) 

/*Mrfftj)  (51) 

Uz=U-^-  (53) 

In  (53),  r«  is  the  radius  of  the  given  streamline  at  the  rotor  exit. 
Other  necessary  values  are  obtained  from: 

H|= >T.  +  2$"  +  i{=^-  (54) 

\*  -£-  (") 

and: 


=  R(-^)T  (56) 


The  values  of  HE,  T^E,  and  Pt£  are  for  the  fluid  entering  the  rotor  along 
streamlines  chosen  for  the  stator  and  leaving  the  rotor  along  the  same 
streamlines  at  radii  assumed  for  the  rotor.  Similarly,  the  values  of 
entropy,  from  which  entropy  increases  through  the  rotor  are  determined 
for  the  several  streamlines,  are  those  for  the  matching  streamlines  at 
the  stator  exit.  These  entropy  values  are  the  increases  in  entropy 
through  the  stator  blading,  as  measured  from  the  constant  (base)  value  of 
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entropy  existent  at  the  stator  inlet.  Hence,  the  rotor  inlet  value  of 
entropy  can  be  considered  known  for  each  streamline  and  the  gradient  of 
entropy  with  respect  to  X2  at  the  rotor  inlet  can  be  found.  Also,  since 
He  is  constant  along  each  streamline  through  the  rotor,  the  gradient, 
dH]£/dX2,  is  known.  It  is  convenient  to  write: 

C*     _   £l~    So        _    S,  —  S,      ,     S.-So  r*       .    r* 


where  s^q  and  dsig/d^  are  known.  Then: 

o^i  _  olSn     _j_   dSio  (58) 

d  Xz       d  X2         d  Xji 

The  unknown  gradient  in  (58),  ds2]/dx2,  is  found  using  (28),  written  in 
terms  of  the  relative  flow  properties. 

The  rotor  loss  coefficient  required  for  the  solution  of  (28)  is 
determined  in  R0T0R1  by  using  the  parabolic  functions  determined  earlier 
in  the  program.  The  independent  variable,  i/is,  is  given  by: 

\=A-A  (59) 

Here,  /j'  is  the  blade  angle.  The  resulting  value  of  i  is  divided  by 
the  positive  stalling  incidence  angle:  (13°4  /57° 3) /radian.  Loss  coeffi- 
cients, both  total  and  profile,  for  the  rotor  are  taken  as  constant  and 
equal  to  the  values  for  i/ig  s  1  and  i/is  =  -1.5  should  the  value  of  i/is 
be  greater  than  1.0  or  less  than  -1.5  respectively. 

Subroutine  R0T0R2  solves  (33)  by  successive  approximations  using  the 
values  computed  in  R0T0R1  and  assumed  values  of  \Va ,  for  the  mean  stream- 
line  and  of  Y2  for  the  other  streamlines.  The  assumed  values  of  Y2  for 
the  first  iteration  are  estimates.  Then,  for  succeeding  iterations,  the 
preceding  iteration's  values  of  Y2  are  used.  For  convenience,  the  non- 
variant  part  of  dsjL/dX2  is  included  with  the  invariant  terms: 
-2tan/5.  d^-  |  sin  -.^  .  Since  the  gradient  of  Y2  at  the  rotor  exit  is  quite 
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severe,  terras  of  (31)  through  the  fifth  power  term  are  used. 

When  either  satisfactory  convergence  of  Y2  values  is  obtained,  or 
the  prescribed  number  of  iterations — 13 — have  been  completed,  thermodynamic 
conditions  at  the  rotor  exit  are  determined  using  (38)  -  (43)  written  for 
equivalent  properties.  The  absolute  tangential  velocity,  Vu_,  is  obtained 
from  (50),  written  for  the  values  at  station  two.  Satisfactory  convergence 
of  Y2  values  is  defined  ass    no  change  greater  than  0.005  in  successive 
iteration  values  of  Y2  for  any  streamline. 

Overall  continuity  is  next  checked  using  subroutine  FLOWREF  with  the 
rotor  exit  properties.  Eq.  (19)  as  used  in  the  computer  program  is  writ- 
ten to  include  the  flow  through  the  tip  clearance  area  as  part  of  the 
flow  between  streamlines  four  and  five.  Written  for  the  trapezoidal 
method  of  integration,  the  integral  of  (20)  iss 


(60) 


When  i  =  4,  the  first  term  in  the  brackets  becomes s 

Here,  k  is  the  tip  clearance  and,  therefore,  k/rm  is  the   AX2  of  the  tip 
clearance.  The  term  (X5-X4)  in  the  denominator  cancels  in  the  solution 
of  (60).  Eq.  (60)  must  yield  the  value  of  *rf/UmZrm)Rotor  giyen  b7  *rf 
determined  in  the  stator  solution.  If  the  rotor  flow  is  less  than 
required  (i.e.,  the  value  from  (60)  is  too  small),  the  assumed  value  of 
axial  velocity  for  the  mean  streamline  is  increased  and,  if  greater,  the 
mean  streamline  axial  velocity  is  reduced.  R0T0R2  and  FLOWREF  are  then 
called  successively  until  continuity  is  satisfied. 

When  overall  continuity  is  satisfied,  subroutine  SLINE  is  called  to 
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check  between-streamline  continuity.  This  subroutine  checks  the  fractions 
of  flow  between  streamlines  one  and  two,  one  and  three,  and  one  and  four 
at  the  rotor  exit  against  the  values  of  flow  fractions  between  the  same 
sets  of  streamlines  at  the  stator  exit.  If  the  rotor  and  stator  fractions 
for  all  three  sets  are  equal  within  tolerance,  the  solution  desired  is 
known.  If  inequalities  exist,  the  radial  locations  of  streamlines  two 
three  and  four  at  the  rotor  exit  are  adjusted  and  successive  approximations 
continued.  Also,  all  variables  required  in  R0T0R2  and  for  the  rotor 
equations  in  FLOWREF  are  adjusted,  as  required,  to  values  appropriate  for 
the  new  streamline  locations.  Then,  R0T0R2,  FLOWREF  and  SLINE  are  called 
as  before  until  the  equation  of  motion  and  both  continuity  requirements 
are  satisfied. 

It  is  appropriate  at  this  point  to  discuss  the  methods  used  to 
determine  required  gradients,  to  correct  radial  locations  of  streamlines, 
and  to  correct  gradients  for  new  values  of  X2.  For  variables  for  which 
parabolic  equations  were  determined  as  f(X)  initially,  the  gradients  are 
given  by? 

^=b+2c0(  (62) 

The  parabola  itself  and  (62)  need  be  corrected  only  for  changes  in  the 
radius  of  the  mean  streamline,  since  the  value  of  the  variable  at  a  given 
radius  does  not  change.  Then,  since: 

where  the  primes  indicate  the  new  values,  the  resulting  correction  factors 
required  to  permit  using  X1  in  the  equations  are: 

Fc,=  ^  (64) 

and: 
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F  -Fz 

fci~  rt"  (65) 

The  constants,  b  and  c,  of  the  parabolic  expressions  of  the  type  (34), 

and  occurring  in  (62),  are  then:  b«  ■  Fc;.bj  and  c*  =  FC2c  The  constant, 
a,  is  unchanged. 

New  values  of  streamline  values  of  X  =  X',  in  terras  of  the  old  value 
of  rm  are  obtained  using  flow  fractions  (see  (47))  and? 

X'-XiKta-^AcJ-^-  (66) 

Here,  the  positive  sign  corresponds  to  a  required  increase  in  X  due  to 
too  low  a  fraction  of  the  total  flow  passing  between  the  hub  streamline 
and  the  streamline  in  question;  the  minus  sign,  then,  corresponds  to  too 
great  a  flow  fraction  and  a  required  decrease  in  X.  The  situation  is 
illustrated  in  Fig.  16a.  Note  that  the  absolute  value  of  flow  fraction 
difference  is  used  in  (66)  and  that  the  gradient,  dX/dWf,  is  always  taken 
for  increasing  X.  Then,  in  the  nomenclature  of  Fig.  16a,  there  is: 

<Ji(   w,(T)-w,(I)  >       J>1  (67) 

The  value,  Xf,  obtained  yields  the  new  radius  to  be  used,  with; 

r  =  )(rr0  (68) 

Since  the  hub  streamline  flow  fraction  is  zero  and  the  tip  value  is  1.0, 
only  the  radii  for  streamlines  two,  three  and  four  need  be  corrected.  New 
values  of  X  are  given  by: 

New  values  of  U«  are  obtained  for  the  new  values  of  r2,  using: 

Ua=  U>  xJfsa-  (70) 
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The  curve  of  Fig.  16b  is  supposed  to  represent  some  arbitrarily- 
selected  f(X)  for  illustration  of  the  method  used  to  determine  gradients 
of  entropy  and  equivalent  enthalpy.  Gradients  at  the  hub  and  the  tip  are 
taken  as  the  straight-line  slopes  between  values  of  f(X)  for  streamlines 
one  and  two  and  streamlines  four  and  five  respectively.  Gradients  at  inter- 
mediate streamlines  (two,  three  and  four)  are  taken  as  the  average  of  the 
straight  line  slopes  between  values  for  these  streamlines  and  the  values 
for  streamlines  on  either  side.  All  slopes  are  computed  for  increasing  X. 

The  manner  in  which  values  of  s*,  HE  and  the  gradients  of  these  two 
properties  are  determined  for  corrected  streamlines  is  illustrated  in 
Figs.  17a  and  17b.  The  curves  of  Fig.  17a  show  the  possible  gradient 
variations  for  curves  which  might  pass  through  the  property  value  for  the 
streamline  of  interest.  Also  shown  are  the  straight  line  slopes  to  the 
property  values  for  the  streamlines  to  the  left  and  to  the  right.  Fig.  17b 
shows  the  increasing  X  portion  of  Fig.  17a  with  curve  segments,  straight 
line  slopes  and  actual  slopes  at  the  streamline  value  of  the  property. 
To  obtain  the  rate  of  change  of  the  gradient,  it  is  assumed  that  the 
straight  line  slope  is  the  slope  of  the  curve  at  the  value  of  X  one-half 
way  to  the  adjacent  streamline.  The  rate  of  change  of  slope  is  then  the 
straight  line  slope  minus  the  slope  of  the  curve  at  the  streamline  point 
on  the  curve,  divided  by  one-half  the  difference  in  X  values  between  the 
two  streamlines.  Note  that  the  positive  case  is  used  in  obtaining  the 
new  value  of  X,  whether  this  value  is  greater  or  smaller  than  the  old 
value.  If  X  must  be  decreased,  the  straight  line  slopes  to  the  left  of 
the  streamline  of  interest  are  considered  to  extend  to  the  right  of  the 
streamline  point.  Then,  once  the  rate  of  change  has  been  determined,  the 
gradient  for  the  new  value  of  X  for  the  streamline  is  the  old  gradient 
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FIGURE  17 


61 


plus  (  AX  =  Knew  "  Xold^  t:"nes  tne  rate  °*   change  of  slope  with  X.  It  is 
evident  that  the  correct  gradient  results,  since  the  sign  of  AX  deter- 
mines whether  the  plot  of  Fig.  17b  or  its  mirror  image,  with  respect  to 
the  old  streamline  value  of  X,is  used.  Similarly,  using  this  new  slope, 
the  value  of  the  variable  at  the  new  X  is  obtained  from  the  old  value 
plus  AX  times  the  new  value  of  the  gradient. 

Once  the  required  solution  of  the  analysis  is  known,  additional  pro- 
perties and  quantities  are  determined  for  each  streamline.  The  change  in 
total  enthalpy  is  given  by: 

Then,  sequentially,  the  following  are  determined: 

\-  \-  ^r  <72> 

%-.  *(¥)*"  (73) 

*.-*&■)*"  (74) 

The  temperatures,  T«.  ,  and  T2  ,  for  isen tropic  expansion  to  p?  from  pt 

and  PtF>  respectively,  are  given  by: 

r-i 

^"^(fc)  (75) 

The  overall  turbine  efficiency  (internal  efficiency)  is  then: 

Additionally,  the  stator  efficiency,  71   ,  and  the  rotor  efficiency,  77R, 
are  computed  using  (78)  and  (79)  below,  even  though  they  are  related  to 
the  input  loss  coefficients  by  (l-£  ).  The  values  obtained  in  the  output 
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provide  a  check  on  program  operation. 


Tt,  -  T, 
Tt0  ~  "His 


ns  =  >  "  ;'  (78) 


n  =  Tt*  ~  T*  (79) 

R   i;  -  t 

Stator  and  rotor  velocity  coefficients,  if     and  (|/ ,  are  obtained  as  the 
square  roots  of  (78)  and  (79),  respectively.  The  theoretical  degree  of 
reaction,  r  ,  and  the  so-called  head  coefficient,  ki<s,  are  given  by: 


■*_  T'is    ^Zis 


r*=^> — £i (80) 

'to  -  ]l;s 

it.    -  Ak  _  2  c,  (Tt .-T2is)  9T  (81) 

In  order  to  obtain  mass-flow-weighted  values  of  moment  and  horsepower, 

the  following  equations  are  used,  where  w^  is  the  flow  fraction  at  the 
rotor  exit,  and  the  subscript  w  indicates  the  mass-flow-weighted  value; 

fciwifoiW&Hjfal"*^  (82) 


2 


Then: 


ID  _  (A  U)*  J  W 

nl  —  — r^; 


550 


(83) 


Mr=  a  (yyo)  (84) 

Referred  values  are  obtained,  following  NASA  practice.  For  T=  1.4: 

(85) 
(86) 
(87) 


W  -  -TFT-  •  (88) 


^re(~ 

Mw  = 

S 

Nref  = 

N 

W      ,    — 

w/9" 

8 
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Here? 


e  = 

6  = 


T*o   _  Tu 


518.4 


-  _& 


r$TP 


14.7 


(89) 
(90) 


Results  and  Discussion 

The  performance  analysis  was  accomplished  for  assumed  stator  exit 

mean  streamline  Mach  Numbers,  M,  ,  between  0.9  and  1.4  in  steps  of  0.05. 

J-m 

The  general  range  of  RPM  used  was  11,000  to  20,000.  However,  at  assumed 
Mach  Numbers  of  1.1  and  above  at  low  RPM,  the  high  values  of  rotor  loss 
coefficients  caused  divergence  in  the  iterations  for  Y„  for  the  tip  stream- 
line. Also,  for  Mn  -  1.4  and  18,000  RPM,  and  for  H-,     =  1.35  and  18,500 
RPM,  the  pressure  ratio  across  the  rotor,  PtF/p2>  indicated  choked  flow 
in  the  rotor  and  overall  continuity  could  not  be  satisfied.  Since  no 
provisions  for  this  contingency  are  included  in  the  program,  17,500  and 
18,000  RPM  were  the  upper  limits  for  M^  =1.4  and  1.35,  respectively., 
Table  I  indicates  the  values  of  M^  and  RPM  for  which  the  analysis  was 
attempted  without  success. 

TABLE  I 

PERFORMANCE  ANALYSIS  INPUT  VALUES  FOR 
WHICH  RESULTS  ARE  NOT  AVAILABLE 


Mean  Streamline 
Mach  Number  - 
Stator  Exit  (M-^) 

RPM  (10~3) 

11 

12 

13 

14 

18 

19 

20 

1.1 

X 

1.15 

X 

X 









1 

1.2 

X 

1.25 

X 

X 

1.3 

X 

X 

X 

1.35 

X 

X 

X 

(18.5] 

X 

X 

1.4 

X 

X 

X 

X    X 

X 

X 
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Stator  inlet  conditions — total  temperature  and  total  pressure — were 

chosen  to  give  a  static  pressure  of  14.7psia   for  the  mean  streamline  at 

the  stator  exit,  except  that  limiting  values  of  660  R  and  43psia  were 

used.  At  the  higher  Mach  Numbers,  M-l  more  than  1.25,  stator  exit  plane 

pressure  was  a  partial  vacuum.  Fig.  18  shows  the  schedule  of  input 

variables  as  functions  of  Mi  . 

-••m 

Representative  output  from  the  program  is  included  in  Appendix  III, 
Since  most  of  the  printout  column  headings  are  self  explanatory,  no 
general  listing  is  given;  however,  App.  Ill  includes  such  description  of 
print-outs  as  is  required  for  clarity. 

Table  II  summarizes  the  results  of  the  performance  analysis  shown  in 
Figs.  19  through  30  for  ease  of  reference. 

The  analysis  showed  values  of  Y^  to  be  essentially  invariant  with 
M^  .  Consequently,  a  single  curve  of  Y^  versus  X^  is  representative  for 
all  values  of  M^  .  Moreover,  since  the  flow  discharge  angle,  #,'  ,  was 
assumed  to  be  constant  for  M±     =  1,  and  since  of-^  values  for  M^  =  0.9 
differ  from  the  M^  -  1  values  only  by  0.05  ,  a  single  curve  of  Vu-^  versus 
X^  is  well  representative  for  all  values  of  Mim. 

Rotor  exit  conditions  vary  with  each  of  the  rotor  inlet  condition 
variables.  Hence,  the  curves  of  Figs.  24  -  26  are  included  to  represent 
the  general  nature  of  rotor  exit  conditions,  as  given  by  the  analysis. 
However,  the  analysis  showed  that  rotor  exit  plane  pressure,  p2,  is 
essentially  constant  with  radius  for  each  set  of  rotor  inlet  conditions. 
The  maximum  difference  in  values  of  p2  for  the  several  streamlines  was 
approximately  0.5  psia  for  given  rotor  inlet  conditions.  It  is  evident 
from  Figs.  24  and  25  that  the  method  used  for  determining  gradients  of 
entropy  and  equivalent  total  enthalpy  is  satisfactory,  since  there  are  no 
marked  inflections  in  the  curves  of  s  and  HE.  It  is  interesting  to  note 
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TABLE  II 

SUMMARY  OF  PERFORMANCE  ANALYSIS 
CONTENTS  OF  FIGS.  19-30 


FIG. 

SHOWS 

VERSUS 

REMARKS 

19 

Referred  Flow  Rate,  w  «.- 

'  ref 

Pt0/Pl 

20 

Stator  Area  Restriction 
factor, ? 

Xl 

curves  for  pt  /p-^  - 
constant 

21 

vai/Ve 

Pt0/Pl 

22 

Vai/(vaira)  &  vUl/(Vai)m 

*L 

curves  essentially- 
independent  of  Mim 

23 

Pf/Pto 

*L 

curves  for  Mi  = 

J-m 
constant 

24 

s^,  s?*  &  their  sum 

X2 

Mx  =  1,  N  =  20,000 

25 

he 

X2 

same  as  above 

26 

Y2 

X2 

«lm  ■  i.o 

27 

(va2)n//e   .  Mlm>  &  Kis 

N/Ve 

28 

Internal  Efficiency,  1\  ^ 

kis 

for  hub,  mean  and 
tip  streamlines 

29 

Theoretical  Degree  of 
reaction,  r* 

kis 

mean  streamline 

30 

Referred  Rotor  Moment 
MR/5 

N//e 

mean  streamline 
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that  the  band  of  variation  of  Y2  with  RPM  for  constant  Mach  Number,  M^  , 
is  relatively  narrow.  This  observation,  although  valid  in  general,  is 
offered  with  qualification,  since,  as  stated  previously,  results  could 
not  be  obtained  for  certain  assumed  inlet  conditions  for  which  the 
iteration  for  Y„  became  divergent. 

Fig.  19  indicates  that  wref  decreases  with  increasing  pressure 
ratio,  Pt0/Pl>  whereas  for  choked  flow  it  could  be  anticipated  that  wref 
would  be  constant.  However,  it  is  evident  from  Fig.  20  that  the  reason 
for  the  decrease  is  that  values  of  the  area  restriction  factor  decrease 
as  Pt0/Pl  increases. 

Considering  that  normally  the  rotor  exit  plane  pressure  is  known, 
selected  performance  parameters  were  plotted  versus  Pt0/P2>  s0  that  by 
cross-plotting,  these  performance  parameters  could  be  shown  as  functions 
of  N/l/~9  and  k^s  for  different  values  of  Pt0/P2  =  constant.  The  resulting 
curves,  presented  in  Figs.  27  -  30,  are  designed  to  permit  performance 
estimations  for  selected  values  of  RPM,  inlet  total  temperature,  pressure 
ratio  Pt0/Po  anc*  either  p^  or  p2.  Figs.  27  -  30,  in  conjunction  with 
the  stator  and  rotor  gas  angle  information  included  in  App.  II,  and  Figs. 
19  and  21  -  23,  provide  sufficient  information  for  determination  of  mean 
streamline  performance.  Recommended  values  of  mean  streamline  radii  for 
use  in  such  a  performance  estimation  are  4.26  in.  for  the  stator  and  4.45 
in.  for  the  rotor.  Variation  of  conditions  with  radius  in  the  stator 
exit  plane  can  be  found  using  Figs.  20  -  23.  Moreover,  since  p2  is 
essentially  constant  with  radius ,  an  approximation  to  variation  of 
properties  in  the  rotor  exit  plane  can  be  obtained,  by  using  p2  =  constant 
and  estimated  values  of  Y2  taken  from  Fig.  26. 

It  might  appear,  on  the  basis  of  the  curves  of  Fig.  28,  that  the  mean 
streamline  values  of  internal  efficiency  are  not  adequately  representative 
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of  turbine  performance,  since  the  mean  streamline  curve  falls  appreciably 
closefr J  to  the  tip  curve  than  to  the  hub  curve.  However,  the  average  of 
the  values  of  t\^   for  all  streamlines  is  generally  only  0.001  to  0.002 
different  from  the  mean  streamline  value.  It  is  significant,  too,  that 
the  values  of  AH  for  all  streamlines  vary  little  from  the  value  for  the 
mean  streamline.  Consequently,  the  mean  streamline  values  of  internal 
efficiency,  as  well  as  of  the  other  performance  parameters,  can  be 
considered  excellent  indicators  of  turbine  performance  as  predicted  by 
the  performance  analysis. 

Unfortunately,  time  available  did  not  permit  thorough  investigation 
of  the  circumstances  leading  to  the  indications  of  rotor  choking.  Outputs 
from  the  computer  for  indicated  choked  flow  in  the  rotor  show  that  the 
manner  in  which  the  area  restriction  factor  and  the  flow  function  vary 
with  pressure  ratio  might  cause  too-early  indications  of  rotor  choking. 
As  the  pressure  ratio  approaches  the  critical  value,  the  rate  of  increase 
for  <£  decreases  while  the  rate  of  decrease  of  f  is  almost  constant,  so 
that  thein  product  .dncreases  only  slightly.  It  appears  that  the  rate  of 
increase  of  this  product  becomes  zero  at  the  critical  pressure  ratio  and 
then  decreases  with  increased  pressure  ratio,  Pt£/P2  ^   this  instance. 
It  will  be  recognized  that  this  decrease  of  the  product,  (j)  f  ,  was  shown 
to  account  for  the  non-constant  referred  flow  rate.  It  is  conceivable 
that,  should  rotor  flow  calculations  show  deficient  flow  rate  at  pressure 
ratios  reasonably  below  the  critical  value,  any  increase  in  axial  velocity 
yields  increased  pressure  ratio  with  no  increase  in  flow  rate.  Hence, 
large  increases  in  axial  velocity  and  pressure  ratio  result  with  succeeding 
iterations.  The  propriety  of  the  effect  of  ^  in  the  foregoing  is  beyond 
question.  However,  whether  the  apparently  overriding  influence  of  the 
area  restriction  factor,  f  ,  is  correct  is  considered  problematical.  It 
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is  pertinent,  though,  that  for  assumed  inlet  conditions  of  Mi  =1.35  and 
18,500  RPM,  for  which  the  analysis  showed  the  rotor  to  be  choked,  if  f 
were  held  constant  after  the  first  determination  of  rotor  flow  rate,  the 
maximum  rotor  flow  rate  given  by  using  $  would  not  have  increased  the 
rotor  flow  rate  to  the  value  required.  In  this  instance,  then,  the  choking 
of  the  rotor  was  not  caused  by  the  decrease  in  values  of  f  .  The  same 
observation  holds  also  for  the  other  choked-rotor  result  of  the  analysis. 
The  divergence  in  Y2  values  which  prevented  solutions  at  low  RPM  was 
probably  caused  in  part  by  the  higher  loss  coefficients  for  the  rotor  due 
to  incidence  angles  approaching  positive  stall  incidence.  However,  the 
loss  coefficients  by  themselves  should  not  have  strong  enough  influence 
to  generate  the  divergence.  Rather,  it  is  more  probable  that  the  initially 
assumed  values  of  Y2  were  enough  in  error  to  start  a  divergence  pattern 
after  the  first  iteration,  since  either  Y2  -ia/^2         or  that  ratio  squared 
is  a  multiplier  of  the  several  terms  of  the  integrand  for  determining  Y2 
(see  (33),  p.  48).  Considering  the  effect  of  J  in  (28),  p.  46  ,  it  is 
apparent  that  the  higher  the  value  of  £  ,  the  greater  is  the  subtractive 
term  in  the  denominator.  Still,  once  the  value  of  £  has  been  introduced 
into  that  term,  later  fluctuations  or  changes  in  the  term  are  due  to  Y| 
(or  to  the  ratio  given  above  squared  as  a  multiplier),  for  axial  velocity 
considered  as  a  constant.  It  is  then  apparent  that  high  rotor  loss  coef- 
ficients can  be  said  to  establish  a  certain  level  of  bias  in  the  entropy 
terms  which,  when  aggravated  by  relatively  large  fluctuations  in  Y2  values, 
can  ultimately  drive  the  subtractive  term  of  (28)  to  values  greater  than 
one.  Of  course,  this  latter  event  is  the  final  manifestation  of  the 
divergence  of  Y2  values,  since  the  logarithm  of  a  negative  number  is 
undefined.  Obviously,  the  subtractive  term  of  the  numerator  of  (28)  is 
affected  by  Y2  value  fluctuations,  but  only  about  one-half  as  much.  In 
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the  instant  analysis,  the  large  gas  efflux  angles  aggravate  the  'bias* 
produced  by  the  loss  coefficients,  since  the  cos"^  in  the  denominator 
serves  to  increase  the  subtractive  term. 

The  method  used  for  this  performance  analysis  appears  sound  and  is 
deemed  relatively  easy  to  use.  Convergence  of  iterations  is  quite  rapid 
except  as  noted  above.  It  was  found  that  the  rate  of  convergence  was 
increased  if  initial  input  RPM,  for  otherwise  constant  input  parameters, 
was  the  middle  value  of  the  expected  RPM  range,  with  subsequent  RPM  inputs 
increased  or  decreased  from  the  initially  assumed  value.  Obviously,  this 
method,  like  all  methods,  can  be  considered  no  more  accurate  than  the 
assumptions  used  in  the  solution.  Undoubtedly,  the  most  critical  values 
assumed  are  those  for  the  loss  coefficients. 

The  program  PERFORM  used  in  this  analysis,  although  apparently  giving 
satisfactory  results,  can  surely  be  improved.  It  was  noted  that  values 
of  a/am  for  the  mean  radius  slipped  somewhat  from  the  desired  value  of 
1  (to  0.9943  in  one  instance).  Since  this  ratio  is  expressed  as  a  para- 
bolic function  of  X,  and  since  slippage  was  due  to  the  correction  factors 
for  the  constants  in  the  parabolic  expression,  the  method  for  correcting 
rotor  exit  streamline  radial  locations  should  be  reexamined  with  a  view 
to  decreasing  the  noted  discrepancy.  Also,  it  is  considered  that  a  more 
appropriate  tolerance  for  flow  fractions  would  be  a  percentage,  say  one 
percent,  of  the  required  flow  fraction.  This  change  would  yield  more 
uniform  requirements  for  between-streamline  flow  rates,  rather  than  the 
presently-used  two,  one,  and  three-fourths  percent  tolerances,  used  in 
order,  for  the  one-fourth,  one-half  and  three-fourths  total  flow  points. 
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4,  Exhauster  Performance  Analysis . 

General 
Analysis  of  the  turbine  during  exhauster  operation  was  undertaken  to 
obtain  preliminary  estimates  of  total-to-total  pressure  ratios  achievable 
across  the  turbine  and  to  predict  shock  conditions  which  may  occur  in  the 
exhauster  piping.  The  analysis  conducted  was  for  one-dimensional  flow, 
based  on  the  momentum,  continuity  and  energy  equations.  Extensive  use 
was  made  of  polytropic  process  relations  and  the  non-dimensional  flow 
function,  $  ,  which  is  obtained  by  combining  the  energy  and  continuity 
equations.  The  relations  for  expansion  and  compression  processes  are 
given  below  without  proof.   For  an  expansion  process: 

=  '$  = 


a*  h 


2Y 


r-i 


(91) 


where : 


n  = 


(92) 


I-Hp(r-I) 

7jp  -  polytropic  efficiency 
r  -  specific  heat  ratio 
w  -  flow  rate  (lbm/sec) 

Tto-  stagnation  temperature  at  the  inlet  (°R) 
A  -  flow  area  at  the  station  of  interest  (in2) 
pto-  stagnation  pressure  at  the  inlet  (psia) 
p  -  static  pressure  at  the  station  of  interest  (psia) 
R  -  gas  constant  for  air  (ft-lb/lbm-°R) 
g  -  universal  gravitational  constant  (lbm- ft/lb-sec2) 
Similarly,  for  a  compression  process,  there  is  J 

^Vavra,  M.  H.,  Problems  of  Fluid  Mechanics  in  Radial  Turbomachincs 
(Rhode-Saint-Genese,  Belgium,  Von  Karman  Institute  for  Fluid  Dynamics,  1965) 
VKI  Course  Note  55a,  pp.  22-35. 
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w  Jii,  ;  r  = 


A%    19 


T-l 


T,  I  P, 


r-(i)n 


(93) 


Here,  the  subscripts  one  and  two  refer  to  inlet  and  station  of  interest 
respectiVffly.  The  polytropic  exponent  is  given  by: 

n  =  ,  _TT?Up)  (94) 

In  (93),  the  quantity  T  is  the  static  temperature  and  the  other  variables 
are  as  defined  and  have  dimensions  as  specified  for  the  expansion  process. 

Fig.  5,  p.  26  ,  shows  the  complete  exhauster  and  its  dimensions  and 
Fig.  13,  p.  37,  shows  the  exhauster  nozzle  and  the  flow  channel  from  the 
turbine  hood*  Fig.  5  also  shows  station  designations  and  subscripts  used 
in  the  analysis  and  the  general  locations  where  the  several  polytropic 
efficiencies  apply. 

Method  of  Solution 

Choked  flow  is  assumed  for  the  flow  through  the  stator  blade  row  and 
the  flow  through  the  exhauster  nozzle.  For  this  condition,  the  pressure 
ratio  in  (91)  is: 


=  (^TT-)n"  <95> 


Then,  (91)  becomes: 

For  choked  conditions,  flow  rates  must  be  determined  using  this  value  of 
^=  <£*•  however,  (91)  remains  valid  for  determining  pressure  ratios  lower 
than  the  critical  value. 

Fig.  31  shows  the  processes  involved  in  the  analysis  as  they  might 
appear  on  a  temperature-entropy  diagram.  Turbine  and  nozzle  flow  processes 
and  the  several  solutions  possible  downstream  of  the  nozzle  exit  plane  are 
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POSSIBLE  EXHAUSTER  FLOW  PROCESSES 
FIGURE  51 
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indicated.  Station  designations  used  in  Fig.  31  correspond  to  those  given 
in  Fig.  5.  At  station  two,  the  turbine  and  nozzle  flows  are  assumed  to 
be  mixed  and  the  combined  flow  is  assumed  to  have  uniform  properties 
perpendicular  to  the  flow  direction. 

As  previously  stated,  the  flow  at  station  one  consists  of  the  flow 
from  the  turbine  hood  plus  the  flow  through  the  exhauster  nozzle.  The 
former  flow  rate  is  the  sum  of  the  flows  through  the  turbine,  through  the 
plenum  labyrinth  seals,  and  through  the  shaft  labyrinth  seal.  Writing 
the  sum  of  these  flows,  w™,  in  the  order  named: 

wT  =  vvb  +  wt  +  Ws  (97) 


The  equations  for  w»   and  wg,  as  derived  in  Section  5,  are: 


£~9'    JX    J   R     >  ■/         4.0791 


vVs=$s4=-(?;        $=0,874 


Ta      v'F 


I-  CWPJZ 

l0+$ir.(P»/O 


(98) 


(99) 


The  flow  rate  through  the  turbine  is  obtained  from  (91),  where  (£  ■  $*  is 
given  by  (96) s 


wk=$*  Ai. 


JTTJ  R 


(100) 


For  the  expansion  process  from  the  hood  to  the  exit  plane  of  the  nozzle, 
it  is  convenient  to  writes 


*=M-ftr-fc) 


ry  +  l 


(101) 


Then,  putting  (98),  (99)  and  (100)  in  terms  of  $»s,  there  is,  with  (97): 


;J  T*.   R 


*t 


P.   Tt 


^-U+Sfc  +  Ut/r 


(102) 


Implicit  in  this  relation  is  the  assumption  that  the  total  conditions  at 
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the  inlet  to  the  plenum  labyrinth  seals  are  the  same  as  those  ahead  of 
the  stator  of  the  turbine.  Combining  the  left  side  of  (102)  with  (101): 

WrJTta 


Wr<l  If. 


nT-i 


f-W«-*A 


Pt« 


*d 


Mb) 


(103) 


or: 


I  Ptd  /  P.T  J 


(104) 


Here,  the  exponent,  i^,  is  the  value  corresponding  to  the  polytropic 
efficiency  assumed  across  the  turbine.  This  expression  can  be  written, 
with  (102): 


AS  [J 

fV 


nT+i 


flWO  "" 


Then,  introducing  the  right  side  of  (101): 


Ac 


(Pto/0 


nT+i  / 


tf-(tF=UHVUt^   (106) 


The  form  of  this  equation  used  in  the  analysis  is  obtained,  after  some 
rearrangement,  as: 


'  zr  /xy* 


2 

ftn-n) 

PlT\n- 

/P.T     |"J        I 

PtJ 

nT  +  l 


P.r 


UHMU^f 


(107) 


In  (107),  there  remain  two  unknowns:  p^_  and  P^,*  Hence,  it  is 
necessary  to  solve  the  equation  by  iteration,  using  desired  values  of  pt 
and  Tt  .  Although  two  solutions  of  (107)  exist,  the  supersonic  solution 
is  not  of  interest  in  the  analysis  and  only  the  subsonic  solution  is 
obtained.  With  an  assumed  value  of  P]_,  Pt,  is  determined  from  (107), 
the  condition  and  velocity  of  the  fluid  flowing  to  the  exhauster  from 
the  hood  are  established  from  the  relations: 
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nT-l 


Tt,  =  t  (^  ^  (108) 


VlT  =  /ZeJcP(TtJ-T)T)  (110) 

Eq.  (96)  with  (91)  gives  the  flow  rate  through  the  nozzle;  with  the 
asterisk  denoting  M  =  1.0  at  the  throat: 


Since  the  flow  rate  at  the  nozzle  exit  is  the  same  as  at  the  throat,  for 
the  exit,  station  1N,  there  iss 


*»1\  [X_    .  £  An 
K  Pt9  /  <? 


-  ^  Ail  -  /   ^     /JVP"  _  /R-f* 
"  ^  AIM  "/    T-l     [?J  [pj 


(112) 


Here,  the  exponent  nN  is  used  to  indicate  that  different  polytropic 
efficiencies  are  assumed  for  flow  to  the  throat  and  from  the  throat  to 
the  exit.  Eq.  (112)  can  be  solved  by  iteration  to  find  P1N.  Although 
both  a  supersonic  and  subsonic  solution  of  (112)  exist,  the  subsonic 
solution  was  not  considered,  since  the  back  pressure  required  for  a  shock 
at  the  exit  is  about  27  psia  for  an  isentropic  nozzle.  With  pn  deter- 
mined,  (109)  and  (110),  written  for  the  nozzle  variables  are  solved  to 
determine  the  temperature  and  velocity  at  station  1^, 

The  momentum  equation,  written  between  stations  one  and  two,  is: 

f^  14-  t K  -  iK- Mn-rt-5  (113) 

The  friction  force,  F^,  can  be  expressed  as  a  pressure  loss? 

F;=-f  Dz2Ap  =  A2AP  (114) 

This  pressure  drop,  in  turn,  can  be  written  ass 


84 


*P=f-5t4v  (115) 

Here,  L  is  the  duct  length  between  the  two  stations,  D2  is  the  duct 
diameter,  V  and  Jo  are  the  average  of  the  values  of  velocity  and  density 
at  stations  Ly  and  two,  and  f  is  the  friction  factor.  Conditions  for 
the  flow  through  the  outer  annulus  at  station  one  are  used,  because  the 
pressure  loss  is  considered  to  be  caused  primarily  by  the  friction  along 
the  duct  walls.  The  factor,  f,  depends  on  Reynolds  Number  and  pipe 
roughness.  The  approximate  range  of  Reynolds  Numbers  obtained  from 
preliminary  calculations  was  1.7  to  2.2  (10^);  and,  for  a  surface  rough- 
ness of  0.00015  ft.  for  commercial  steel  pipe,  a  surface  roughness  to 
diameter  ratio  of  0.00024  was  obtained.  With  these  values,  f-=  0.017  was 
determined  from  the  Moody  Diagram  given  by  Eschbach.  [l] 
Letting  pif  =  p^  -  Ap,  (115)  becomes: 

H*H  -(?k-PfK = wti-  ^  (U6) 

The  total  temperature  of  the  fluid  stream  after  mixing  is  given  by: 

Tt;,=   "'  V  ^  j         »=^»T  (117) 

Then,   from  continuity,   there  is: 
Solving  (118)  for  p2  and,  with  (117): 

Next,   combining  (119)  and  (116)  to  eliminate  p2: 

»»+»t  u  _  ktu  -JttU  =  a  p  _  J^flJ^kjLJyL %^\       (120) 
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Multiplying  (120)  by  V2,  expanding  the  last  term,  combining  like  terms 
and  transposing,  the  form  used  is  obtained: 


< 


wM-t-wT     w  R 


9 


-v, 


f%+f  ^t+A2pJ^RVO  (121) 


Lets 


B,=fO-,-y 


B,= 


tV..+  f.r+P^ 


and; 


B3=wR\ 


(122) 


(123) 


(124) 


The  solution  of  (121)  is  then: 


\l,= 


B, 


2B,  ± 


m2 


B, 


(125) 


This  equation  will  yield  both  a  subsonic  and  a  supersonic  solution.  With 
the  velocity,  V2,  known,  the  fluid  properties  at  station  two  are  obtained 
from  (119)  ands 


\  =  \  - 
K-  *  ft 


29rCf> 
r 

T-l 


(126) 
(127) 


For  the  subsonic  solution  of  (125),  the  fluid  is  compressed  in  the 
diverging  portion  of  the  exhauster  and  (93)  applies  ass 


=  JL  = 


2T 


T- 


m*  -ghT* 


(128) 


This  equation  is  solved  by  iteration  to  determine  p^«  Since  the  flow  must 
adjust  to  atmospheric  pressure  at  the  exit  plane  of  the  diverging  portion, 
if  the  value  of  p-j  obtained  does  not  equal  atmospheric  pressure,  the  value 
assumed  for  p-i  is  in  error.  For  P3<pa,  a  higher  value  of  p^  is  assumed 
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and  for  p-j}  pa,  a  lower  value  is  taken.  Iteration  is  then  continued 
until  the  value  of  piT  for  which  p3  ■  pa  is  found.  Values  of  fluid 
properties  and  the  velocity  at  station  three  are  then  determined  from: 


n,-l 


T3  -  T,  (f)  *  (129) 

l/3  =/ZqJc,(Tti-T,)  (130) 

The  supersonic  solution  of  (125)  gives  the  velocity  at  station  two. 
The  temperature  is  obtained  from  (126),  the  pressure  from  (119),  and  the 
total  pressure  from  (127).  Knowing  these  properties,  the  location  of  a 
normal  shock  can  be  estimated.  The  possibilities  which  must  be  considered 
in  the  one-dimensional  analysis  are: 

1.  a  normal  shock  at  the  inlet  to  the  diverging  portion,  followed 
by  subsonic  compression; 

2.  supersonic  expansion  from  the  inlet  to  the  diverging  portion  to 
some  axial  distance,  X,  downstream,  with  a  normal  shock  at  that  point  and 
subsequent  subsonic  compression;  and 

3.  supersonic  expansion  to  the  exit  plane  of  the  diverging  portion, 
followed  by  an  oblique  shock  outside  the  exhauster. 

Not  listed  in  the  above  are  supersonic  expansion  to  an  exit  pressure 
either  equal  to  or  less  than  the  full  expansion  pressure.  In  the  analysis, 
this  consideration  is  not  of  practical  interest;  however,  consideration 
of  these  possibilities  is  included  in  the  solution  procedure  for  complete- 
ness. 

Since  the  flow  must  adjust  to  ambient  pressure,  possibility  (1.)  is 
checked  first,  because  it  provides  maximum  pressure  recovery.  If  the 

resulting  value  of  P3  is  less  than  pa,  there  can  be  no  normal  shock  from 
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station  two  downstream.  Then,  any  shock  would  have  to  be  between  stations 
one  and  two.  Due  to  the  assumption  of  uniform  properties  at  station  two 
—  necessary  for  application  of  the  momentum  equation  between  stations 
one  and  two,  where  mixing  is  assumed  complete — no  definitive  information 
can  be  obtained  regarding  shock  location  upstream  of  station  two.  It  is 
worth  noting  that  should  mixing  be  completed  prior  to  the  assumed  station 
two,  the  pressure  recovery  would  be  less  than  that  for  the  situation 
considered  in  the  analysis.  For  the  purposes  of  this  analysis,  then,  the 
supersonic  solution  was  considered  non-existent  if  possibility  (1.)  did 
not  provide  sufficient  pressure  recovery. 

The  well-known  normal  shock  relations,  written  for  a  shock  at  station 
two  ares 


JL  -  ztM2z  -(t-i) 

Pi  T+l 


(132) 


T,   "~  /T  +  l  r+i  h.z\  1-L^.J 


,2.  _   M*  +  y^T 


M*  =  2T  J'  -  (134) 

pt2  -  p^_M^+i^i    r+i      ;  st*** 


Here,  the  subscripts,  *2*  and  fj2  ',  refer  to  upstream  and  downstream  of 
the  shock  respectively.  For  subsonic  compression  after  the  shock  at  two, 
there  iss 


w  (ft   pR        x    '  zr 


=  $, 


^j-Pj   /9     *3~ '  /  T"' 


Wtf-Qf* 


(137) 
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Again,  this  equation  must  be  solved  by  iteration  to  yield  P3  and  the  value 
obtained  is  compared  with  atmospheric  pressure.  As  noted  above,  the  solu- 
tion— no  shock  from  station  two  downstream — is  kaown  if  P3<Pa«  If  P3  =  Pa» 
the  solution  is  known  and  the  properties  of  the  fluid  at  station  three  are 
determined  from  (129)  -  (131),  with  appropriate  subscripts. 

For  P3>  p  obtained  from  the  shock-at-station-two-check,  the  next 
check  made  is  for  conditions  in  the  exit  plane  with  supersonic  expansion 
to  that  station.  The  expansion  equation  is  written  as? 


W    At    /_R_     ._    *  2T       (P3  \"e       /P3    ^    n 


--2-       ,  „     xHfiJlL 


(138) 


Here,  pt   is  a  fictitious  stagnation  pressure  determined  for  stagnation 
conditions  which  would  exist  for  a  converging-diverging  nozzle  with  a 
Mach  Number,  M2,  at  the  point  where  the  area  is  A2.  Then: 

With  P3  determined  from  iteration  of  (138),  T3  is  obtained  from  (109),  V3 
from  (110),  written  with  appropriate  subscripts,  and  stagnation  pressure 
from: 


\-*ftf 


r-i 

(140) 


Also,  there  is: 

M3=   /   ^t  (141) 

The  value  of  p,  obtained  corresponds  to  the  value  for  full  expansion.  If 
P3  .  pa,  the  solution  is  known  and  the  properties  and  velocity  have  been 
obtained.  If  P3  >  pa,  the  flow  is  underexpanded  and  the  flow  adjustment 
required  takes  place  outside  the  exit  plane  in  the  form  of  a  Prandtl-Meyer 
expansion.  In  this  instance  also,  the  solution  is  known  and  the  properties 
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at  the  exit  are  taken  as  those  for  full  expansion,  with  comment  that 

the  noted  expansion  occurs.  For  P3<p_,  the  practical  consideration,  two 

possibilities  must  be  considered:  an  oblique  shock  at  the  exit;  and  a 

shock  in  the  diverging  portion.  The  highest  pressure  for  occurrence  of 

an  oblique  shock  at  the  exit,  p    is  obtained  from  (132),  with  the  up- 

3 

stream  Mach  Number  equal  to  that  for  full  expansion.  If  the  pressure 
downstream  of  the  shock  is  less  than  or  equal  to  pa,  the  solution  is  known 
as  a  shock  at  the  exit.  Exit  conditions  are  then  taken  as  those  for  full 
expansion,  but  with  comment  that  an  oblique  shock  occurs  outside  the  exit 
plane. 

If  Pn _ < p  ,  it  is  necessary  to  determine  the  location  of  the  normal 
shock  in  the  diverging  portion.  The  process  of  solution  is  one  of  itera- 
tion, using  supersonic  expansion  to  an  assumed  distance  X  from  station  two, 
followed  by  a  normal  shock  and  subsequent  compression.  All  equations 
required  have  been  given  previously,  except  the  one  used  to  determine  the 

area  at  'X'.  From  the  geometry  of  the  diverging  portion  there  iss 

2 

The  iteration,  using  (138),  (132)  -  (136)  and  (137)  in  that  order,  and 
with  appropriate  subscripts,  is  continued  until  the  distance  X  is  found 
such  that  P3  ■  pa.  Exit  plane  conditions  are  then  determined  using  (129) 
-  (131),  with  the  subscript  9 H  v  substituted  for  two. 

Computer  Program 
Program  EJECTOR,  included  in  Appendix  IV,  performs  the  analysis 
described  in  the  foregoing  paragraphs  in  the  manner  and  in  the  order 
given.  Inputs  to  the  program  are  the  values  of  polytropic  efficiencies 
assumed,  the  number  of  different  inlet  conditions  for  which  solution  is 
desired,  and  the  sets  of  values  of  inlet  total  temperature  (°R)  and  total 
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4U  =  vri  =   tt 


rz  +  -f(r3-r2) 


(142) 


pressure  (psia).  The  program  obtains  the  subsonic  and  supersonic  solutions 
for  both  stator  profiles  (i.e.,  the, converging  and  the  converging-diverging 
nozzle  types)  for  each  set  of  inlet  conditions.  The  several  subroutines 
are  listed  below  in  the  order  called,  and  the  function  of  each  is  given. 

1.  NOZZLE.  This  subroutine  determines  the  conditions  at  the  exit 
of  the  exhauster  nozzle  by  iteration,  using  (111),  (112),  (109)  and  (110). 

2.  TURBINE.  For  an  assumed  value  of  pressure  at  station  lj>   for  the 
flow  from  the  turbine  hood,  this  subroutine  solves  (107)  for  pt  by  itera- 
tion, using  also  (98),  (99)  and  (100)  to  determine  required  flow  rates  (or, 
flow  functions,  ({)  ).  Eq.  (97)  is  used  to  find  the  total  flow  rate  from 
the  turbine  hood  to  the  exhauster  piping. 

3.  MMENTUM.  Eq.  (125)  is  solved  in  this  subroutine  for  the  subsonic 
conditions  at  station  two.  In  order  to  get  values  of  velocity  and  density 
at  this  station,  for  use  in  determining  average  velocity  and  density  re- 
quired for  calculation  of  the  pressure  loss,  (125)  is  solved  for  zero 
friction.  Pressure  loss  is  then  determined  using  (115),  and  (125)  is  again 
solved  for  the  new  value  of  j^  .  Eqs.  (117),  (119),  (126)  and  (127)  then 
yield  the  additional  properties  at  station  two. 

4.  SUBC0MP.  This  subroutine  solves  (128)  to  determine  exit  pressure. 
The  main  program  contains  the  required  check  for  P3  ■  pa,  assumes  a  new 
value  of  plT  as  required,  and  again  calls  TURBINE,  MMENTUM  AND  SUBC0MP. 
This  sequence  of  events  continues  until  the  required  condition  is  met. 

The  subsonic  conditions  at  station  three  are  then  determined  from  (129)  - 

(131). 

5.  MMENTUM.  After  the  subsonic  solution  is  known,  this  subroutine 
is  again  called  and  the  supersonic  solution  at  station  two  is  determined. 

6.  SHOCK.  This  subroutine  solves  the  shock  equations  for  the  super- 
sonic solution  properties  at  station  two — (132)  -  (136). 
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7.  SUBCOMP.  This  subroutine  is  again  called  for  the  properties 
determined  by  SHOCK  and  the  pressure,  p-j,  is  determined.  If  P3<Pa>  the 
solution  is  known  as  no  shock  at  station  two  or  in  the  diverging  portion. 
If  p-j  ■  pa,  the  solution  is  known  also  and  the  conditions  at  station  three 
are  those  given  by  (129)  -  (131).  In  the  event  that  p-j>  pa,  the  program 
proceeds  through  the  subsequent  steps. 

8.  SUPEXP.  Using  the  MMENTUM  solution  conditions  at  station  two, 
this  subroutine  determines  conditions  at  station  three  for  full  expansion, 
using  (139),  (138),  (109),  (110),  (140)  and  (141).  The  main  program  then 
checks  for  p-^  ■  pa.  If  the  equality  is  satisfied,  or  if  p-A  pa,  the 
solution  is  known  as  full  expansion  or  under-expansion  and  the  exit  plane 
conditions  are  taken  as  those  for  full  expansion.  The  output  includes  a 
Hollerith  print  of  the  condition  found  to  identify  the  type  solution 
obtained,  i.e.,  either  full  or  under-expansion.  For  P3<pa,  the  subsequent 
steps  are  taken. 

9.  SHOCK.  To  determine  the  maximum  back  pressure  for  a  shock  at 
the  exit  plane,  this  subroutine  is  again  called  for  the  conditions  found 
for  full  expansion  to  station  three.  The  main  program  then  checks  post- 
shock  pressure  against  ambient  pressure.  The  solution  is  known  as  an 
oblique  shock  at  the  exit  if  this  post-shock  pressure  is  equal  to  or 
greater  than  ambient  pressure.  Such  solution  is  specified  by  a  Hollerith 
statement  and  the  exit  plane  conditions  are  taken  as  those  upstream  of 
the  shock.  If  the  post-shcck  pressure  is  less  than  pa,  the  program 
continues  with  the  following  calls. 

10.  SUPEXP.  This  subroutine  is  called  for  an  assumed  shock  location, 
X,  and  with  an  area  given  by  (142),  to  determine  conditions  for  supersonic 
expansion  to  that  location. 

11.  SHOCK.  Conditions  downstream,  for  a  shock  at  X,  are  determined. 
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12.  SUBCOMP.  This  subroutine  is  called  next  to  find  the  exit  plane 
pressure  for  the  conditions  downstream  of  the  shock  at  X.  Comparison  of 
this  pressure  with  ambient  pressure  is  accomplished  in  the  main  program 
as  before.  If  p^  =  pa,  the  solution  is  known,  the  shock  occurs  at  X,  and 
conditions  at  station  three  are  given  by  (129)  -  (131).  If  p3>  p  ,  the 
distance  X  is  increased  and  if  P3<pa,  the  distance  assumed  is  decreased. 
Iteration  varying  X  continues  until  the  required  equality  is  obtained. 

It  must  be  noted  that  subroutine  NOZZLE  is  called  only  once  in  the 
above  series  of  steps.  Hence,  it  is  tacitly  assumed  that  any  slight 
expansion  or  shock  at  the  exit  plane  of  the  nozzle  is  insignificant. 
Additionally,  the  pressure  at  station  one  is  assumed  constant  and  equal 
to  the  pressure  for  the  flow  from  the  turbine  hood,  since  this  pressure 
is  used  in  the  solution  of  the  momentum  equation. 

Results  and  Discussion 

Solutions  of  the  exhauster  analysis  were  obtained  for  inlet  total 
pressures  from  44  through  36  psia  in  two  psia  increments,  and  inlet  total 
temperature  was  varied  between  660  and  760° R  in  25°  steps  for  each  value 
of  inlet  total  pressure.  The  values  of  polytropic  efficiencies  assumed 
are:  for  the  turbine,  f\-  =   0.75;  from  the  turbine  hood  to  station  L_, 
fld  =  0.925;  for  flow  to  the  nozzle  throat, r?N  =  0.975;  from  the  nozzle 

throat  to  the  nozzle  exit,7lM  =  0.9;  for  supersonic  expansion  in  the 

•U - -  -S s  4-*  ^"^ 

diverging  portion,  7?  ■  0.8;  and  for  subsonic  compression  in  the  diverging 

e 

portion,  72c  ~   0*7. 

Results  of  the  analysis  showed  process  pressures  to  be  essentially 
independent  of  inlet  total  temperature.  Of  primary  interest  is  the  vari- 
ation of  total  to  total  pressure  ratio  across  the  turbine  with  inlet  total 
pressure.  Such  variation  for  both  the  converging  and  the  converging- 
diverging  nozzle  stator  profiles  is  shown  in  Fig..  32. 
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VARIATION  OF  TURBINE  TOTAL  TO  TOTAL  PRESSURE  RATIO 
WITH  INLET  TOTAL  PRESSURE-EXHAUSTER  ANALYSIS 
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Values  of  velocities  and  temperatures  given  by  the  analysis  varied 
with  both  inlet  total  pressure  and  total  temperature.  However,  values 
of  supersonic  Mach  Numbers  obtained  for  station  two,  the  inlet  to  the 
diverging  portion,  were  independent  of  inlet  total  temperature.  Pressure 
in  the  exit  plane  for  a  shock  at  station  two  (maximum  pressure  recovery) 
was  also  independent  of  inlet  total  temperature.  Maximum  values  of  p,  of 
13.88  and  13.83  psia  for  the  converging  and  converging-diverging  guide 
vanes,  respectively,  were  obtained  for  the  maximum  inlet  total  pressure 
of  44  psia    Hence,  the  analysis  performed  indicates  that  any  shock 
adjustment  will  occur  in  the  constant  area  duct  upstream  of  station  two. 
This  result  would  seem  to  substantiate  the  expectations  from  general  two- 
dimensional  considerations  that  a  series  of  oblique  shocks  followed  by  a 
normal  shock  will  occur  in  the  duct. 

Values  of  referred  flow  rates,  w/e"/8  ,  obtained  from  the  analysis 
were:  3.535  lbm/sec  through  the  nozzle 5  1.438  lbm/sec  through  the  turbine 
for  the  converging  guide  vanes;  and  1.022  lbm/sec  through  the  turbine  for 
the  converging-diverging  guide  vanes.  The  totals  for  the  two  types  of 
guide  vanes,  in  the  order  given  above,  are  4.973  and  4.557  lbm/sec. 

Although  at  present  an  upper  limit  of  660°R  exists  for  turbine  inlet 
total  temperature,  during  future  tests  with  the  exhauster  operating  it 
may  be  found  desirable  to  increase  that  limit  to  prevent  undue  condensation 
and/or  freezing  in  the  turbine  discharge  flow.  Consequently,  provided  in 
Fig.  33  are  curves  of  minimum  inlet  total  temperatures  required  to  main- 
tain turbine  hood  temperature  ten  degrees  above  freezing  and  to  maintain 
turbine  hood  flow  temperature  at  the  nozzle  exit  plane  at  492°R.  These 
values  are  plotted  versus  inlet  total  pressure.  The  curve  for  T^  =  32°F 
is  common  to  both  sets  of  inlet  guide  vanes. 
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5.  Determination  of  Flow  Rates. 

General 

In  Section  2,  the  several  required  flow  rates  were  enumerated.  In 
addition,  the  method  for  determining  the  flow  rate  through  the  exhauster 
nozzle  was  specified  (see  p.  21).  To  estimate  this  flew  rate,  (111)  given 
in  the  preceding  section  can  be  used.  In  this  section,  the  equations  for 
use  in  determining  the  flow  rates  through  the  flow  measurement  nozzle,  the 
plenum  labyrinth  seals  and  the  shaft  labyrinth  seals  are  developed.  Equa- 
tions developed  for  the  flow  nozzle  and  plenum  labyrinth  leak  rates  are 
based  on  experiments,  but  the  equation  for  shaft  leakage  provides  flow 
rate  estimation  only. 

Test  Procedures 

The  configuration  for  flow  nozzle  calibration  tests  is  indicated  in 
Fig.  1,  p.  21.  The  removable  standard  orifice  assembly  was  installed, 
the  outlet  from  tank  number  two  to  the  turbine  and  the  two-inch  pipe  were 
blocked,  and  the  valve  for  controlling  the  flow  to  the  exhauster  nozzle 
was  closed.  The  standard  square-edged  orifice  of  the  removable  assembly 
is  fitted  with  both  vena-contracta  and  flange  taps.  Upstream  pressures 
were  measured  against  atmospheric  reference  in  inches  of  mercury,  while 
differential  pressures  were  measured  in  inches  of  water  using  a  portable 
manometer  with  a  range  of  70  in.,  and  with  0.1  in.  graduations.  Flow 
nozzle  instrumentation  is  shown  in  Fig.  3,  p.  24.  The  accuracy  of 
pressures  and  differential  pressures  was  taken  as  ±  0.03  in.  and  the 
accuracy  of  temperature  measurements  was  assumed  to  be  ±  1.0  R. 

Flow  nozzle  calibration  tests  were  made  at  essentially  constant 
supply  pressures  of  either  30  or  45  psia.  Flow  rate  was  controlled  by 
using  the  valve  at  the  exit  of  the  removable  standard  orifice  piping  to 
set  the  differential  pressure  accross  the  vena-contracta  taps.  Supply 
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pressure  was  maintained  by  adjusting  the  main  discharge  valve  of  the 
compressor  installation.  Differential  pressureif.  was  increased  in 
increments  of  ten  inches  of  water  for  increasing  flow  rates  to  the 
maximum  differential  pressure  obtainable,  then  decreased  by  five  inches 
of  water,  followed  by  incremental  reductions  of  ten  inches  of  water  to 
minimum  flow.  Approximate  ranges  of  flow  rates  of  the  tests  were  1.44  to 
3.67  lbm/sec  at  30  psia,  and  1.8  to  4.6  lbm/sec  at  45  psia. 

Instrumentation  of  the  two-inch  pipe  orifice,  which  is  used  for  deter- 
mining plenum  labyrinth  seal  leak  rates,  is  basically  the  same  as  that  of 
the  removable  standard  orifice  assembly.  Plenum  labyrinth  leak  rate 
tests  were  made  by  varying  the  supply  pressure  (approximately)  from  24  to 
44  psia.  Pressures  were  set  by  observing  the  turbine  plenum  reading. 
Supply  pressure  was  varied  in  increments  of  about  five  inches  of  mercury 
for  increasing  pressure \   then,  after  an  initial  reduction  of  2.5  in. Hg 
from  the  maximum  supply  pressure,  incremental  decreases  of  five  inches  of 
mercury  were  used.  Tests  were  made  without  the  exhauster  and  associated 
piping  installed.  Hence,  the  maximum  ratio  of  seal  inlet  total  to  static 
discharge  pressure  was  limited  to  approximately  3. 

Formula  Development 

The  basic  flow  equation  for  a  square-edged  orifice  or  flow  nozzle, 

n 

in  lbm/hr,  iss 


W=  359.1  CfrfKY.J/o.k,  (143) 


Here: 


D«  ■  orifice  diameter  or  nozzle  throat  diameter  (in.). 

CX  =  coefficient  of  thermal  expansion  based  on  D^  (dimensionless). 

7Stearnes,  R.  F.,  et.  al.,  Flow  Measurement  with  Orifice  Meters 
(New  York?  Van  Nostraad  Co.,  Inc.,  1951)  p.  6. 
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K  =  discharge  coefficient  (dimensionless). 
Y^  =  expansion  coefficient  (dimensionless). 
sO,  ■  fluid  density  at  upstream  tap  (lbm/ft  ). 
hw  -  differential  pressure  across  taps  (in.  H«0). 
The  equation  of  state  for  an  ideal  gas  isi 

A  -  W-  (144) 

Here,  p,  is  in  psia,  R  -  53.35  ft-lb/lbm-°R,  and  T-,  is  in  °R.  For  p^  in 
inches  of  mercury,  (144)  becomes? 

/*  ~         (53.35)  T, "  'L328  T  (144a) 

Inserting  (144a)  in  (143)  and  converting  to  lbm/secs 


W  =  0.115  DjrfKY.jSS-  (145) 

The  coefficient,  C**,  is  determined  from  the  linear  coefficient  of  thermal 
expansion,  &L  .  For  D~  -  orifice  diameter  at  elevated  temperature: 

l^j   =  [n-of^At^zzH-Z^fAO+tAttV  (146) 

Neglecting  the  higher  order  final  term  and  letting  (dI)  =  (XD?  : 

O(0l=[\  -HZ^(At)]D2Z  (147) 

From  (147)  it  is  evident  thats 

0(=   |  +  Zofe(At)  (148) 

For  both  standard  orifices,  made  of  type  304  stainless  steel, 
CX-  (9.5)10~6  °F"i,  giving? 

0(o  =  \   +  O.OOI9(iT^')  (149) 

g 
(The  value,  0.00193  was  used  in  calculations.)   For  the  flow  nozzle,  made 

8Ibid.,  p.  257. 


of  2024-T4  aluminum,     CX  =  (12.6)  10~6     °F~1,   and: 

«n  =  1  +  0.0  02  S2(1f^i-)  .  (150) 

The  discharge  coefficient,  K,  is  dependent  on  the  ratio  of  orifice 
or  nozzle  diameter  to  upstream  pipe  diameter  (i.e.,  j3  -   D^/D-^),  Reynolds 
Number  based  on  D^,  and  the  type  measuring  taps  used.  To  simplify 
computations  for  square-edged  orifices,  the  following  relation  may  be 
used:9 

K  K* 

(I  +  A/Ra)   =  (I  +  A/Ri)  (151) 

where: 

A  =  D,(830  -5000/?  +  9000/?2-  420^+  B)  (152) 

with: 

B  =  7==  for  flange  taps  (153a) 

B  -  -4= 100  for  vena-cont.  taps  (153b) 

y   Di 

The  Reynolds  Number  expression  can  be  written: 

Re  =  ijilfi-  (154) 

where  z*  is  absolute  viscosity  in  centipoises  at  flow  temperature  and 
pressure.  For  z*  independent  of  pressure  and  assuming  linear  variation 
of  z'  in  the  temperature  range  from  50  to  300° F,  there  is: 

a  =100*'=  1.9  +  0.24  (-^q-  -l)  (155) 

Putting  (155)  in  (154),  with  W  -  3600w,  and  defining,  X  =  Re(10-6): 

V  _  2.275*  (156) 

9Ibid.t  p.  65.   10Ibid.,  p.  212.   nibid.,  p.  64.   12Ibid.t  p.  335. 
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Let  i  ■  1  +  A/Rp.  Thens 


r(>\ 


t  =  I  +  ^P  (157) 

Eq.  (151)  then  becomes? 

-^  =  ^r  (158) 

For  flange  taps,  values  of  IT  =  K^,  corresponding  to  Re  =  00,  are 
tabulated.13  Eq.  (158)  for  flange  taps  is  then: 

K  =  !X  (158a) 

For  vena-contracta  taps,  this  shortcut  method  is  not  authorized  and  values 
of  Kqq  are  not  tabulated.  With  K  tabulated  as  a  function  of  Reynolds 
Number,1  an  iterative  procedure  using  these  variables  is  normally 
required,  with  flow  rate,  w,  as  a  third  variable.  However,  by  judicious 
choice  of  Re  ,  and  hence  K  ,  (148)  becomes  for  vena-contracta  taps: 

K  =  Mr-  U58t>) 

For  the  removable  standard  orifice  assembly,  Re*  was  chosen  as  (10  ), 
giving  K*  -  0.698.  For  a  Reynolds  Number  of  2(105),  the  value  of  K  deter- 
mined from  (158b)  differs  by  0.114  percent  from  the  tabulated  value,  while 
the  range  of  Reynolds  Numbers  for  tests  was  approximately  3.9(10  )  to 
1.3(10^).  Similarly  for  the  two-inch  pipe  square-edged  orifice,  Re  was 
selectedA(105),  giving  K*  -  0.6103.  At  a  Reynolds  Number  of  5(104),  the 
computed  value  of  K  differs  from  the  tabulated  value  by  about  0.1  percent, 
and  at  Re  =  2(105)  the  difference  is  about  0.13  percent.  The  range  of 
Reynolds  Numbers  in  the  tests  conducted  to  date  was  approximately  5(10  ) 
to  1.2(10^).  Tests  yet  to  be  conducted  with  the  exhauster  installed  will 

13Ibid.,  Table  19.    14Ibid.,  Table  20. 
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be  at  somewhat  higher  Reynolds  Numbers  than  1.2(105).  However,  since  it 
is  anticipated  that  the  labyrinth  seal  flow  will  become  choked  at  pressure 
ratios  (plenum  total  to  seal  discharge  static)  only  slightly  higher  than 
3:1,  Re  =  2(10$)  is  considered  a  realistic  upper  limit. 

The  expansion  factor,  Y^>  compensates  for  compressibility  in  the 
flow  through  the  orifice  or  nozzle  and  is  defined  ass15 

Y,  =  1-0-  (0.41  +  0.35/0^  (159) 

Here,   Ap  and  p^  must  be  in  consistent  units,  and  f   is  the  specific  heat 
ratio.  For  T=  1.4,  Ap  =  h^  (in.  H20),  and  px  =  ?1   (in.  Hg),  this 
equation  becomes: 

y  _  I  Q   0.41  +  0.3S/f 4    hw  nfi  v 

Eq.  (145)  can  be  written  in  the  form: 

w=Co<§Y,/5j&  (161) 

Here,  the  constant,  C,  includes  D^  and  either  Koo  ,  for  flange  taps,  or 
&*/"§*>  for  vena  contracta  taps.  Values  of  C,  along  with  values  for  A  of 
(152),  and  the  several  constants  in  the  expressions  for  Y^  and  X,  were 
determined  by  using  the  CDC  1604  computer.  Table  III  lists  the  values  to 
be  used  for  each  standard  orifice.  In  solving  (161),  it  is  necessary  to 
assume  initially  that  isl,  With  the  resulting  value  of  w,  X  is  deter- 
mined from  (156),  and  the  required  value  of  £  is  obtained  from  (157). 
Eq.  (161)  is  then  solved  using  this  new  value  of  f . 

Since  the  object  of  the  flow  nozzle  calibration  was  the  determination 
of  the  discharge  coefficient,  Kn,  (145)  is  written: 

K"  =   o.iis  tf^Y.tf.h..  (162) 

15Ibid.,  p.  71. 
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TABLE  III 

PLOW  EQUATION  CONSTANTS  -  SQUARE-EDGED  ORIFICES 
TRANSONIC  TURBINE  TEST  RIO 


'«* 


Orifice 
Installa- 
tion 

^ff 

Flange  Taps 

Vena  Cont  Taps 

Ox 

Oy 

0 

A 

0 

A 

EQ(l6l) 

EQ(152) 

EQ(l6l) 

EQ(152) 

EQ(156) 

EQ(l60 

Removable 

0.6995 

1.4257 

0.00218 

1.4405 

0.001756 

0.5559 

0.026 

2"  pipe 

0.5991 

0.0476 

0.000505 

0.04761 

0.0002249 

2.7561 

0.022 

EQUATIONS 
(161) 

(152) 
(156) 

i 

t 

(160) 
(155)     , 

;  I 

(149) 
(157) 


*  =  o^sy,  ih  E    s 

V     *i!      ' 

A  r  D2(650  -  5000^+    9000/ 2  -  4200/ 5+  B) 
x  ..    (6.516)   (5600)  IP"4  w    c  Q    jL 


Do  Z 


— 

■» 

] 

Yj  s  1.0  - 

0.41  +O.55/04 

(1.4)  (15.59) 

Z  s  1.9  +  0.24 

100 

<x=  1.0  +  0.00195 

r-t  -68i 

1(X 

) 

feir    -  •  _    ~ ..      kw 


■  1.0  -  c. 


J  rl 


S«i.o  +    MiQ"6) 

X 
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Then,  with  the  flow  rate,  w,  determined  from  (161),  crtn  given  by  (150), 
and  Y1  by  (159a),  KQ  is  found  from  (162). 

The  basic  method  used  for  estimating  leak  rates  through  labyrinth 
seals  is  that  developed  by  Egli.  [7]   The  equation  of  continuity  combined 
with  the  energy  equation  for  an  ideal  labyrinth  seal  (i.e.,  all  kinetic 
energy  is  destroyed  after  each  throttling)  is: 

w  =  o<q  4^-/4-  (l63) 

where : 


CP  = 


&    ,  ,  -A 


i-  (p/ptor 


n+  4in(>t./p) 


l-r2 


n+Jr^C/r) 


r 


(164) 


0<  =  discharge  coefficient  for  a  single  throttling  (dimensionless). 

A  =  cross  sectional  area  of  seal  passage  (in^). 

pt  =  total  pressure  at  seal  entrance  (psia). 

p  =  static  pressure  at  seal  discharge  (psia). 

Tt  =  total  temperature  at  seal  inlet  (°R) ;  Tt  =  constant  across  the 
seal  for  assumed  adiabatic  flow. 

R  =  gas  constant  for  air  ( f t-lb/lbm-°R) 

n  =  number  of  throttlings. 

r  =  overall  seal  pressure  ratio  (p/pt0)» 
To  account  for  a  non-ideal  labyrinth  seal,  a  velocity  carryover  factor, 
T*j  is  introduced  to  adjust  for  carryover  of  kinetic  energy  from  one 
throttling  to  the  next.  Eq.  (163)  then  becomes  % 

/Tt0  V  K 
For   6=  tooth  clearance,  t  =  tooth  width,  and  s  -  chamber  width,  (see 
Fig.  4,  p.  25)? 

q£=   f(S/t,  n)  ~  constant 
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T  =  f (  6/s,  n)  =  constant 

The  dimensionless  seal  pressure  ratio  function,  cp  ,  was  developed 
by  series  expansion  with  a  lower  limit  of  pressure  ratio  across  any  one 
throttling  of  0.6.  Since  preliminary  analysis  of  exhauster  performance 
showed  an  overall  pressure  ratio  across  the  shaft  labyrinth  seal  of  about 
0.55,  this  assumption  was  considered  valid  for  the  ten  throttlings  of 
that  seal.  However,  since  plenum  labyrinth  seal  leak  rates  were  deter- 
mined experimentally,  use  of  (165)  was  limited  in  this  instance  to 
attempts  to  correlate  experimental  data  with  the  method  used  and  to  deter- 
mine whether  test  data  for  overall  seal  pressure  ratios  of  0.34  to  0.6 
could  be  extrapolated  to  lower  pressure  ratios  of  0.16  to  0.34.  For  the 
two  identical  labyrinths  on  either  side  of  the  plenum,  (165)  can  be 
written,  introducing  the  non-dimensional  flow  function,  ^  : 

\=\^[j  =*«.«  (166) 

If  (X   and  f     are  in  fact  constant,  K  =  constant,  and  experimental  data 
for  higher  values  of  r  can  be  extrapolated  to  lower  values  of  r.  Since 
it  is  anticipated  that  the  flow  through  the  seal  will  be  choked  at 

ty  max  =  0»2798  ('  =  0.225),  the  range  of  extrapolation  required  would  be 
relatively  small. 

The  curves  of  o<  versus  5/t  given  by  Egli  [7]  make  determination  of 

0(  for  the  geometry  of  the  shaft  seal  somewhat  uncertain.  Consequently, 
the  value  of  o(  =  0.76  was  taken  from  experimental  data  presented  by 
Jerie,  [8]  where  (X   is  plotted  as  a  function  of  S/t  and  S/b,  Here, 
b  =  chamber  heighth  (see  Fig.  4,  p.  25).  Plots  of  T*  versus  $/s   given 
by  Egli  are  based  on  experimental  data  for  8 > 0.01  in.  Egli»s  citation 
of  Friedrich*s  tests  of  straight-through  labyrinths  with  S  =  0.006  to 
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0.01  in.,  and  with  two  and  three  throttlings,  indicates  that  somewhat 
smaller  values  of  T  than  those  given  by  Egli  would  be  proper  for  the 
geometry  of  the  shaft  seal  of  the  TTTR.  Considering  also,  ten  rather 
than  two  or  three  throttlings,  and  data  for  multiple  throttlings  given 
by  Jerie,  Friedrich»s  value  of  1.15  was  selected  for  T*  in  preference  to 
Egli's  value  of  1.25.  With  these  values  of  <X  and  T*,  the  final  formula 
for  estimating  shaft  seal  leakage  is,  with  A  =  TT  (1.2) (0.005): 

^  =  <#28t^U?s  (167) 

v  't. 

Then,  introducing  the  non-dimensional  flow  function,  $  s 

Since  (p  =  f(r),  the  shaft  seal  leakage  can  be  represented  by  a  single 
plot  of  $  versus  overall  seal  pressure  ratio,  r.  This  plot  is  shown  in 
Fig.  34. 

Results 

Since  the  tolerance  for  vena-contracta  taps  is  less  than  that  for 
flange  taps,  plotted  and  tabulated  data  included  herein  are  vena-contracta 
values  unless  otherwise  indicated. 

Fig.  35  shows  nozzle  discharge  coefficient,  Kn,  versus  Reynolds 
Number  as  determined  from  tests  of  8  September  and  13  September  1965.  It 
is  evident  from  this  plot  that  Kjj  can  be  taken  as  constant  at  Reynolds 
Numbers  greater  than  about  6(105) — (w  ^  2  lbm/sec).  Values  of  Kft  = 
constant  were  determined  by  the  method  of  least  squares  for  the  following 
data  groups:  8  September  tests  at  30  psia;  13  September  tests  at  30  psia; 
these  two  tests  combinedj  13  September  tests  at  45  psia  5  and  all  listed 
tests  combined.  Resulting  minimum  least-square-error  values  of  K^  are 
tabulated  in  Table  IV.  In  addition,  values  of  Kn  =  constant  for  vena- 
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6     7     8     9     10     II 
Reynolds  Number,  Re(lO^) 
FLOW  NOZZLE  DISCHARGE  COEFFICIENT 

TRANSONIC  TURBINE  TEST  RIG 
FIGURE  35 

TABLE  IV 

Flow  Nozzle  Discharge  Coefficients 

Transonic  Turbine  Test  Rig 


Test 
Date 

Nominal 
Pressure 

Data 
Points 

Flange 

Taps 

Vena  Cont  Taps 

Error  # 

Kn 

Error  4 

Kn 

8  Sep 

30  psia 

9 

.00828 

.9994 

.00599 

1.0038 

13  Sep 

30  psia 

12 

.01055 

1.0026 

.00684 

1.0012* 
1.0024 

8/13  Sep 

30  psia 
30  psia 

21 
21 

.01531 

1.0012 

.01071 

8/13  Sep 

45  psia 

12 

.01614 

1.0011 

.01633 

1.0012 

13  Sep 

45  psia 

12 

.00497 

1.0008 

.00864 

.9992* 

it     Error  figures  show  least  squares  average  difference  of  data  points 

from  Kn  for  minimum  error 
*  Recommended  values 
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contracta  data  for  the  individual  tests  are  indicated  in  Fig.  35.  Recom- 
mended values  of  K^  for  30  psia  and  45  psia  are  indicated  by  an  asterisk 
in  Table  IV.  Some  difficulty  was  encountered  with  water  collecting  in 
orifice  pressure  tap  lines  during  the  8  September  tests  (30  psia).  Although 
these  lines  were  continually  drained  during  such  tests,  the  data  are 
considered  inferior  to  those  obtained  at  30  psia  on  13  September,  when 
the  water  collection  problem  did  not  exist. 

Experimentally  determined  values  of  fy   and  $u  for  the  plenum 
labyrinth  seal  (see  (166))  are  plotted  versus  overall  seal  pressure  ratio 
in  Figs.  36  and  37  respectively.  It  is  apparent  from  Fig.  36  that  extra- 
polation of  Kj£   to  lower  pressure  ratios  would  be  of  questionable  accuracy. 
Hence,  subsequent  tests  with  the  exhauster  operating  will  be  required. 
To  provide  a  mathematical  expression  for  fo  ■  f  (r),  the  data  plotted  in 
Fig.  37  were  fit  with  a  parabola  by  the  method  of  least  squares.  The 
general  equation  used  in  the  curve  fitting  was: 

j  =  a  +  b$t  +  c$  (169) 

In  addition  to  experimental  data,  the  known  condition  that  at  r  -  1.0, 

$.  ss  0,  was  introduced.  Since  all  values  of  &  were  of  like  sign  (  +  ), 
symmetry  about  the  r  axis  was  assumed  by  setting  b  =  0.  The  tests  of 
20  October  1965  showed  considerable  variation  between  values  of  K«   and  h 
computed  from  data  for  increasing  supply  pressure  and  decreasing  supply 
pressure.  Consequently,  especial  care  was  taken  to  assure  flow  stabili- 
zation prior  to  recording  data  during  the  26  October  tests.  Although 
the  noted  discrepancy  was  not  eliminated,  it  was  reduced  materially. 
Hence,  the  20  October,  decreasing  supply  pressure  values  of  Qp  (corres- 
ponding to  Kjg   points  to  the  upper  right  in  Fig.  36)  were  neglected  in 
the  curve  fitting  process.  The  resulting  least-squares  second  order 
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polynomial  obtained,  for  which  the  least  squares  error  is  0.02079,  is: 

0  =  1.0-4  079,1;  (m) 

Solved  for  the  normally  dependent  variable,  &  *   (170)  becomes s 


^e  "  /  4.07  91  <170a> 

This  equation  is  plotted  in  Fig.  37  along  with  the  experimental  data. 

Also  shown  in  Fig.  37  are  segments  of  three  curves,  $»s  2KIP,  for  values 

of  Kjg   appropriate  to  three  different  ranges  of  seal  pressure  ratio. 

Data  reduction  required  and  the  solution  of  (164)  were  accomplished 

by  using  the  CDC  1604  computer.  Listed  below  with  their  purposes  described 

are  the  various  programs  used.  The  programs  themselves  are  included  in 

Appendix  V,  along  with  the  description  of  inputs  and  outputs  for  each 

program.  Since  these  programs  use  the  equations  of  this  section  in  a 

straight-forward  manner,  no  detailed  explanation  is  given. 

1.  ORIFICE.  Computes  the  constant  C  of  (161),  the  value  of  A  of 
(152)  and  the  constants  in  the  expressions  for  Y-^  and  X. 

2.  FLOCAL.  Determines  the  values  of  Kn  of  (162),  based  on  experi- 
mental data. 

3.  FL0CAL2.  Finds  the  least  squares  value  of  K  =  constant. 

4.  LABRINT.  Solves  (164)  for  pressure  ratios  from  0.025  to  0.975. 

5.  LABLEAK.  Determines  flow  rates  and  values  of  fy  and  fy  of 
(166),  using  experimental  data. 

6.  MINSQ2.  Determines  the  least  squares  parabola  for  plenum 
labyrinth  leakage  (169). 

Discussion  and  Recommendations 
Based  on  the  results  of  nozzle  flow  rate  experiments  conducted,  it 
is  recommended  that  a  value  of  0.9992  be  used  for  the  nozzle  discharge 
coefficient  in  (145).  However,  for  supply  pressures  of  30  psia  or  lower, 
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a  value  of  1.0012  would  be  more  appropriate. 

Eq.  (170a)  is  recommended  for  determination  of  plenum  labyrinth 
leak  rates  in  the  pressure  ratio  ranges  r  -  0.03  to  1.0.  Since  plenum 
labyrinth  leak  rates  (0.04  to  0.1  lbm/sec)  are  very  much  smaller  than 
turbine  flow  rates  (2.0  to  3.9  lbm/sec),  the  slight  errors  in  leakage 
flow  rates  introduced  by  the  use  of  this  equation  should  hare  negligible 
effect  on  the  accuracy  of  turbine  flow  rates.  Furthermore,  this  equation 
should  give  satisfactory  flow  rates  at  pressure  ratios,  r,  higher  than 
those  used  in  its  generation. 

For  exhauster  operation  with  a  supply  pressure  of  45  psia  to  the 
exhauster  nozzle  and  the  turbine,  the  estimated  shaft  seal  leak  rate 
given  by  (167)  is  about  0.003  lbm/sec,  while  flow  rates  through  the  - 
exhauster  nozzle  and  the  turbine  are  estimated  as  3.38  and  3.87  lbm/sec, 
respectively.  Using  Egli*s  [7]  values  of  discharge  coeffioient  and 
carryover  factor  (  ck  »  0.9  and  T  *  1«25),  shaft  seal  leak  rate  estimates 
would  be  about  50  percent  greater  than  estimates  given  by  (167).  However, 

ML 

plenum  labyrinth  leak  rates  predicted  using  values  of  ex  and  nf  deter- 
mined from  Egli  are  approximately  50  to  70  percent  greater  than  the  leak 
rates  determine  experimentally.  Hence,  (167)  is  recommended  for  predicting 
shaft  seal  leakage  flow  rates. 

The  recommended  flow  rate  equations  are  given  in  Table  V  for  ready 
reference. 

Tolerance  in  Flow  Measurement 

The  overall  error  tolerance  in  percent  in  flow  measurements  using 
square-edged  orifices  is  given  by  the  square  root  of  the  sum  of  the  squares 
of  "exponents  times  tolerances,  in  percent"  ascribable  to  the  individual 
items  in  the  flow  equation.  Hence s1" 

l6Ibid.t  p.  24. 
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TABLE  V 

SUMMARY  OP  FORMULAS  FOR  DETERMINING  FLOW  RATES 
TRANSONIC  TURBINE  TEST  RIG 


Desired 
Flow  Rate 


Nozzle 
(lbm/sec) 


Basic  Formula 


w  ■2.07?g(BTlj/Pl  h* 


Auxiliary  Formulas 


c<nr  l.o  +O.OO252 


tx-6^ 


100 


II  as  1,0  «  0.025 


Plenum  Labyrinth 
Seal  Leak  Rste 
(lbm/sec) 


w  s  0.116 


I2S 
7°^ 


*  -  /~iTo  - 1 
?  \/  4.0791 


«;■  p 


H, 


Shaft  Labyrinth 
Seal  Leak  Rate 
(lbm/seo ) 


w  a  0.01466  JkJip. 


/ 


Tt, 


$s  0.874  cp 


■s 


1.0  -  r* 
T 


16  +  ^  In  (1/r) 


« 


*t, 


;    r= 


.  °p 


I    '    r. 
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AM=  {[ZAfDz^+CAf*)]1*  [A(K)]l+  L  A  (  Y, )  }2 

4iA(hJ]2  +  [iA(Fl)]2  +  [iA(T()]2j!/2  (171) 

Tolerances  of  the  individual  items  for  the  Transonic  Turbine  Test  Rig 
installation  are  given  in  Table  VI. 

The  error  in  values  of  discharge  coefficients  of  the  flow  nozzle  and 
the  square-edged  orifice  must  be  approximately  the  same.  It  follows  that 
the  tolerances  in  flow  nozzle  flow  rate  measurements  must  be  the  same  as 
those  for  the  square-edged  orifice  used  in  calibration.  Hence,  vena- 
contracta-based  values  are  preferable,  due  to  a  lower  error  tolerance 
than  that  for  flange  taps. 

Tolerances  in  measured  labyrinth  leakage  rates  also  should  be  de- 
pendent on  the  tolerance  of  the  calibration  orifice.  However,  it  is 
apparent  from  Fig.  35  that  there  is  considerable  scatter  in  the  data. 
Hence  a  tolerance  arbitrarily  set  at  ±  1056  would  appear  reasonable. 

Error  tolerances  in  estimated  shaft  seal  leak  rates  are  difficult 
to  specify  exactly.  Here,  a  reasonable  tolerance  could  be  set  at  ±  20/6, 
based  on  the  qualifications  set  forth  in  the  literature  cited. 
6.  Calibration  of  the  Force  and  Moment  Measurement  System 

The  axial  force  and  moment  about  the  turbine  axis  of  the  stator- 
plenum  assembly  are  measured  by  Wiancko,  two-arm,  variable  reluctance 
bridge  force  pickups  using  a  phase  shift  frequency  modulated  oscillator. 
The  linear  range  of  the  oscillator  is  from  ten  to  12kc.  The  permissible 
loadings  of  the  two  capsules  are  50  lb.  and  600  in-lb  for  the  force  and 
moment  pickups,  respectively. 

Tests  were  conducted  to  establish  the  linearity  conformance  of  the 
measurement  systems  and  to  determine  whether  or  not  cross-coupling  exists 
between  the  two  systems.  In  addition,  readings  of  force  and  moment  were 
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TABLE  VI 

PERCENT  TOLERANCES  OF  FLOW  EQUATION  VARIABLES 
AND  PARAMETERS  -  SQUARE-EDGED  ORIFICES 
TRANSONIC  TURBINE  TEST  RIG 


a.     Independent  of  the  Orific 

e«    . 

A(cx.)  s  ±  0.05 

\ 

A(YX)  =  0  to±0.05 

for  0.0  <--l  <   0.015 
*1 

s  ±0.5  to±1.0 

for  0.015  <  5*. 
*1 

A(hJ  *  ±0'°?    (100) 

^P  )  n   ±o.o?   (100) 

Pl  i    ; 

b.  Removable  Standard  Orifice  Assembly; 

A(D2)  =   ±°~f||-  (100)«  ±0.0118    > 

A(K)  s;  ±  1.0  for  flange  taps 

s  ±  0.4  to  -  0*6  for  vena  oontracta  taps 

o.  Two-inch  Pipe  Orifice  I 

d(D2)  «  ±  2i2S2L  (ioo)  =  ±  o.o6o<5 

0.825 
A(K)  ■  +  1.0  for  flange  taps 

s  ±  0.4  to  ±  0.5     for  vena  oontracta  taps 
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taken  during  plenum  labyrinth  leak  tests  to  check  for  net  force  and/or 
moment  due  to  the  flow  through  the  plenum  labyrinth  seals. 

Fig.  38  shows  schematically  the  configurations  for  the  force  and 
moment  calibration  tests.  Standard  weights  were  used  to  obtain  known 
forces  and  moments.  The  following  description  of  procedures  applies  to 
the  force  system;  however,  the  same  procedure  is  used  for  moment  calibra- 
tion, except  that  maximum  load  is  600  in-lb  and  increments  of  load  are 
100  in-lb  rather  than  the  corresponding  values  of  50  lb  and  12.5  lb  for 
force  capsule  calibration.  The  zero  setting  was  obtained  with  no  load 
as  ten  kc  and,  with  the  full  load  of  50  lb,  the  bandwidth  was  adjusted 
to  set  the  reading  at  12kc.  Linearity  was  checked  with  one-half  full 
load  applied,  with  the  desired  reading  as  like.  If  this  reading  was  not 
obtained,  the  linearity  was  adjusted  so  that  the  new  reading  was  differ- 
ent from  like  by  three  times  the  original  difference  from  like,  and  in 
the  same  direction  as  that  difference.  The  zero  point,  bandwidth  and 
linearity  were  then  adjusted  in  order  repeatedly  until  the  desired  read- 
ings were  obtained  without  further  adjustment  (i.e.,  ten,  12  and  like 
for  zero,  full  and  one-half  loads,  respectively).  Starting  with  zero 
load,  the  force  was  increased  in  increments  of  12  Jit  to  the  maximum  load 
and  then  decreased  by  the  same  increments  to  zero  load,  with  readings 
recorded  for  each  load.  The  differences  in  readings  for  loading  and 
unloading  at  each  given  force  indicate  the  hysteresis  (and  non-linearity) 
of  the  system.  The  maximum  departure  observed  for  the  force  capsule  was 
eight  cycles  per  second  (0.2  lb)  at  the  one-half  load  point j  maximum 
departure  for  the  moment  system  was  five  cycles  per  second  (1.5  in-lb). 

Cross-coupling  between  the  two  systems  was  checked  by  applying  a 
fixed  load  to  one  system  aaad  then  applying  various  loads  to  the  other, 
while  watching  for  changes  in  reading  of  the  former.  Although  it  was 
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concluded  that  no  significant  cross-coupling  exists,  there  remained  some 
doubt  as  to  whether  the  moment  affected  the  force  reading.  The  uncertain- 
ty arose  because  unpredictable  changes  in  force  reading  were  noted.  These 
changes  in  reading  were  consistent  neither  in  magnitude  nor  direction; 
also,  by  jiggling  the  moment  calibration  arm  the  reading  could  often  be 
returned  to  the  proper  value.  It  was  concluded  that  the  variation  of 
force  readings  was  caused  by  the  indeterminate  friction  force  developed 
between  the  shaft  to  which  the  moment  calibration  arm  is  attached  and  the 
bearings  through  which  the  shaft  passes.  It  is  significant  in  this  regard 
that  full-load  displacement  of  the  capsules  used  is  only  0.05  in. 
Pertinent  also  is  the  fact  that  no  changes  could  be  detected  in  the  force 
reading  while  applying  moments  with  a  wrench,  and  with  the  aforementioned 
shaft  disengaged. 

During  plenum  labyrinth  leak  test,  no  significant  changes  in  either 
force  or  moment  readings  were  observed.  It  was  concluded,  therefore, 
that  the  flow  through  these  seals  produces  no  net  force  or  moment. 
6.  Turbine  Tests. 

General 

Tests  of  the  turbine  were  conducted  both  with  and  without  the  rotor 
installed.  Only  the  converging  nozzle  profile  stator  has  been  tested  and 
tests  with  the  rotor  installed  were  done  with  this  same  stator  and  the 
full-size,  uncut  rotor.  Axial  clearance  for  tests  conducted  to  date  was 
0.243  in.  The  gas  path  for  these  tests  was  the  basic  path  shown  in  Fig. 
10,  p.  34.  Rotor  tip  clearance  was  0.025  in.  Instrumentation  for  the 
tests  was  as  shown  in  Fig.  6,  p.  28,  except  that  the  total  pressure  probe 
line  apparently  became  disengaged  inside  the  machine.  Hence,  subsequent 
to  the  third  rotational  test,  the  total  pressure  upstream  of  the  stator, 
pto,  was  not  measured.  Values  of  pt()  were  determined  by  applying  a 
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correction  based  on  previous  rims,  to  the  static  pressure  at  the  probe 
location,  p0.  It  bears  repeating  that  the  pressure  on  the  downstream 
side  of  the  cover  plate  inside  the  stator  was  measured  during  stator  tests 
but  not  during  tests  with  the  rotor  installed. 

Prior  to  each  test,  the  force  and  moment  measurement  systems  were 
calibrated  for  proper  readings  at  zero,  full  and  one-half  loads,  as 
described  in  the  preceding  section.  Also,  prior  to  each  test  with  the 
rotor  installed,  the  rotor  moment  measuring  system  was  calibrated  for: 
minimum  reading  of  zero  at  50  in-lb $  and  maximum  reading  100  at  700  in-lb. 
Then,  the  applied  moment  was  increased  from  50  to  700  and  then  decreased 
to  50  (in-lb).  Readings  obtained  were  plotted  versus  applied  moment  to 
yield  the  calibration  curve.  During  all  tests  it  was  necessary  to  main-^ 
tain  a  tare  loading  in  the  axial  direction  to  keep  force  readings  in  the 
linear  range. 

During  stator  tests,  inlet  total  pressure  was  varied  from  minimum 
values  to  the  maximum  achievable  ahead  of  the  stator,  about  42.5  psia. 
Data  were  recorded  at  increments  of  supply  pressure  of  about  five  inches 
of  mercury. 

Since  it  was  planned  to  present  turbine  performance  data  in  the  form 
of  curves  for  constant  pressure  ratio  (total  inlet  to  static  discharge), 
for  tests  with  the  rotor  installed  the  dynamometer  was  set  initially  for 
maximum  power  absorption.  Supply  pressure  was  increased  until  the  value 
for  the  desired  test  pressure  ratio  was  obtained.  With  desired  pressure 
set,  and  with  maximum  power  absorption,  minimum  RPM  for  the  test  pressure 
ratio  resulted.  Then,  RPM  was  increased  by  reducing  power  absorption  and 
data  was  obtained  at  desired  values  of  RPM,  while  turbine  inlet  total 
pressure  was  maintained  essentially  constant.  Inlet  total  temperatures 
were  kept  below  200°F,  while  the  minimum  inlet  total  temperature  used — 
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computed  in  advance— was  that  value  for  a  static  temperature  of  42°  F  at 
the  rotor  discharge. 

Experience  showed  that  supply  air  pressure  and  temperature  and  RPM 
each  affected  the  others  to  a  considerable  extent.  Hence,  every  effort 
was  made  to  record  RPM,  rotor  moment,  stator  force  and  moment,  and  the 
several  pressures  and  temperatures  as  simultaneously  as  possible. 

Data  Reduction 

The  turbine  flow  rate  is  determined  by  subtracting  the  plenum  laby- 
rinth leak  rate  from  the  flow  rate  through  the  flow  measurement  nozzle. 
The  former  is  obtained  from  (170a)  and  the  latter  from  (145),  with  L, 
the  nozzle  discharge  coefficient,  equal  to  0.9992. 

As  stated  in  Section  2,  the  instrumentation  of  the  sta tor-plenum 
assembly  was  designed  to  permit  use  of  the  momentum  and  the  moment  of 
momentum  equations  to  determine  stator  performance.  In  vector  form,  the 
general  momentum  equation,  for  the  assumed  steady  flow,  is: 

(cM  -f  dm&  =  Jkpjfc,  +  f-n,ptdse  +F  (172) 

Similarly,  the  moment  of  momentum  equation  can  be  written? 

j  r,x,V,  dm,  -  f   rex\7e  dme  =  f-f,  xn.p.ds,  -t-|-rexnepedse  +  M      (173) 

The  integral  term  for  shear  forces  between  fluid  particles  in  the  exit 
plane  of  the  stator  has  been  omitted  from  each  of  these  equations.  Such 
forces  act  perpendicular  to  the  turbine  axis  and  would  have  no  effect  in 
the  momentum  equation  for  the  axial  direction.  Also,  their  effects  in  the 
moment  of  momentum  equation  would  be  small,  since  the  forces  at  either 
side  of  the  between-blade  flow  channels  oppose  each  other  and  should  be 
nearly  equal.  In  (172)  and  (173),  the  terms  F  and  M  represent  the  force 
and  moment,  respectively,  exerted  on  the  flow  by  the  blades  and  walls,  and 
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are  due  to  pressure  and  shear  forces. 

The  control  volume  used  consists  of  the  inner  flow  passages  for  the 
fluid,  with  openings  at  the  inlet,  station  »e»,  and  at  the  stator  exit, 
station  one.  Fig.  39  shows  these  stations,  along  with  the  several  pres- 
sures, dimensions  and  forces  pertinent  to  the  reducing  (172)  and  (173) 
to  the  forms  used  in  data  reduction.  Fig.  6,  p.  28  ,  shows  the  stator- 
plenum  assembly  in  more  detail. 

It  is  evident  from  Fig.  39  that  integrals  for  the  entrance  station 
contribute  no  axial  force,  since  the  pressure  acts  perpendicular  to  the 
axial  direction,  and  since  the  velocity  of  the  fluid  entering  the  control 
volume  has  no  axial  component;  too,  provided  that  the  pressure  and  velocity 
at  the  inlet  opening  are  symmetrical  with  respect  to  the  turbine  center- 
line,  no  net  moment  about  the  centerline  acts  on  the  flow  at  the  inlet. 
Hence,  the  integrals  containing  pressure  and  velocity  at  the  inlet 
reduce  to  zero  in  the  equation  for  the  moment  of  interest  also.  The 
expression  for  the  axial  force  then  becomes s 

Fa=mVa(+P|As  (174) 

and  that  for  the  moment  about  the  centerline  reduces  tos 

Ma=mrmyUl  (175) 

Here,  V„  and  V„  are  the  values  at  the  mean  radius,  r-,  . 
*1      ul  m 

The  force  and  moment  pickups,  affixed  as  shown  in  Fig.  39,  measure 
the  reaction  force  and  moment  acting  on  the  stator-plenum  assembly. 
Hence,  these  pickups  measure  the  negative  of  F&  and  Ma.  Letting  R^  and 
M-  denote  the  reaction  force  and  moment,  (174)  and  (175)  becomes 

Ri  =  -*Va,-P,As  (17*) 
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Mi  -  -  m  fm,VU|  (177) 

Included  in  the  force  measured  by  the  axial  force  capsule  is  any 
differential  force  acting  on  the  exterior  of  the  assembly.  As  shown  in 
Fig.  39,  ambient  pressure  acts  to  the  right  over  the  entire  case,  while 
acting  in  the  negative  direction  are  the  pressure  on  the  cover  plate, 
P]t»  the  pressure  plT,  acting  on  the  area  (Ac-Ap);  and,  presumably, 
ambient  pressure  acting  on  the  annular  area  (Ac»  -Ac) ,     Let  R  be  the 
reaction  force  measured  by  the  force  pickup,  taken  as  positive  in  the 
positive  ia  direction.  For  R  positive,  the  force  reading  given  by  the 
force  measurement  system  will  be  positive  (i.e.,  greater  than  tenfkc). 
Then,  since  the  sum  of  all  forces  in  the  axial  direction  is  zeros 

R*mVa,-PiA5  +  paAc-RT(Ac-AT)-piiAH  =  0  (its) 

Solved  for  the  desired  quantity,  Va, ,  and  substituting  w/g  for  m,  there 


is: 


Va,  =  -HR+  paAc-  rA5-  Rt(Ac-At)-  P|iA„1  Uts) 

Considering  that  a  positive  moment  acting  on  the  stator-plenum 
assembly  results  in  a  negative  reading  from  the  moment  measurement  system 
(less  than  ten  kc),  for  Mg  defined  as  the  measured  moment  for  the  stator: 

Ms=mr„,Vu,  (180) 

Solved  for  Vu- ,  and  with  w/g  for  ms 

V„-3--£  .  USD 

With  an  imposed  axial  load  (tare),  R  in  (181)  must  be  the  force  given  by 

the  reading  minus  the  tare  force. 

With  Va,  and  V„,  known  from  the  solutions  of  (179)  and  (181),  the 
*x     ui 
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momentum  mean  flow  discharge  angle  at  station  one  is  obtained  from: 

Of,  =  tern1  ^  (182) 

Then,  there  ares 

V,  =  ^  (183) 

T.*\--&  (184) 

VUf-f^  (185) 

and,  stator  efficiency  is  given  by: 

7l,=5^P (186> 

't0    '"is 

For  choked  flow  conditions  through  the  stator,  an  area  restriction 
factor,   f  ,  can  be  obtained  using  the  non-dimensional  flow  function,  <J? . 
This^area  restriction  factor  indicates  the  fraction  of  the  throat  area 
which  would  be  required  to  pass  the  actual  flow  rate  under  isentropic 
conditions  and  is  analogous  to  the  factor,   f ,  used  in  the  turbine 
performance  analysis  (Section  3).  This  restriction  factor  is  defined  as: 

£  _  w_  =  A* X  /q  =  JL  (187) 


Here,  $  ,  is  the  choked-flow  isentropic  value  of  f>  for  the  throat  area, 
A  ,  and  w  is  the  isentropic  flow  rate  for  choked  flow.  Then: 


The  non-availability  of  the  throat  pressure  is  the  reason  why  the  re- 
striction factor  can  not  be  found  for  non-choked  flow,  since,  for  such 
condition,  $.  must  be  found  from: 
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(189) 


In  addition  to  the  foregoing,  an  effective  exit  plane  area  was  deter- 
mined using  the  axial  velocity  given  by  (179)  and  continuity  considerations, 
The  continuity  equation,  solved  for  this  area  iss 


/u 


(190) 


When  compared  with  the  geometric  area  of  the  exit  plane  annulus,  this  area 
gives  an  approximation  for  the  blockage  of  flow  area  due  to  both  boundary 
layer  effects  and  trailing  edge  thickness. 

To  be  calculated  for  turbine  tests  ares  the  momentum  mean  values  of 
absolute  and  relative  velocities,  angles  and  properties  after  the  stator 
and  rotor;  turbine  horsepowers  and  the  efficiencies  of  the  stator,  rotor 
and  the  complete  turbine. 

Since  the  pressure,  pi  ,  is  not  measured  during  tests  with  the  rotor 
installed,  use  of  the  momentum  equation  (179)  must  be  considered  approx- 
imate only.  For  comparison  purposes,  however,  p-i  is  estimated  so  that 
V   can  be  obtained  from  (179)  for  use  in  (182).  Due  to  viscosity,  the 
fluid  between  the  rotor  and  the  cover  plate  will  have  some  mean  value  of 
angular  velocity,  which  can  be  expressed  as  a  fraction,  f ,  of  the  rotor 
angular  velocity.  Any  value  of  f  assumed  is,  of  course,  quite  conjectural. 
However,  it  is  known  that  the  fluid  adjacent  to  the  rotor  has  an  angular 
velocity  equal  to  that  of  the  wheel,  while  the  fluid  adjacent  to  the  cover 
plate  is  at  rest.  Fig.  40  shows  a  fluid  particle  of  unit  depth  into  the 
paper.  The  forces  acting  on  the  particle  are  the  centrifugal  force  due 
to  the  particles  angular  velocity,  OJ£  =  fa/  ,  and  the  pressure  force. 
The  pressure  forces  a vt 
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Fp=(P"|£r^-df)de-(p+|£^)(r-,f)de  +  2p^dr 


(191) 


Since  the  centrifugal  force  per  unit  mass  is  ox~,   r,  the  centrifugal 
force  acting  on  the  particle  is: 

FC<F=  u$l  r^ordBdr) 


(192) 


For  radial  equilibrium  of  the  particle,  the  pressure  forces  must  equal  the 
centrifugal  force.  Expanding  (191),  equating  (191)  and  (192)  and  dividing 
by  rd9  dr? 


$-=/"*' 


(193) 


Integrating  from  r  =  r^  ,  where  p  «  p-i  ,  to  r,  where  the  pressure  is  p: 

H  H 

P-R  =^¥-(rlO  (194) 


ors 


p=Pi-^K-r>) 


(195) 
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The  differential  axial  force  acting  on  the  cover  plate  is  2 

dFa  =  pdA  =  [Fi^/^«-ri)]2i.rdr  (196) 

This  equation,  integrated  from  r  =  0  to  r  »  rlu>   an(*  divided  by  the  area 
of  the  cover  plate  yields  the  expression  for  determining  the  average 
effective  pressure  to  be  used  as  PiTs 

r      2-nj  rdr 

o 

Integrated,  (197)  iss 

K,r jri (198) 


or: 

%=^-^f-jk>  (199> 

For  p^  in  psia,  yO  in  lb-sec */tlfif   and  r-.   in  inches,  the  factor,  144 
in.4  converts  the  second  term  of  (199)  to  psia.  Expressing  ur~   as  a 
fraction  of  the  rotor  angular  velocity,  (199)  becomes: 

Using  assumed  values  of  f  =  0.5  and  yO  -   0.002378,  the  pressure 
given  by  (200)  is  used  in  (179)  to  obtain  an  approximate  value  of  the 
flow  angle  at  the  stator  exit  and  the  fluid  properties  at  station  one. 
However,  the  values  obtained  are  not  used  in  subsequent  calculations. 
Rather,  the  conditions  at  one  are  obtained  from  continuity  considerations 
and  the  value  of  VUl  given  by  (181).  The  isentropic  expansion  temperature 
at  station  one  is  determined  from  (186).  An  iteration  process  is  then 
used.  The  velocity  at  one  is  given  bys 


V,  -  J  2qJcp(Tto-T,)  (201) 
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Values  of  T^  are  assumed  by  adding  increments  to  T^.  .  The  velocity  is 
determined  from  (201),  and  the  flow  angle  at  one  is  given  by: 

*,=  sin~'^j  (202) 

Then,  a  value  of  T-,   is  computed  for  comparison  with  the  value  assumed: 

-JL 

/o,  R 


T.'  =  ^rb—  (204) 


When  a  match  for  T^  =  T-j  is  obtained,  the  velocity  and  temperature  are 
known,  and  Va,  is  found  from: 

Vq,  =  V,  COS  of,  (205) 

The  solution  for  the  relative  flow  velocities  and  angles  ahead  of 
the  rotor  is  obtained  using  the  velocity  triangle  relationships: 

WUl=Vu,-U,  ;   U,=  arr;m/iz  (206) 

A   =  rcm'^  (207) 

W'  =  To^-  <208) 

Equivalent  enthalpy  is  determined  from: 

He  =  cf T,  +  $j  +  4i#  •>      u^=  ■»*.'«  <209> 

With  p9  known  as  ambient  pressure,  conditions  after  the  rotor  are 

obtained  by  the  iterative  process7~based  on  continuity?-- used  for  the 

stator,  but  with  appropriate  relative  flow  variables  and  equivalent 

enthalpy.  The  required  relative  peripheral  velocity,  Wu  ,  is  determined 

2 

from  the  measured  rotor  moment,  M^.  The  change  in  total  enthalpy  across 
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the  rotor  is  given  by: 

AH  =  H0-H2=i^^^  '     (210) 

Here,  AH  is  in  ft2/sec2.  Then,  from  Euler's  turbine  equation,  written 
for  positive  AH: 

V-^V.,-^  (211) 

The  required  value  of  Wu  is  then  found  from: 

WU^VU,-UZ  (212) 

The  isentropic  temperature  used  as  a  base  for  assuming  values  of  T2  is 
that  value,  1~  t   for  isentropic  expansion  from  Tt  and  ptF  to  the  pressure 
P2.  There  is  2 


JkU 


'-I 


Pt«=n(-fi  (213> 

and: 

r-i 

T^=Tt>fi)T  (214) 

When  the  assumed  and  computed  values  of  T2  agree  within  tolerance, 
the  relative  velocity,  flow  discharge  angle,  and  temperature  at  station 
two  are  known.  The  following  equations  are  used  to  complete  the  solution: 

0(,  =  tan-,f-Ss»)  (215) 


Vz  -  fe  (216) 

r 


ir 


Referred  values  are  obtained,  as  is  done  in  the  performance  analysis, 
using: 
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6  =5184  (218) 

&  =  -pb-  (219) 

MRre+  =  -^  (220) 

Nre^J^  (221) 

^  =  ^g  (223) 

k|s  =  "Ufft  '  Ahts  in  ft/W                                                                 (224) 


and: 


-*_  iC-hz 


Ah;. 


(225) 


The  quantity,  k-  ,  is  defined  as  the  head  coefficient  and  is  related  to 
the  often-used  velocity  ratio,  U-,/C0,  where  C2/2  =  ^h^s  by: 

-u-  -  El  (226) 


Co     \l   kis 

The  theoretical  degree  of  reaction,  r*,  is  the  ratio  of  isentropic  enthalpy 

drop  from  p^  to  po  (determined  along  the  line  of  constant  entropy  given  by 

stator  inlet  conditions)  to  the  isentropic  enthalpy  drop  from  p^  to  po. 

*-o 

The  several  efficiencies  desired  from  the  rotor  solution  are  found 
from: 

tfR=  tE  "  y  (227) 

for  rotor  efficiency;  and: 


*i=ThT=  T  -T  (228) 

""*      't0     '2iS 

for  internal  efficiency.  The  velocity  coefficients  are: 
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V! 


fc--t.-/ii  (229) 


,vZth 

Computer  Program 

The  data  reduction  computer  program,  TTTRPGS,  included  in  Appendix 
VI,  proceeds  through  the  various  computations  required  in  the  solution 
of  the  foregoing  equations  by  using  subroutines  described  below.  Inputs 
to  the  program  are;  the  number  of  test  data  sets  to  be  reduced;  the  date 
of  test  (month  and  day);  ambient  pressure  and  temperature;  the  number 
of  data  points  for  each  data  set;  a  value  for  the  variable, TYPE,  which 
indicates  whether  the  data  is  for  stator-only  (TYPE  =  0)  or  full  turbine 
(TYPE  ^   0)  tests;  and  all  experimental  data.  All  experimental  data  inputs 
are  the  values  as  recorded  during  the  test,  except  that  temperatures  are 
in  °F  rather  than  in  millivolts,  and  rotor  moments  are  in  in-lb  rather 
than  in  percent  of  maximum  moment. 

Subroutines  are  called  in  the  following  orders 

1.  FLORATE.  Determines  the  flow  rate  measured  by  the  flow  nozzle 
(145)  and  the  plenum  labyrinth  leak  rate  (170a)  and  (166).  The  turbine 
flow  rate  is  then  found  as  the  former  flow  rate  minus  the  later. 

2.  STATOR.  Solves  the  stator  equations  and  determines  the  proper- 
ties at  the  stator  exit.  Values  based  on  the  solution  of  the  momentum 
equation  are  computed  and  printed,  but  values  used  in  subsequent  calcu- 
lations are  those  obtained  from  iteration  using  continuity  considerations. 
Also  determined  by  this  subroutine  are  the  values  for  the  area  restriction 
factor,  stator  efficiency  and  stator  velocity  coefficient ,  Mean  stream- 
line radius  ii   taken  as  4.26  in. 

3.  ROTOR.  This  subroutine  converts  absolute  flow  conditions  at 
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one  to  relative  flow  quantities,  determines  AH,  and  performs  the  iteration 
using  continuity  considerations  to  determine  relative  gas  efflux  angle 
and  velocity,  and  temperature  at  the  rotor  exit.  The  other  thermodynamic 
properties,  flow  angles,  rotor  efficiency  and  rotor  velocity  coefficient 
are  computed  also.  A  value  of  4.45  in.  is  used  for  the  mean  streamline 
radius . 

4.  REFERD.  This  subroutine  converts  data  to  their  referred  and/or 
non-dimensional  forms. 

Outputs  from  TTTRPGS  are  the  thermodynamic  properties,  velocities, 
flow  angles  for  both  the  absolute  and  relative  flows  as  appropriate, 
along  with  the  non-dimensional  values  (•<•_,  efficiencies,  r  ,  etc.),  and 
the  actual  and  referred  values  of  rotor  moment,  horsepower,  RPM  and  flow 
rates. 

Results 

Three  tests  without  the  rotor  installed  were  conducted  on  27  December 
1965  and  on  5  and  6  January  1966.  Large  fluctuations  in  axial  force 
readings  were  noted  during  the  first  test  and  the  data  obtained  were 
unusable.  The  turbine  hood  was  removed  to  eliminate  the  backflow  caused 
by  impingement  of  discharge  air  on  the  hood.  All  subsequent  testing, 
including  tests  with  the  rotor  installed,  was  done  with  the  hood  removed. 
Data  from  the  second  stator-only  test  were  not  considered  valid,  since 
an  open  turbine  plenum  drain  line  was  found  at  the  conclusion  of  the  test. 
Consequently,  only  the  data  from  the  third  test  form  the  basis  for  results 
presented  for  tests  without  the  rotor  installed. 

The  turbine  was  rotated  for  the  first  time  on  12  January  1966,  and 
five  additional  tests  were  made  between  that  date  and  16  February.  Rotor 
moments  measured  during  the  first  two  such  test  were  obviously  too  high, 
since  unbelievable  turbine  efficiencies  resulted.  Results  of  these  two 
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tests  were  considered  valid  only  for  stator  performance  considerations* 
Turbine  tests  were  discontinued  during  the  fifth  test  run  on  28  January, 
due  to  suspected  bearing  failure  in  the  air  dynamometer  unit.  Upon 
disassembly  of  that  unit,  the  bearings  were  found  to  be  satisfactory. 
However,  the  dynamometer  rotor  had  been  scraping  against  its  housing* 
Moreover,  the  spline  at  the  dynamometer  end  of  the  quill  shaft  connecting 
the  rotor  to  the  dynamometer  unit  was  badly  fretted,  indicating  possible 
misalignment*  After  reassembly  and  reinstallation  of  the  dyanmometer, 
the  turbine  was  realigned  and  testing  was  resumed*  It  was  necessary  to 
seat  the  dynamometer  spline  (male)  as  deeply  as  possible  in  the  quill 
shaft  to  maintain  adequate  undamaged  spline  bearing  and  shear  areas* 
Dynamometer  bearing  failure  was  again  suspected  on  18  February  during 
preparations  for  the  seventh  test.  At  this  time,  the  dynamometer  was 
returned  to  the  manufacturer  for  disassembly,  inspection  and  necessary 
repair*  Hence,  testing  was  again  interrupted  pending  return  of  the 
dynamometer  and  receipt  of  the  new  quill  shaft  which  had  been  ordered* 
The  turbine  became  operable  again  on  15  April*  However,  it  was  too  late 
for  further  turbine  performance  tests  for  inclusion  in  this  study. 
Consequently,  zero  speed  and  zero  torque  tests  were  made  as  the  final 
tests  for  consideration  herein. 

It  is  pertinent  to  consideration  of  the  test  results  that  the  total 
pressure  probe  pressure  line  parted  inside  the  stator-plenum  assembly  at 
the  conclusion  of  the  third  test.  Values  of  p*.  used  in  data  reduction 

''O 

for  subsequent  tests  were  determined  by  applying  a  correction,  deduced 
from  previous  tests,  to  the  static  pressure  readings,  p0,  obtained  at  the 
probe  location. 

Results  of  the  sta tor-only  and  the  turbine  tests  are  shown  in  Figs. 
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41  through  49.  For  ease  of  reference,  Table  VII  Indicates  the  contents 
of  these  figures.  Also  shown  in  Figs.  41  through  49  are  selected 
results  of  the  performance  analysis. 

Discussion  and  Conclusions 
As  is  evident  from  Fig.  41,  stator-only  tests  yielded  values  for 
stator  efficiency  of  0.9,  based  on  the  continuity  equation  iteration,  and 
0.93,  based  on  the  momentum  equation.  The  difference  between  values  of 
the  gas  efflux  angle  (  o(, )  given  by  the  two  methods  of  solution  (see 
Fig.  42)  reflects  the  difference  in  efficiency  values,  since  common 
values  of  tangential  velocity  were  used  in  computations  performed  by  the 
two  methods.  The  values  for  OC-.  used  in  the  performance  analysis  are 
seen  to  definitely  favor  the  continuity-based  values,  whereas  the  esti- 
mated stator  efficiency  used  for  the  performance  prediction  was  0.9095. 
Any  error  in  the  continuity-based  values  can  be  due  only  to  the  stator 
exit  area  used  in  computations.  This  area  was  the  stator  exit  annular 
area  minus  the  product  of  blade  heighth,  number  of  blades,  and  the 
dimension  for  the  projection  of  blade  trailing  edge  thickness  on  the 
stator  exit  plane.  Although  accuracy  of  this  exit  plane  area  is  dependent 
on  manufacturing  tolerances,  the  area  must  be  considered  quite  accurate. 
The  momentum-based  values  are  very  critical  with  respect  to  the  pressure 
terms  and  also  depend  on  the  axial  force.  For  a  typical  data  point 
(Pt  /Pi  "  2.61),  the  fluctuation  in  the  latter  quantity  over  a  period  of 
about  ten  seconds  was  ±  1.25  lbs  from  the  average  value  of  24.65  lbs. 
With  the  flow  rate  of  3.44  lbm/sec,  a  variation  of  ±9.4  ft/sec  results. 
Considering  that  the  values  of  axial  velocity  obtained  for  this  test 
point  were  404  ft/sec  from  the  momentum  equation  and  334  ft/sec  from 
continuity  considerations,  the  fluctuation  in  axial  force  measurements 
appears  to  be  relatively  unimportant.  It  must  be  considered  then  that 
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TABLE  VII 


SUMMARY  OF  THE  TURBINE  TEST  RESULTS 
PRESENTED  IN  FIGURES  41  THROUGH  49 
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either  one  or  more  of  the  pressures  used  in  the  momentum  equation 
computations  is  in  error  or  that  one  or  more  of  the  initial  assumptions 
concerning  the  areas  upon  which  the  several  pressures  act  is  not  valid. 

Subject  to  the  limitation  imposed  by  the  assumption  in  the  performance 
analysis  of  the  loss  coefficients  being  constant  with  radius,  it  can  be 
considered  that  the  pressure  distribution  from  hub  to  tip  given  by  the 
analysis  would  indicate  that  the  average  of  hub  and  tip  valuesof  pressure 
should  yield  a  representative  value  for  the  pressure  in  the  stator  exit 
plane,  p^.  In  addition,  the  pressure  pattern  given  by  the  analysis — low 
at  the  hub  and  high  at  the  tip—was  observed  during  the  tests.  For 
instance,  for  a  mean  streamline  stator  exit  Mach  number  of  one,  the 
performance  analysis  gave  values  (in  psia)  of  13.03  and  16.21  at  the  hub 
and  the  tip,  respectively,  while  matching  test  values  are  13.97  and  15.84. 
Considering,  however,  that  the  performance  analysis  average  pressure 
occurs  at  a  smaller  radius  than  that  of  the  mean  streamline,  an  average 
effective  pressure  obtained  by  integration  of  pressure  times  differential 
area  would  be  higher  than  the  average  of  the  hub  and  tip  pressure.  It 
should  be  noted  that  such  higher  value  would  reduce  the  value  of  axial 
velocity  obtained  from  the  momentum  equation,  which  is  the  effect  required 
for  sta tor-only  tests.  The  entire  stator  exit  plane  annular  area,  18.38 
in.  ,  is  acted  upon  by  the  pressure  p^.  For  a  change  of  70  ft/sec  in 
axial  velocity,  the  difference  between  values  given  previously,  a  value 
of  pi  only  about  0.42  psia  greater  than  the  average  of  hub  and  tip 
pressure  would  be  required.  Unfortunately,  a  pressure  survey  in  the 
stator  exit  plane  would  be  required  to  accurately  determine  the  radial 
pressure  gradient  which  would  be  needed  for  finding  the  force-average 
pressure.  The  performance  analysis  pressure  distribution  could  be  used 
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as  an  approximation,  however.  Further  consideration  of  the  pressure 
terms  which  occur  in  the  momentum  equation  indicates  that  it  would  be 
desirable  to  know  the  pressure  acting  on  the  flow  director  ring  down- 
stream of  the  stator  exit  plane,  on  both  the  annular  area  at  the  sharp 
break  at  the  stator  exit  plane  and  the  annular  area  at  the  downstream 
end  (see  Fig.  39,  p.  122). 

Any  comparison  of  stator  efficiency  and  gas  efflux  angle  values 
given  by  the  momentum  equation  and  from  continuity  considerations  for 
tests  with  the  rotor  installed  would  be  subject  to  the  considerations 
given  above  for  the  rotor.  However,  it  is  evident  from  Fig.  42  that  for 
the  rotor  tests  the  values  of  axial  velocity  obtained  from  the  momentum 
equation  are  too  low  (i.e.,  oC-.  too  large).  Of  some  interest  is  the 
fact  that  the  radial  gradient  of  pressure  obtained  with  the  rotor 
installed  was  somewhat  more  severe  than  that  for  stator-only  tests.  For 
the  conditions  considered  in  the  preceding  paragraph,  the  values  of  hub 
and  tip  pressures  measured  during  tests  with  the  rotor  were  13.5  and 
16.41  psia,  respectively.  These  considerations  serve  to  emphasize 
previously-made  comments  that  the  pressure  acting  on  the  cover  plate 
must  be  known. 

The  parabolic  nature  of  the  equation  for  the  radial  distribution  of 
pressure  acting  on  the  cover  plate  indicates  that  pressure  taps  should 
be  installed  at  the  center  of  the  cover  plate  and  at  a  radius  of  0.707 
times  the  cover  plate  outer  radius.  The  redesigned  thrust  and  torque 
flexure  will,  of  course,  yield  the  force  acting  on  the  cover  plate,  as 
previously  discussed.  Even  though  this  flexure  and  pressure  taps  in  the 
cover  plate  effectively  give  the  same  information,  installation  of  both 
is  considered  desirable. 

It  should  be  noted  that  the  spread  of  stator  efficiency  values  for 
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tests  with  the  rotor  installed,  Fig.  41,  is  not  reflected  in  values  of 
Of-,,  Fig.  42.  This  situation  is  considered  to  be  caused  by  the  difficulty 
noted  previously  of  obtaining  all  readings  rapidly  enough  to  preclude 
effects  of  variation  of  test  conditions.  Since  provisions  for  automatic 
data  logging  are  nearing  completion  at  the  time  of  this  writing,  this 
problem  should  be  alleviated  in  the  near  future. 

Of  considerable  interest  is  the  increased  stator  efficiency  with 
the  rotor  installed,  which  can  be  estimated  as  0.93  to  0.94, over  that 
without  the  rotor.  This  result  belies  the  assumption  that  conditions 
downstream  do  not  affect  upstream  conditions.  Since  it  seems  unlikely 
that  stator  blade  profile  losses  are  reduced  by  the  presence  of  the  rotor, 
it  is  conjectured  that  the  rotor  acts  to  reduce  the  effects  of  secondary 
flow  losses.  Considering  the  sharp  break  in  the  gas  path  in  the  stator 
exit  plane,  it  seems  probable  that  the  presence  of  the  rotor  may  provide 
sufficient  flow  guidance  in  that  area  to  reduce  losses  due  to  vortex 
formation  as  the  flow  spills  over  into  this  relatively  large  region. 
The  increased  severity  of  pressure  gradient  in  the  stator  exit  plane  for 
tests  with  the  rotor  installed  over  that  observed  for  stator-only  tests 
indicates  that  the  flow  pattern  in  the  area  is  altered  by  the  presence 
of  the  rotor.  Installation  of  pressure  taps  in  the  flow  director  ring 
at  this  break  could,  perhaps,  give  some  insight  into  the  nature  of  the 
flow  alteration,  as  well  as  being  desirable  for  application  of  the 
momentum  equation. 

The  curves  of  Figs.  43-46  illustrate  the  performance  of  the  turbine 
as  measured  in  tests  three  through  six  with  the  rotor  installed  and  with 
an  axial  clearance  of  0.243  in.  Some  evidence  of  the  effects  of  non- 
simultaneous  readings  exists  in  the  scatter  in  values  of  internal  effi- 
ciency, Fig.  43,  whereas  the  moment  curves  of  Fig.  46,  which  depend  only 
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on  flow  rate  and  measured  torque,  do  not  exhibit  undue  scatter.  Constancy 
of  referred  flow  rate  is  indicated  by  Fig.  44,  with  wref  =  1,375  lbm/sec* 
The  variation  of  theoretical  degree  of  reaction  indicated  in  Fig.  45  is 
considered  quite  reasonable  for  an  impulse-type  turbine.  Figs.  48  and 
49  are  included  here  for  illustrative  purposes  and  discussion  in  the 
next  section,  with  respect  to  comparison  of  performance  prediction  and 
test  results. 

Both  internal  efficiency,  Fig.  43,  and  rotor  efficiency,  Fig.  47, 
vary  with  incidence,  angle  in  a  manner  characteristic  of  an  impulse-type 
turbine,  i.e.,  with  a  narrow  range  of  incidence  angles  for  high  effi- 
ciencies. Except  for  the  differences  between  maximum  rotor  efficiency 
values  for  the  different  pressure  ratios,  the  curves  of  Fig.  47  could 
be  replaced  by  a  single  curve. 

The  values  of  rotor  gas  efflux  angle  are  not  depicted  graphically, 
since  this  angle  varied  over  a  small  ranges  -67.1°  to  -67.8°,  with  the 
higher  value  applicable  to  higher  pressure  ratios,  pt  /p2»  In  addition, 
the  stator  area  restriction  factor,  f  ,  was  essentially  constant,  vary- 
ing from  0.944  to  0.952,  and  with  the  variations  indicative  of  no  dependence 
on  stator  pressure  ratio. 
7.  Discussion  and  Recommendations. 

Comparison  of  experimental  and  performance  prediction  results 
indicates  that  the  most  serious  deficiency  in  the  performance  analysis 
is  in  the  variation  of  the  area  restriction  factor,  Fig.  20,  while 
experiment  yielded  an  essentially  constant  value.  The  effect  is  indicated 
in  the  reduced  referred  mass  flow  rates  after  the  flow  is  choked  given 
by  the  performance  analysis,  Fig.  19.  The  effect  is  evident  also  in  the 
agreement  between  the  predicted  curve  of  referred  axial  velocity  and 
experimental  values  at  lower  pressure  ratios,  pt  /pn,  with  subsequent 
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departure  at  higher  pressure  ratios,  Fig.  48.  Apparently  £  should  be 
taken  as  constant  after  the  flow  is  choked,  so  that  w  f  will  be  constant. 

Quite  obviously,  any  performance  prediction  is  no  more  accurate 
than  the  values  used  for  stator  and  rotor  efficiencies.  The  differences 
in  magnitude  between  predicted  and  observed  values  of  internal  efficiency, 
rotor  efficiency,  and  stator  efficiency  with  the  rotor  installed  serve  to 
emphasize  this  fact.  The  high  deflection  through  the  rotor,  129° ,  makes 
estimation  of  7Z„  quite  arbitrary,  since  very  little  is  known  about 
turbine  losses  for  rotors  with  deflections  of  greater  than  about  110- 
115°.  General  practice  in  gas  turbine  design  limits  the  turning  angle 
to  105°.  Predicting  the  positive  stalling  incidence  angle  is  subject 
to  the  same  difficulties  attending  the  choice  of  rotor  efficiency.  From 
Fig.  47,  the  stalling  incidence  angle  for  the  transonic  turbine  can  be 
estimated  as  about  five  or  six  degrees,  whereas  the  predicted  value 
used  was  13.4  .  The  effect  of  this  difference  is  evident  in  Figs.  43, 
47  and  45.  The  much  greater  low-loss  ranges  for  the  prediction  curves 
of  Figs.  43  and  47  are  clear  indications.  However,  the  manner  in  which 
the  theoretical  degree  of  reaction  is  influenced  (Fig.  45)  requires 
further  explanation.  One  must  first  imagine  that  the  predicted  curve 
of  Fig.  45  is  moved  upward  so  that  its  position  corresponds  to  the  higher 
efficiency  (over  predicted)  obtained  from  the  tests.  Then,  assuming 
that  the  loss  coefficient  remains  constant,  whea  positive  stalling 
incidence  is  reached,  the  smaller  the  value  of  this  incidence  angle, 
the  higher  will  be  the  RPM  at  which  the  loss  coefficient  becomes  constant. 
This  consideration  means  that  the  negative  reaction  portions  of  the 
prediction  curves  of  Figs.  29  and  47  would  flatten  out  (towards  zero 
reaction) . 

The  gas  efflux  angles  used  in  the  performance  analysis  proved  to 
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agree  very  well  with  those  determined  experimentally.  Fig.  42  indicates 
almost  exact  agreement  for  the  stator,  while  the  predicted  values  for 
the  rotor  angle  of  67.76°  to  69  agree  favorably  with  the  values  of  67. ]? 
to  67.8  observed  experimentally.  Further,  the  exit  Mach  Number  for 
which  the  value  of  67.8  was  obtained  was  about  0.6,  and  the  corresponding 
value  of  /& 2   used  in  the  performance  prediction  was  67.85  . 

The  discrepancy  between  predicted  and  measured  values  of  referred 
rotor  exit  axial  velocity,  Fig.  49,  can  be  attributed  largely  to  the 
differences  noted  previously  between  predicted  and  experimentally 
determined  values  of  the  stator  area  restriction  factor. 

It  is  considered  that  the  addition  of  instrumentation  recommended 
earlier  in  this  paper  is  necessary  prior  to  further  testing.  In 
addition,  it  is  recommended  that  attempts  to  correlate  the  performance 
prediction  program  as  given  herein  (or  a  similar  program)  be  continued, 
by  assuming  different  loss  coefficients  and  different  radial  variation 
of  such  loss  coefficients.  It  is  pertinent  in  this  regard  that  a  thesis 
reporting  on  cascade  tests  of  the  converging  guide  vanes  and  the  double 
circular  arc  rotor  blades  will  be  published  shortly. 
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APPENDIX  I 
FORMULA  DEVELOPMENT  FOR  TURBINE  PERFORMANCE  PREDICTION 
A.    Equation  of  Motion  for  Relative  Flows. 

The  equation  of  motion  for  relative  flow  in  vector  form,  (2), 
p.  39,  is: 


VHR  =  Wx(VxW  +  zTJj)  +  TVs 


(2) 


Expanding  the  several  terms  of  (2)  gives; 


VxW= 


ie  Ilk  4.  i    lik  +.  i  lik 


r   3^ 


9§ 


rW„ 


is 

r 


d 


Wx(\7xW)  = 


-     i        (Wr 

-   le  la* 


dr  ) 


7~vae         a?    / 


lr 
Wr 


(231) 


(232) 


/awr      awa\         i  /a(rw«) _  awr>     j_ /a_wk  _ afrv/uV 

la*  dr  )  r\$r  B9  J        r  Ue       52     J 


KB^-SP] 


Wu-F 


awa 


ae 


MMJl) 


+ 


ir{»tf[|^-|^]-«.[^-|^]} 


(233) 


WX  2  or  =(LelA/u+iJWa+irWr)x(UZar)  =  ir(2ar|/Vu)-Le(2arWr) 


(234) 


TVs  =  T(Uf  ff-  +  UH-  +  irif-) 


(235) 
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Combining  (231)  through  (235)  and  writing  the  equations  for  each  unit 
vectors 


d*L 


ae 


Wr  La*       ar    J     T~ 


a? 


_d(rlVu)l       \A/r  [a(rlVu)      aA/r]      ,      w  T    ds  z,,^ 

'd2~"J-  "l^Lar  -Zurwr       —go-         (236a) 

(236b) 


Je~     a?    J  "*"  '  "ar 


Lr'    Br     ~    r 


B(rWu)      dWr]       u/ 


a  wr     d\A/0 


Id*      arj 


4-ZoaWu-r-T-§f 


(236c) 


For  assumed  axisymmetric  flow,  all  derivatives  with  respect  to  0  are 
zero.     Then: 


■  Wr  a  ?      r  a  r 


dz 


r   32 


az 


:  .  aMR  _  vvu  a(riVu)    tt/  awr  ,  u,  awa  ,  ,,  rlA/  ,  T  a^ 


Eq.  (237a)  can  be  rewritten: 
Substituting  this  expression  in  (237b): 


Tr-^aT—^a^+TTi^a?      r2aJ  "wT  +  '  a? 


WuWr 


(237a) 
(237b) 
(237c) 


(238) 


or: 


ang_   wuwr  d(rM*)    w  aWr.u/  a^«    r, ,  ^^    Tas  f2vn 


Multiplying  (237c)   by  Wr  and  (239)  by  Was 


<9H*_  Wal/Vr  c9(rWu) 


aw* 


^^=^^-^^,§^+^^1^  +  20,^^^^^ 


v#= 


^^^^^l^.^^a^,^^^^^^^ 
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Adding  these  two  equations; 

Wrf^  +  K  ft  =  T"(wr|f  +  Wa|f)  (240) 

For  adiabatic  flow,  total  relative  enthalpy  is  constant  along  any  given 
streamline.  Hence,  the  stream  derivative  is  zero.  With  B/dd   =  0: 

Wd±-VHR=0  =  Waj^+\A/r|^-  (241) 


or: 


jyk=_^aH*  (241a) 

Then,  with  (241),  (240)  becomes? 

.3s. _  _  y^il  (242) 

Substituting  (241a)  and  (242)  in  (239) s 

Wa^r  "  rWA   dr        "rvvr5z   VVr  dr   *  *"  Wa    Wa  '  3r 
Then,  multiplying  this  expression  by  (-Wa/Wr)  gives : 


-3T-7ar Kj^  +  WA^-+  Z*K+  \-jr  (243) 


Comparison  of  (243)  and  (237c)  shows  the  expressions  to  be  identical. 
Hence,  this  expression  is  the  equation  which  must  be  solved.  Rewriting 
(243)  and  with  (Wa§^  =  if^)  there  iss 

^.2^  +  ^|^+4caW1|.£|&+2r#=0  (244) 

Eq»   (4),  p.  40,   can  be  written? 

HR=HE-4"  (4) 
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Also,  (5),  p.  40,  can  be  written: 

,2 


T*  = 


He 


W, 


Eqs.  (4)  and  (5  )  are  now  introduced  into  (244).  The  nature  of  the 
equations  relating  total  relative  and  equivalent  enthalpies  dictates 
that  (244)  be  applied  to  the  rotor  discharge  with  this  substitution. 
Then: 


(5f) 


-I-  2 


2c 


(245) 


The  multiplier  of  the  final  term  can  be  expressed: 

ir  [z HE  -  «+  <+<)]=  £  [2H,-<-<  (l  +  $■] 
Eq.  (246)  becomes,  with  (tan2  A  =  W2/W2)  and  (1+tan2  A  =  l/cos2*   ): 


(246) 


2H, 


-^]  = 


zUr-w: 


Cp  C052Ai 


(247) 


Eq.  (245),  combined  with  (247),  and  with  some  rearrangement,  gives: 


dr2 


a>r2         art  Cp 


2HE-WU^  =  0 


There  is  also; 


a(u22)     d(<»7<?) 


dr, 


dr, 


Brz 


=  2uAif£  =  ^r, 


Eq.   (248)  and  (249)  combine  to  give: 


(248) 


(249) 
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-f-Z^r.+  l^H.-W^ 


=  0 


(250) 


To  make  this  equation  non-dimensional,   it  is  multiplied  by  (r^/W^  ) 
where  the  second  subscript,  m,  denotes  the  mean  streamline  value.     Then: 

r2rn  atWg'J        H*  %»      afrfe)     \fitlrM ds>  ,  ?lMUwa&J&Sfc) 


44 


taSm  Wu2  % 


K. 


Wa2     V^^I  ^  cr 


I?     0        (251) 


3fa 


Here,  it  will  be  observed,  the  subscripts,  f2f,  can  not  be  neglected  un- 
less the  equation  is  to  be  applied  to  that  station  (rotor  exit).  Since 
the  station  of  interest  is  station  two,  the  subscript  will  be  dropped. 
Let: 


Wa 


W, 


a.m 


and: 


U 

x= 

Eq.  (251)  can  now  be  written  in  terms  of  these  non-dimensional  quan- 
tities: 


(252) 
(253) 

(254) 


d(YZ)       ?  t  r   dWr       V2       ds*  ,   ?y  tan/f  diXYtan/?) 


A  LLYtan/?      ?   ^Mg  ,  9  Um  U, 


<ax     •  % 


Om 


2M* 


-ftan2/? 


ds< 


dX 


=  0 


(255) 
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But: 


oy  tan/? 


fy(XYtan/?)=2V^ 


ax 


XY 


i^  +  Xton^+Y^f* 


=  ZY2tan/?|^  +  ZYtan^^4-^tan> 


ax 


Also,   there  is: 


and: 


d  tan/?  =  J M 

3X  "  C05^/P  dX 


4-     ~2z>  ^Y         4-2/,  c)(Y    ) 


Eq.    (255)   then  is 


T—  (l-Htan/j-^TTT-r^^--— r^   +2Y    777^  +  7- *<V 


ax 


Wa  "d*    "Vos2a8X 


■  4UnYtqn/<?      2    6>HE   .  7LUUl  1 

T  ~ TTTI  T7T3V  t  £    ,7/1  "  "» 


w, 


atn 


tiLM 


'am 


IN 


am 


*&-**«» 


ds*_ 


ax 


=   0 


Multiplying  by  (cos2/?  /Y2),  and,   since  (1+tan2/?     =  1/cos2^      ): 

1      dill)  -i-m^J     l^md%       I        35Vn    .^iliZciAj. 


(256) 


f£  =  0  (257) 


The  second  term  of  this  equation  accounts  for  the  effects  of  streamline 

curvature,  as  represented  by  the  cos2  A  and  (-  rp£   £^r  ).  Fig.  15, 

wa  &  * 

p.  42,  shows  the  sign  convention  for  streamline  curvature  (i.e.,  positive 
for  the  radius  of  curvature  in  the  positive  coordinate  direction).  Fig. 
15  defines  also  the  streamline  displacement,  $T,   and  the  length  over 
which  this  displacement  occurs,  L.  Streamline  curvature  can  be  approxi- 
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mated: 


17 


COSzX- 


Lz  +(8r)2 


(258) 
(259) 


It  must  be  noted  that  (258)  and  (259)  are  defined  here  for  curvature 
across  the  rotor  only,  because  of  the  small  blade  heighth  of  the  stator 
and  because  of  the  large  flare  area  available  across  the  rotor.  In 
(258),  k  is  positive  for  positive  curvature.  Putting  (258)  and  (259)  in 
(257): 


Mfi=-^ 


:Urmf)-{^%]-2tans&-f&- 


2  Mr  cos  V     e-mr 


M^--™? 


"flrr. 


te  (26°> 


^amY  IACYS 

It  is  evident  that  (260)  is  non-dimensional,  except  for  the  terms  contain- 
ing equivalent  enthalpy  in  BTU/lbm;  lets 

C,  =  2  <j  I  (ft2-lbm/BTU-sec2)  (261) 


Then,  (260)  becomes  completely  non-dimensional  as: 

l(±^_(i^)^_2tan/^_|.5in>_ 


dx 


4U...5ia^co5/?  _  2LUUcosfr?   ■  l.cosk?  dH£    [c,  HE  cos^     ,:njJ_ds* 
«fc-Y  CY1  ^Y2  dX     [vC  Y1  H"^ 


(262) 


Eq.  (262),  the  final  result,  is  included  in  the  basic  paper  as  (26). 


-"-^Vavra,  M.  H.  Aero  Thermodynamics  and  Flow  in  Turbomachines 
(New  York,  London:  John  Wiley  and  Sons,  Inc.,  1960)  p.  453. 
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B.  Theoretical  Basis  for  the  Area  Restriction  Factor,  s  . 

Fig.  50  shows  typical  conditions 
at  the  exit  of  a  blade  row,  with  boun- 
dary layer  on  both  sides  of  the  flow 
channel.  Outside  the  boundary  layer 
the  velocity  is  very  nearly  the  the- 
oretical velocity.  Within  the  boun- 
dary layer,  turbulent  flow  can  be  as-  qq  / 
sumed,  so  that  the  velocity  ratio  can 


18 


be  taken  as: 


(263) 


Illustration  of  Boundary- 
Layer  Effects  at  the  Exit 
of  a  Blade  Row 

FIGURE  50 


The  mass  flow  rate  leaving  the  blade  row  can  be  written: 


h  =  /An\/rHcos^(s-|5?^-f^)+  2[  <y*y 


m 


(264) 


Here,  j^.h   and  Vj^  are  the  values  for  isentropic  expansion  from  inlet  to- 
tal pressure,  pt  ,  to  the  pressure  at  the  discharge,  pj,  which  is  assumed 


constant  across  the  spacing,  s. 
Fig.  51  depicts  the  thermodynamic 
process  for  the  flow  outside  the 
boundary  layer,  and  for  a  point  in- 
side the  boundary  layer.  It  is  evi- 
dent that  the  first  term  of  (264) 
gives  the  flow  rate  outside  the  boun- 


YrlM 


^Tcp 


rtd 


/?» 


\ 


2qJcr 


Thermodynamic  Process 
through  a  Blade  Row 


FIGURE  51 
dary  layer,  while  the  sum  of  integrals  is  the  flow  rate  within  the  boun- 


18 


Vavra,  M.  H.  Problems  of  Fluid  Mechanics  in  Radial  Turbomachines, 


Parts  III  and  IV  (Rhode-Saint-Genese,  Belgium:  Von  Karman  Institute  for 
Fluid  Dynamics,  1965)  VKI  Course  Note  55b,  pp.  G46-50, 
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dary  layers.  The  discharge  angle  is  given  byt 

«d  =  cos"1  (5  °L  t/c0Sfltf) 

With  (265),   (264)  becomes; 


(265) 


■ i/       a  cos**      {scosxf  -t       «£  c/scos^a  -t\    ,    <      (/    a  c  f^  w  j. 


Cancelling  terms,   there  is: 


m=AA„a|l-i.U'-!^U>)) 


For   the  constant  pressure  assumed  at  the  discharges 


/£__  JVjl_    Tto  -  (Tt.  -  TTH)  I  ~Xc 

where : 

*— ft) 

Eq.  (266)  can  then  be  written,  with  (263),  ass 


T-l 

r 


m 


:A> 


VTHa{.-^|[l-(.-Xe)(^^dn]) 


Defining  the  so-called  displacement  thickness s 


s*=s 


'-('-x«HOT*)dn 


Eq.  (269)  becomes s 


™=aX  a 


i- 


£S* 


(266) 


(267) 


(268) 


(269) 


(270) 


(271) 


For  V,,  the  actual  average  discharge  velocity,  the  loss  coefficient 
of  the  blade  row  iss 


<--«$ 


(272) 
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This  loss  coefficient  can  also  be  expressed  in  terms  of  the  average  loss 
of  kinetic  energy.  Hence: 

S  =  -^p-  =  \--^r-  (273) 

where  E  is  the  actual  kinetic  energy  of  the  flow: 

E=A,X  (a  -1  S)4-  +  I  (  W'd*  <274> 

Again,  the  first  term  represents  the  value  for  the  extra-boundary  layer 
flow  and  the  second  term  the  value  for  flow  in  the  boundary  layers, 
Eq.  (274),  transformed  in  the  manner  used  to  obtain  (269),  is: 

&=/v«  \^-&\>-"-™t£h?*M\  <275) 

Then,  the  so-called  energy  thickness  is: 

g~_  *[,-(,_*> j^py*]  (276) 

The  loss  coefficient  can  now  be  expressed,  using  (271),  (275),  and  (276), 
as: 


5  =  i- 


/^H-^a 


FTFI 


or: 


i  -Is- 

1=1-      1%  <277> 

I   ^  a 

The  area  restriction  factor,  f  ,  represents  the  fraction  of  the  flow 
area  which  would  be  usable  for  uniform  theoretical  velocity.  The  defi- 
nition of  f  is  then: 

=  I  -  ^—  (278) 
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Also,  the  energy  parameter  is  defined  as: 


H***=   ±r-  (279) 


Using  (277)  and  (279)? 


(280) 


,*** 


*   H  -I  +t 

Here,  ^  is  the  profile  loss  coefficient. 
C,  Development  of  the  Energy  Parameter, 

The  integrands  of  (270)  and  (276)  can  be  evaluated  by  expanding 
the  denominator  by  means  of  the  binomial  theorem: 

The  integral  of  (276)  is  then: 

)o[7l      +Xe7l   +  Xe?l   +/e7l   +Xen   +....JdTl 
Integrating: 

(      ?*  dn  -         !         -I-        Xe     -u       Xe         ,        Xe  ,       Xe         . 

J0    !-XeTl*m    uri  "      3m+l    ^  5m  +  l     ^7m  +  l  9m  +  l      +   llm-H        •••• 

Eq.  (276)  can  now  be  expanded,  giving: 

S*"_  i     /     I Xe       ,      XeZ      .       X/       .      XeJ\ ,       Xe  Xe*     ,     xl  Xe4     ,     x/ 

g      —   '      ^3m  +  l  5m+l     T  7m+t      +  9m4  I  ilwfll       3m  +  l         5m  +  l         7m  +  l       9mH  T  Hm+i 

or: 


Then: 


Q_  _  I  _  /  y    _  a/      I         .       Xe        .        Xe        .        Xa       ,      Xe     ^ 
C      -I        lAe     "HSm  +  l     ■*"    ?m  +  /  7m+l  9m+l      T  Wm  +  I   / 


5**'_/y        ft/      I  , L      i         Xe  Xl  XJ       ,        Xe     ^  (281) 


~r 
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By  similar  expansion  and  integration,  there  is  for  /^     : 

-f  =(^-|)(w  +  T^ti^r+5Jr-  +  iJr+  d£r]       <282> 

Then,   from  (281)  and  (282),   the  energy  parameter  is: 

I  +  I  ±         Xe        +        Xe         ■         Xe        +        Xe 

«***_         Xe-I  3rvt  +  l  5m+l  7m+l  9m  +  |       '    llrrn-l  ,00^ 

W   ~        i.i       .      Xc     ,.     X8     .      xl     .      xi  (283) 

X«-l  toi+l  3m  +  l  5mfl  7m -H  9m+l 

Fig.  52  shows  H***  as  a  function  of  pressure  ratio. for  f  -   1.4  and  m  -  0.15. 
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APPENDIX  II 
COMPUTATION  OF  FLOW  ANGLES  AND  LOSS  COEFFICIENTS 

A.  General. 

1.  The  basic  method  used  to  estimate  gas  outlet  angles  and  loss 
coefficients  is  based  on  formulas  and  graphs  developed  by  Ainley  and 
Mathieson.  [4,5]   Where  computations  are  based  on  information  from  other 
sources,  footnotes  will  so  indicate.  Fig.  53  shows  definitive  blade 
geometry  pertinent  to  the  calculations. 

B.  Formulas  Used. 

The  General  formula  for  gas  out- 
let angle  prediction  for  Mach  Numbers 
of  0.5  or  less  is? 

<X0=tan"[(l-A)tano<j  +  Atari  <]       (284) 
Where : 


a  _C  A_  COS  0(i 

H        h  cos<xj 


(285) 


BLADE  GEOMETRY  USED  IN 
GAS  ANGLE  PREDICTIONS 

FIGURE  53 


Here;  F  is  a  factor  which  accounts  for  the  type  clearance  (1.35  for  the 
unshrouded  blades  of  the  TTTR  rotor);  h  is  the  blade  heighth  (in.);  CK  . 
is  the  inlet  blade  angle;  and  o(,   is  the  gas  outlet  angle,  uncorrected 
for  tip  clearance,  given  by: 


*  -  •  5 


*d'  =  o(6   +  4f 


(286) 


In  (286),  s  is  the  blade  spacing  and  e  is  the  curvature  of  the  back  of 
the  blade;  These  two  quantities  are  given  by  (see  Fig.  53): 


,*    -if  a 


(287) 
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and: 

e=iif-  (288) 

When  a  blade  row  has  zero  tip  clearance  the  gas  outlet  angle  is  given 
by  (286). 

For  a  Mach  Number  of  one,  the  mean  streamline  gas  outlet  angle  is 
determined  from: 

(«dl=,=  coffjf-)  (289) 

where  Aj.  is  the  throat  (exit  opening)  area  and: 

^d  =  T-(^-DHl)-^^rT-r„);   Dr-  =DT+ZK  (290) 

It  is  evident  that  (289)  is  basically  (287).  Flow  angles  for  other  Mach 
Numbers  (between  0.5  and  1.0)  are  estimated  by  drawing  a  smooth  curve 
between  the  values  for  0.5  and  1.0,  but  with  an  inflection  point  at 
M  =  0.75. 

The  basic  equation  for  the  profile  pressure  loss  coefficient  for 
aero  incidence,  (Yp)^Q,  is: 

w{  w. + eW[w~i-  V.JK-s^J  (291) 

Here;  t  is  the  maximum  blade  thickness,  and  the  ratio,  t/c,  is  restricted 

to  a  maximum  of  0.25;  and  the  equality  subscripts  indicate  which  data 

curves  are  to  be  used.  The  values  of  profile  pressure  loss  coefficient 

for  other  then  zero  incidence  angle  are  obtained  from  curves  of  Y^/fif  ) .  _ 

to  P'  x  p  i=0 

versus  i/ie,  where  i/i.  is  the  ratio  of  incidence  angle  to  positive 
stalling  incidence  angle.  This  latter  value  is  defined  as  the  positive 
incidence  angle  at  which  the  profile  loss  is  two  times  its  minimum  value. 
For  the  calculations,  is  is  determined  from  curves  based  on  a  solidity  of 
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0.75,  with  corrections  applied  for  the  actual  solidity. 

Secondary  and  tip  clearance  losses  are  determined  using: 

COS3- 


^Wx+Bi-M2-^  (»») 


where:  Of  .  the  mean  of  the  gas  inlet  and  outlet  angles,  is  a  function  of 
incidence  angle;  B  is  a  type -clearance  correction  factor  (0.5  for  radial 
clearance);  and: 

|fc-=  2 (tan  o(i  +  tan  «d)  COS  *m  =  f(i)  (293) 

X=f  [(Ad/Ai//(l+  I.D./O.D.)]  (294) 

Here: 

A()=?hnan  (295) 

The  terms,  I.D.  and  O.D.,  are  the  annulus  inner  and  outer  diameters  in 
the  exit  plane,  respectively.  Values  of  X  are  obtained  from  a  curve 

based  on  the  bracket  quantity  of  (294). 

1 Q 
Leaving  losses  are  found  from: 

\,„=  C0Sl«i  ^j^  (296) 

where,  see  Fig.  50,  p.  156  : 

S  _  8  +tc a  "*"  a  (297) 


and: 


«S    a  +  te      ill 

1  ^  a 


a    H^-IH-«I+Skl+M 


(298) 


19yavra,  M.  H.  Aero-Thermodynamics  and  Flow  in  Turbomachines 
(New  York,  London:  John  Wiley  and  Sons,  Inc.,  1960)  p.  89. 
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The  energy  parameter,  H***,  is  obtained  from  Fig.  52,  p.  161. 

The  pressure  loss  coefficients  given  by  (291),  (292)  and  (296)  are 

defined  as: 

Y  m   Loss  in  Total  Pressure 

Outlet  Total  Pressure  -  Outlet  Static  Pressure 

For  use  in  the  performance  prediction  program,  values  obtained  must  be 

converted  to  the  loss  coefficient,   5  .  For  incompressible  flow,  assumed 

for  these  calculations  of  loss  coefficients,  the  required  conversion  is: 

1>  =  777-  (299) 

This  conversion  should  not  be  made,  however,  until  the  total  pressure 
loss  coefficient  has  been  determined,  since: 

I^t)^f^  (300 

Nevertheless,  in  these  calculations,  it  was  necessary  to  convert  to  £  for 
the  sum  of  profile  and  secondary  and  tip  clearance  values  of  Y  in  order 
to  obtain  the  leaving  losses.  No  appreciable  error  results  due  to  the 
relatively  small  magnitude  of  such  leaving  losses. 
C.  Calculations . 

Calculations  for  determining  stator  gas  outlet  angles  are  given  in 
Table  VIII,  while  stator  loss  coefficients  are  computed  below.  Table  IX 
contains  the  computations  for  the  rotor  gas  outlet  angles  and  those  for 
the  rotor  loss  coefficients  are  tabulated  in  Table  X. 

Figs.  54  and  55  show  the  variations  of  predicted  stator  gas  efflux 
angles  with  Mach  Number  and  Radius,  respectively;  Figs.  56  and  57  show 
similar  information  for  predicted  rotor  gas  outlet  angles.  Fig.  58  shows 
rotor  total  loss  and  velocity  coefficients  versus  incidence  angle  ratio, 
i/is. 
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Estimation  of  Stator  Loss  Coefficients 
Losses  will  be  taken  as  constant  with  radius  and  equal  to  those  for 
the  mean  radius. 

Y  =f(s '0,0(1   ,  of,  ): 

s'c=0.^03;  &\=0\    cXf=7*>.2f>° 

Y  =0.C5  (from  curves) 
P 

Y  =A(Ct/(s/c))2(cos2<Vcos3gCh  ) 

3       L 

<*m=tan_1(tan  <*i /2)=63.95° 
CL/(s/c)=2tan^  cosc^  =2(4.09)  (0.439)  =3.  59 
(CL/(s/c))2=12.9 

=f((AH,/A.)2/(l+I.D./0.D.)) 

tn   i 

Ath=aZh=(0.197)(31)(0.69)=4.2U,in.2 
A  «(TT/4)(9.  172-7.792)=18.38in.2 

(A  , /A.)2/ (1+I.D./0.D.)=(4. 216/18. 38)2 (1+7. 79/9. 17)  =0.02855 

tn   1 

=0.0054;   008^=0.05^.7;   cos 3^^ -0.0848 


Y  =(0.0054) (12. 9) (0.0567/0. 0848) =0.0467 
s 

Y  +Y  =0.09CJ 
a  p 

f  +f  =(0. 09^.7/1. 0967)=0. 0883 

mm  X 


*** 


H 


=2.28  (p/p  =0.843);   t  =0. 


024 


S  /a=(  "§  +£  V(H  '  -1+  i  +$   )=(0.0883)/(l. 28-K). 0883)=O. 065 
s   p  s   p 

^;s=(  S*'a+t  /a)/(l  +  t  /a)=(0.065-K).1218)/(1.1218)=0.1663 
>mix=(0.0563)(0.1663)2/(0.8337)2  =  0.00225 

Smix  =  °-°02245 
* total  =  °-09°54 
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TABLE  VIII 
CALCULATION  OP  PREDICTED  STATOR  GAS  OUTLET  ANGLES 
ITEM  HUB  MEAN  TIP 

M  <  0.5 
radius  (r)  -  in.  5.895  4.24  4.585 

radius  ratio,  X  »  r/rm  0.919  1.0  1.08 

chord  (c)  -  in.  1.425  _ — — — 

spaoing  (s)  -  in.  0.7895      0.8594      0.9295 

exit  opening  (a)  «.  in.         0.177       0.197       0.225 

blockage  (t+  )  -  in. 0.071  — — — 

e  =  }2/Qb;   bro.l6;  j-0.76      ' 4.52 

s-tt  0.7185      0.7884      O.8585 

o 

cos"^ | )  -  dog,  75.75  75.50  74.95 

tfj  =  f(cos-l(— i ))  75.75  75.50  74.95 

s  He 
4(s/e)  O.698  0.76  0.822 

o<izo({  4-  4(s/e)  76.45  76.26  75.75 

Ms  1.0 

hh  =  (51)(0.197)(0.69)  «  4.216  in? 

*1  -fr/4)(9.172  -  7.792)  -  (0.071)(51)(0.69)  r  16.86  in? 
*lK  -  1  =  00»"1   <AtlA>  s  cos"1  (4.216/16.86) 
*1>JC.1  "  75.5 

^=^')msi-^')m<0.5s-°^ 

^l'li-1  75.67  75.5  74.99 
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TABLE  IX 
CALCULATION  OF  PREDICTED  ROTOR  GAS  OUTLET  ANGLES 
ITEM  HUB  MEAN  TIP 

radius  (r) >  in.  5.825  4.55  4.875 

radius  ratio,  X  a  r/rm  0.88  1.0  1.12 

chord  (c)  -  in. — 0. 751— — -— — 

spacing  (s)  -  in.  0.4006      0.4555     0.5224 

exit  opening  (a)  -  in.         0.117       0.157       0.166 

blockage  <t+  )  -  in.  — 0.0597— — — — — 

(s-t+  )  0.5409     0.5958     0.4627 


'9 


cos"1(~Ar )  -  deg.  -69.95      -69-7      -69.0 

s-t^ 

A*  s  f(cos"l( * ))  -  deg.   -69.45      -69.2      -68.5 

,  s"tte 

/^2  s/^2   (flinoe  W«  *  0)  -69.45      -69.2       -68.5 

tip  clearance  (k)  -  in.  •— ——.-0.025—————— 

blade  heighth  (h)  -  in.  1.05 — - 

k/h  0.0258 

type-clearance  factor  (P) ————1,55  -————. 


cos/tfj     (^1  r  62°)  -0.469 — — 

cos/2  0»551  0.555  O.566 

tan/^2  -?«67  -2.64  -2.54 

A  =  PCkAJCoos^/cos/^^)  °*°45  0.0425  0.04125 

(1-A)  0.957  0.9575  0.95875 

(l-Ajtan^  -2-555  -2.525  -2.455 

Aitan^j)  0.0809  0.0799  0.0775 

0  =  ((l^Un^^+AiUn^j))        -2.474  -2.445  -2.558 

/^2  *  tan"1^)  -68.0  -67.76  -67.0 
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TABLE  IX  (continued) 
ITEM  HUB  MEAN  TIP 

The  foregoing  calculations  are  for  M  <  0.5* 

&)Sl»    tean  Radius 
At  =  A^(l-kA)   +   DTk 
A^  ■  aZh  =  (0.157)(6O)(l.O5)  s  9.45 
At  «  (9.45)  (0.9762)    4-   (9.75)  (0.025)  «  8.45+0.0766 
At  s  9.196  in? 
*n2  =(Tr/4)((DT4-2k)2.D|)  -  tte2h 

■Crr/4)((9.8)2  -  (7.65)2)  -  (0.0579)(60)(1.05) 
A^  =  29.2  -  5.645  *  25.655  in? 
/*2  s  oos-1   (A^/A^)  s  oos-1  (9.196/25.655) 
/&2  -  deg.  -69.0 

/*2     "  d6«'  -69.24  -69.0  -66.24 
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TABLE  X 
CALCULATION  OP  PREDICTED  ROTOR  LOSS  COEFFICIENTS 

a.  Profile  Loss  Coefficient,  Zero  Incidence*  *  . 

t/c  =  0.296/0.751  s  0.594        use  0.25 
/0X  =  62° 

/^2  =  -69.2 

(A^2)2  =  0.805 
s/c  =  0.4555/0.751  «  0.607 

"  0.897 

Tp(i=0)  S  (0.058  +  0.605(0.1*  "  0.0J8))(0.25/0.2) 

Ip(1.  0)  =  0.1*54 

b.  Stalling  Incidence  (is) 


s  1.05 


/4(s/c*0.75) 
/2  =  .67.76° 

/2(s/c  =  0.75)  r-64'5° 
4*8  s  *e  -  i8(s/c  =  0.75)  *    4'4° 

/V/2(e/c«0.75)s-°-962 

is/flyc_0  yc\  =  9.O0  (value  ie  conjectural  due  to  extra. 
~"    '        polation) 

i.  =  15-4° 
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TABLE  X  (continued) 
c.   Secondary  and  Tip  Clearance  Losses. 

*8  +  Yk  =l>  *  B^/h)l  t  Ol/<»A)32  (cos^/cos^/J 

B  :  0.5 

k/h  =  0.0258 

oos^2  =  0.14J5 

"Ar  f((A2/A1)V(l  +  I.D./O.D.)) 

A2=  Zhgag  r  (6o)(1.05) (0.157)  ■  8.65  in? 

A2=  Zh1a1  s  (60) (1.05) (0.194)  =  8.85  in? 

(l  +  I.D./O.D.)  s  1+7.65/9.60  =  1.76*5 

(<VA1>2/(1  +  XJ>./0J>.)>  =  <8.65/8.85)2(l/1.765)=  0.558 

Xr  0.05  (extrapolated) 

Bk/h  s  0.0119 

X+BkA  s  0.0419 

(A+BkAJcos^g  "  0.00601 
Y8  +  Yk  r  0.00601  (0I/(s/c))2(l/cos^m) 

i/is                                  -5.0         -1.5         -0.8  0.0        0.5  1.0 

!                                   -40.2       -20.1       -10.72  0           6.7  15.4 

tan/^j                              0.4        0.866     1.247  1.88  2.565  2.84 

(tan^  +  taxy^)/2     -1.022    -0.7875  -0.599  -0.2825    0.6  0.1975 

^a                       *         -45.6       -58.5       -50.95  -15.8        5.45  H.2 

oos/fja                              0.7        0.785      0.857  0.962  0.997  0.98 

tan/^j  -  tan/f2          2.845      5.5H      5.692  4.525  5.0IO  5.285 

0L/(s/o)                        5.98        5.20         6.52  8.55  10.0  10.54 

(Oj/ia/o))2                    15.9        27.0        40.0  69.5  100.  107.5 

(l/oos^a)                   2.92         2.07         1.588  1.122  1.008  1.06J 

Yg+  Yk                           0.279      0.558      0.582  0.468  O.606  0.689 
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d. 


TABLE  2 

(continued) 

i/i. 

-5.0 
0.279 

-1.5 
0.558 

-0.8 

0.0 

0.5 
O.606 

1.0 

*s  +  *k 

0.582 

0.468 

0.689 

Total   Loss   Coefficient. 

VYp(i  =  0) 

4.28 

2.1 

1.42 

1.0 

1.22 

2.08 

y 
p 

0.626 

0.5075 

0.208 

0.1464 

0.1787 

0.5045 

VV^ 

0.905 

0.6455 

0.59 

0.6144 

0.7847 

0.9955 

?P+V^k 

0.475 

0.5925 

0.571 

0.58 

0.459 

0.498 

H***-i  (y2sc 

.5) 

1.28 

1.28 

1.28 

1.28 

1.28 

1.28 

H^-l+V*. 

*'k 

1.755 

1.6725 

I.65I 

1.66 

1.719 

1.778 

§Va 

0.271 

0.255 

0.225 

0.229 

0.255 

0.280 

6*A+Va 

o  1824 

__...qo  rr **          

\J  *  A.  Utl 

0.4554 

0.4174 

0.4074 

0.4114 

0.4574 

0.4624 

S/i 

0.5655 

0.555 

0.545 

0.548 

0.570 

0.5915 

(1-  Vs) 

0.6165 

0.647 

0.655 

0.652 

O.650 

0.6085 

(  *A)2 

0.147 

0.1245 

0.119 

0.121 

0.157 

0.155 

(1-    S/8)2 

w 

0.58 

0.419 

0.429 

0.425 

0.596 

0.57 

(  $/s)2/(l- 

S/b)2 

0.587 

0.2975 

0.278 

0.285 

0.546 

0.414 

^mix 

0.0555 

0.0426 

0.0599 

0.0409 

0.0496 

0.0594 

Omix 

0.0526 

0.0408 

0.0584 

0.0595 

0.0475 

0.056 

$  total 

0.5276 

0.4555 

0.4094 

0.4195 

0.4865 

0.554 

-n=l-  tror 

0.4724 

0.5667 

0.5906 

0.5807 

0.5157 

0.446 

v^\fn 

0.688 

0.752 

0.77 

0.765 

0.717 

0.669 

Take  the  profile  loss  coefficient  as  j 
%*  (Yp/(WYk»ftotal 
f  0.5295     0.187       0.1509    0.0906     0.10         0.1529 
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APPENDIX  III 
COMPUTER  PROGRAM  FOR  TURBINE  PERFORMANCE  PREDICTION 
A.  Inputs . 

Listed  below  are  the  inputs  required  for  Program  PERFORM  in  the  order 
of  the  READ  statements. 

1.  Stator  mean  streamline  gas  outlet  angles  (ALFAM — radians)  for 
stator  mean  streamline  exit  Mach  Numbers  of:  0.5,  0.7,  0.75,  0.8,  and  1, 

2.  Same  as  (1.),  but  for  the  rotor  (BETAM). 

3.  Radii  (RS — in.)  assumed  for  streamlines  one  through  five  at  the 
stator  exit  plane. 

4.  Same  as  (3.)  but  for  the  rotor  (RR). 

5.  Rotor  loss  coefficients,  ZETAR,  for  incidence  angle  ratios  ofj 
0.0,  0.5,  1.0,  -0.8,  and  -1.5. 

6.  Rotor  profile  loss  coefficients,  ZETAPR,  for  Incidence  angle 
ratios  of:  0.0,  0.5,  1.0,  -0.8,  and  -1.5. 

7.  Differences  in  stator  gas  outlet  angles  (DALF— radians)  from 
the  mean  streamline  angle,  for  streamlines  one,  three  and  five. 

8.  Same  as  (7.)  but  for  the  rotor  (DBET). 

9.  Stator  blade  exit  opening  (AS — in.)  for  streamlines  one,  three, 
and  five. 

10.  Same  as  (9.)  but  for  the  rotor  (AR), 

11.  The  number  of  sets  of  assumed  conditions  for  which  computations 
are  to  be  made  (NSETS). 

12.  Assumed  conditions  for  sets  one  through  NSETS:  stator  mean 
streamline  exit  Mach  Number  (AMS);  inlet  total  pressure  (PTO— psia); 
inlet  total  temperature  (TTO — °R);  rotor  exit  mean  streamline  Mach  Number 
(AMR);  wheel  speed  (RPM);  and  an  indicator  (NIND)  which  denotes  whether 
assumed  conditions  except  RPM  are  the  same  (NIND  /  0)  or  not  (NIND  *  0). 
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B.  Outputs . 

Sample  output  is  included  at  the  end  of  the  program.  Since  symbols 
used  to  label  output  are  generally  those  used  in  the  analysis,  printouts 
are  self-explanatory  for  the  most  part.  Printouts  requiring  clarification 
of  meaning  are  given  below  in  the  order  of  the  program  output. 

1.  I-(   )  denotes  the  integrand  of  (   ). 

2.  XI  =  f 

3.  PHI  =  $ 

4.  ARAT  »  a/^ 

5.  The  SUMS  1-4,  WSUM  line  gives  the  values  of  flow  rate  integrals 
between  streamlines  (in  order)  one  and  two,  two  and  three,  three  and 
four,  and  four  and  five,  followed  by  the  sum  of  these  values. 

6.  STREAMLINE  FLOW  FRACTIONS,  M  ■  .  Here,  M  denotes  the  stator 

(M  ■  1)  or  the  rotor  (M  -  2).  Numbers  printed  are  the  value  of  X  =  r/rm 
followed  by  the  flow  fraction  passing  between  the  hub  streamline  and  the 
streamline  at  X  for  streamlines  two  three  and  four  in  order. 

7.  W-FRAC  denotes  flow  fraction. 

8.  The  CONSTANT  INTEGRAND  1-5  line  of  the  ROTOR  SOLUTION  gives  the 
values  of  the  integrand  for  the  angle  part  of  Eq.  (33)  plus  the  entropy 
part  ds*/dX2  for  streamlines  one  through  five. 

9.  INT2  is  that  part  of  the  integrand  of  (33)  for  the  remainder 
of  the  entropy  term. 

10.  INT3  is  that  part  of  the  integrand  of  (33)  due  to  the  equivalent 
enthalpy  term. 

11.  INT4  is  that  part  of  the  integrand  of  (33)  due  to  coriolis 
acceleration  and  rotor  peripheral  velocity. 

12.  INT  is  the  total  integrand  of  (33)*'        "<;  iVossm  ,•«.._•...  •  ♦ 

13.  DSDX1-2  denotes  ds|/dX2. 
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Units  of  printed  values  are:  temperatures — °K;   pressures — psia; 
angles — radians;  velocities — ft/sec;  radii — in.;  and  enthalpy  terms— 
BTU/lbm. 
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APPENDIX  IV 
COMPUTER  PROGRAM  FOR  EXHAUSTER  ANALYSIS 

A.  Inputs, 

Listed  below  are  the  inputs  required  for  Program  EJECTOR  in  the 
order  of  the  READ  statements. 

la  The  several  values  of  polytropic  efficiency  assumed  for  the 
flow:  through  the  turbine — ETAT;  through  the  stator — ETAA;  from  the 
turbine  hood  to  the  ejector  nozzle  exit  plane — ETAD;  to  the  ejector 
nozzle  throat — ETAN;  from  the  ejector  nozzle  throat  to  the  nozzle  exit 
plane — ETAI;  in  the  diverging  portion  (subsonic) — ETAC — and  (supersonic) 
— ETAU;  after  a  shock  in  the  diverging  portion — ETAL.  NSETS  is  the 
number  of  inlet  conditions  for  which  the  analysis  is  to  be  done. 

2.  Inlet  conditions  one  through  NSETS;  total  pressure  (PTO — psia); 
and  total  temperature  (TTO — °R). 

B.  Outputs. 

Sample  output  is  included  at  the  end  of  the  program.  Labels  for 
the  printouts  are  the  symbols  from  the  equations  for  the  analysis. 
Units  are:  temperatures — °R;  pressures — psia;  velocities — ft/sec;  and 
flow  rates — lbm/sec. 
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APPENDIX  V 
COMPUTER  PROGRAMS  PERTINENT  TO  FLOW  RATE  MEASUREMENTS 

A.  Program  ORIFICE. 

Inputs  for  this  program  are:  the  upstream  pipe  diameter  (Dl — in.); 
the  orifice  diameter  (D2 — in.);  the  discharge  coefficient,  Koo  t   for  the 
flange  taps  (AKF);  the  discharge  coefficient,  K  ,  for  vena  contracta 
taps  (AKV);  and  the  Reynolds  Number  for  which  the  value  of  K  applies 
(REV). 

Outputs  from  the  program  are  (see  Table  III,  p.  102):  the  diameter 
ratio  (BETA);  the  constant  in  the  expression,  X  =  Cxw/z,  (XC  =  Cx) ;  the 
constant  in  the  expression   ...  for  Y,  (YC  ■  Cy) ;  the  value  of  A  for  the 
flange  taps  (AF)  and  for  the  vena  contracta  taps  (AVC);  and  the  constant 
C  for  the  flange  taps  (CFL)  and  for  the  vena  contracta  taps  (CVC). 

B.  Program  FLOCAL. 

Inputs  for  this  program,  in  the  order  of  the  READ  statements  are: 

1.  The  number  of  data  points  (N);  atmospheric  pressure  (PA — in.Hg); 
reference  manometer  reading  (TARE) ;  and  the  base  supply  total  pressure 
for  the  test  (PUS— psia). 

2.  For  data  points  one  through  N:  temperature  (TN)  and  pressure 
(PN)  upstream  of  the  flow  nozzle;  differential  pressure  across  the  nozzle 
(HWN) ;  temperature  upstream  of  the  orifice  (TOR) ;  pressures  upstream  of 
the  orifice  for  flange  taps  (PFL)  and  for  vena-contracta  taps  (PVC);  and 
the  differential  pressures  across  the  flange  taps  (HWFL)  and  the  vena- 
contracta  taps  (HWVC) •  Temperatures  are  in  F,  pressures  are  manometer 
readings  (in.Hg)  and  differential  pressures  are  in.H^O. 

Sample  output  is  included  at  the  end  of  the  program.  Column  head- 
ings are  self-explanatory. 
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C.  Program  FL0CAL2. 

Inputs  for  this  program  in  the  order  of  the  READ  statements  are: 

1.  The  total  number  of  values  of  nozzle  discharge  coefficient  to 
be  read  (N) ;  the  number  of  values  of  discharge  coefficient  to  be  read 
for  test  number  one  (Kl)  and  for  test  two  (K2). 

2.  Values  of  discharge  coefficient  one  through  N  for  flange  taps 
(AK)  and  for  vena-con tracta  taps  (AK2). 

Sample  output  is  included  at  the  end  of  the  program.  The  BEGIN  and 
END  values  in  the  output  specify  which  of  the  N  values  of  discharge 
coefficient  were  used  in  the  computations  for  the  DATA  GROUP.  The 
column  labeled  POINT  gives  the  number  of  the  program's  calculation 
subsequent  to  initial  passage  of  the  minimum  error  value  of  discharge 
coefficient.  The  desired  value  of  discharge  coefficient  is  obtained 
from  the  calculation  which  shows  minimum  error. 

D.  Program  LABRINT. 

No  inputs  are  required  for  this  program.  Outputs  are  values  of 
overall  seal  pressure  ratio  (BETA)  and  the  corresponding  values  of  the 
flow  function  (PHI  =  (p  ) . 

E.  Program  LABLEAK. 

Inputs  for  this  program  are: 

1.  The  number  of  data  points  (N);  ambient  pressure  (PA — in.Hg); 
and  the  number  of  throttlings  of  the  labyrinth  seal  (Z) • 

2.  Test  data  for  points  one  through  N:  reference  manometer  reading 
(TARE);  temperature  upstream  of  the  square-edged  orifice  (T);  pressures 
upstream  of  the  orifice  for  flange  taps  (PFL)  and  the  vena  contracta 
taps  (PVC);  differential  pressures  across  the  orifice  for  flange  taps 
(HWFL)  and  vena-contracta  taps  (HWVC);  turbine  plenum  temperature  (TP); 
and  turbine  plenum  total  pressure  (PPT).  Temperatures  are  in  F, 
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pressures  are  manometer  readings  (in.Hg)  and  differential  pressures 
are  in  in.^O. 

Sample  output  is  included  at  the  end  of  the  program.  Column  head- 
ings are  self-explanatory. 
F.  Program  MINSQ2. 

Inputs  for  this  program  are: 

1.  The  number  of  sets  of  values  of  pressure  and  flow  function  to 
be  read  (N);  and  initial  assumed  values  for  the  constants  A  and  B. 

2.  Values  of  pressure  ratio  (R)  and  flow  function  (W)  for  sets  one 
through  N. 

Sample  output  is  included  at  the  end  of  the  program.  Values  of  B 
and  associated  errors  during  final  iterations  are  included  in  the  output. 
Values  of  the  flow  function  (W)  as  computed  using  the  constant  B,  along 
with  the  experimental  values,  are  printed  for  each  test  pressure  ratio 
(R).  In  addition,  values  of  the  flow  function  computed  from  the  para- 
bolic equation  are  given  for  pressure  ratios  from  0.15  to  0.95. 
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APPENDIX  VI 
COMPUTER  PROGRAM  FOR  TURBINE  TEST  DATA  REDUCTION 

A.  Inputs . 

Listed  below  are  the  inputs  for  Program  TTTRPGS  in  the  order  of  the 
READ  statements. 

la  The  number  of  data  sets  (NSETS). 

2.  For  each  data  set,  the  number  of  data  points  (NPTS),  the  run 
number  (NRUN),  the  date  of  the  run  (NMOM)  and  (NDAY),  ambient  pressure 
(PA — in.  Hg  abs.),  and  ambient  temperature  (TA —  R) . 

3.  For  each  data  point,  one  to  NPTS,  values  of:  the  reference 
manometer  reading  (PTARE);  the  pressure  upstream  of  the  flow  nozzle  (PN); 
plenum  pressure  (PP);  turbine  inlet  total  pressure  (PTO)  and  static 
pressure  (PO);  the  stator  exit  plane  pressures  at  the  tip  (PIT),  the 

hub  (PlH),  and  in  the  vicinity  of  the  cover  plate  (P1I);  and  the  turbine 
discharge  pressure  (P2).  These  pressures  are  manometer  readings  (in.  Hg) . 
If  P1I  has  not  been  measured,  the  value  of  PlH  is  used  for  P1I  also. 

4.  For  each  data  point,  one  to  NPTS,  values  of:  flow  nozzle 
differential  pressure  (HW — in.HoO);  the  temperatures  upstream  of  the 

flow  nozzle  (TN)  and  in  the  plenum  (TP);  the  inlet  total  temperature  (TTO); 
the  axial  tare  force  (AXTARE — lb);  the  axial  force  reading  (AX)  and  the 
stator  moment  reading  (AMOMS),  both  in  kc.;  rotor  moment  (AMOMR — in, lb.); 
and  wheel  speed  (REV)  in  RPM.  Temperatures  are  in  °F. 

B.  Outputs . 

Sample  output  from  TTTRPGS  is  included  at  the  end  of  the  program. 
Column  headings  are  the  standard  symbols  for  the  turbine  stations  and 
are  self  explanatory.  Temperatures  are  in  R,  pressures  are  in  psia.; 
angles  are  in  degrees,  and  velocities  are  in  ft/sec.  The  column  labeled 
SPEED  gives  rotor  mean  streamline  peripheral  velocity,  U-, . 
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